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PROJECT  SUMMARY 


Phis  report  describes  the  work  performed  at  Physical  Sciences  Inc.  (PS1)  as  a  Phase  II  Small  Business 
Innovative  Research  (SB1R)  program  sponsored  by  the  Office  of  Naval  Research.  The  general  goal  of  this 
program  was  the  continued  development  of  the  flow-  measurement  technique  known  as  the  Vorticity  Optical 
Probe  iVOPf.  The  specific  obiectives  of  the  program  were  to  demonstrate  that:  1)  automated,  time-resolved 
vorticity  vector  measurements  can  be  accurately  performed  at  a  single  well-defined  point  within  a  tlowficld-  2) 
the  VOP  measures  the  correct  vorticity  statistics  in  a  turbulent  flowfield;  and  3)  the  VOP  can  be  made  to  work 
in  water  with  a  new  type  of  water-compatible  probe  particle. 

fhese  three  goals  have  been  met  during  this  program.  An  engineering  prototype  VOP  electro-optical  system 
and  a  data  reduction  and  analysis  algorithm  were  developed  and  tested.  This  VOP  system  was  employed  in  the 
study  ot  several  precision  Coueue  Mow  regimes,  including  laminar,  turbulent  and  periodic  Mow.  The  measured 
vorticity  statistics  compared  favorably  with  expected  values.  Finally,  a  new  type  of  vorticity  probe  particle, 
made  with  acrylamide  gel.  was  invented  in  this  program.  This  invention  permits  the  VOP  technique  to  be  used 
with  water  as  the  working  Muid.  In  many  respects,  the  new  water-compatible  particles  are  superior  probes  to 
those  used  previously  in  non-aqueous  working  fluids.  The  use  of  these  particles  for  vorticity  vector 
measurements  in  laminar,  periodic,  and  turbulent  Couette  flow  has  been  demonstrated. 

The  remainder  of  this  report  details  these  accomplishments  in  six  chapters:  Chapter  1  is  a  general  introduction 
to  the  VOP.  Chapter  2  describes  the  fundamentals  of  vorticity  measurement.  Chapter  3  is  a  detailed  account  of 
the  full  VOP  system  in  its  current  form.  Chapter  4  descnbes  water-compatible  probe  particle  properties  and 
manufacture.  Chapter  5  ts  the  presentation  and  interpretation  of  the  new  precision  Couette  flow  vorticity  vector 
data.  Lastly,  Chapter  6  summarizes  the  progress  on  the  VOP  to  date  and  examines  possible  future  directions  in 
VOP  application  and  development. 
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1.  INTRODUCTION 


'The  imp  Ttance  of  vorticity  for  the  understanding  and  description  of  the 
turbulent  flow  or  high  Reynolds  numoer  motion  of  a  uniform,  incompressible 
fluid  cannot  be  over-emphasized. " 

-P.G.  Saffrm 

The  observation  of  Saffman  quoted  above,  and  similar  sentiments  from  countless  other  investigators  of  t 
flow  phenomena,  have  resounded  through  the  fluid  dynamics  community  for  many  decades.  It  is  only  v 
the  past  decade,  however,  that  a  systematic  exploration  of  the  vorticity  field  of  turbulent  flow  has  beco: 
possible.  This  recent  progress  is  the  result  of  two  significant  developments.  First,  the  explosive  growti 
the  past  decade  of  high-speed  computers  (approaching  10^  operations/ sec)  capable  of  resolving  features 
fundamental  interest  in  turbulence  theory,  and  second,  new  experimental  methods  for  vorticity  measurer 
which  have  been  recently  devised,  of  which  the  vorticity  optical  probe  is  one.  (For  reviews  of  these  mt 
see  Wallace.  1986  or  Foss  and  Wallace,  1989.) 

The  interest  in  vorticity  measurement  reflects  the  growing  recognition  in  the  fluid  dynamics  community 
dominance  of  coherent  structure  in  turbulent  shear  flow.  Hussain  (1986)  has  dubbed  these  'coherent  vc 
to  emphasize  the  rotational  nature  of  these  objects,  and  to  distinguish  them  from  irrelevant  large  scale  n 
the  background  flow.  This  was  also  noted  specifically  by  Laufer  (1975).  in  his  discussion  of  the  kinem. 
what  he  calls  the  "quasi -ordered'  structure  of  turbulence;  ”...  vorticity  measurements  suggest  themselve 
most  promising  method  for  the  quantitative  study  of  the  ordered  motion.  Unfortunately,  there  are  sever 
senous  difficulties  with  this  method:  I)  direct  measurement  of  vorticity  has  not  yet  been  successfully 
accomplished  with  sufficient  accuracy;  2)  single-point  measurements  are  inadequate;  and  3)  following  tf 
t'orticity-contawing  fluid  in  a  Lagrangian  sense  is  intrinsically  difficult*. 

The  first  difficulty  has  been  substantially  overcome  by  the  vorticity  optical  probe  (VOP,  Frish  and  Web 
1981).  The  three  components  of  the  vorticity  vector  are  measured,  any  two  of  them  quite  accurately. 
Multipoint  measurement  arrangements  have  been  devised  (Ferguson  and  Webb,  1983).  And  lastly,  the 
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optical  probe  has  'quasi-Lagrangian'  character  in  that  the  vorticitv  is  measured  at  the  location  ot  probe 
particles,  as  they  convect  passively  through  the  sampled  region. 

The  vorticitv  vector,  defined  in  Subsection  1.1  below.  It  describes  the  rotation  of  a  material  point  in  the 
tlewfield.  The  methods  ot  vorticitv  measurement  in  the  flow  interior  are  divided  into  two  categories;  those 
which  estimate  vorticitv  or  vorticitv  components  from  finite-difference  approximations  of  the  local  velocity  field 
derivatives,  and  direct  measurements  (i.e..  not  implicitly  dependent  upon  velocity  measurements)  of  the  local 
rotation  rate  of  fluid  elements. 

In  the  first  category,  finite-ditference  based  methods  place  a  restrictive  upper  limit  upon  the  accessible  Reynolds 
number  of  the  flow  in  the  laboratory  like  their  computational  counterparts  on  grids  of  finite  resolution. 
Alternately  expressed,  the  dimension  of  the  finite-difference  element  sets  a  limit  of  spatial  resolution.  Free 
-.hear  tlows  and.  generally,  non-parallel  flows  place  even  greater  burdens  upon  the  capabilities  of  finite- 
difference  techniques  due  to  the  questionable  nature  of  Taylor's  hypothesis  applied  to  the  rapidly  changing 
length  scales  and  velocities  in  these  flowfields. 

The  second  category,  direct  measurement  of  vorticity,  to  our  knowledge  has  only  two  members:  a  recently 
developed  forced  Rayleigh  scattering  technique  (Agui  Garcia  and  Hesselink.1990),  and  the  Vorticity  Optical 
Probe  (VOP).  The  VOP  technique  was  developed  by  Frish  and  Webb  (1981)  specifically  to  address  the 
problem  of  resolving  the  smallest  scales  of  intense  turbulence,  even  tn  the  presence  of  powerful  mean  shear.  In 
its  present  form,  the  VOP  measures  three  components  of  vorticity  in  virtually  any  geometry,  in  liquid  flows 
including  water. 

1 . 1  Motivation 

Why  measure  vorticity?  The  vorticity  field,  u(X,t),  which  has  been  called  the  "sinews  and  muscles  of  fluid 
motion"  (Kuchemann.  1965),  is  defined  as  the  curl  of  the  velocity  field  UX.t); 
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u i  =  Vxu 

Taking  the  curl  of  the  vorticity  field,  we  have 

Vxoi  =  Vx(Vxu)  =  V(V-u)  -  V‘u 
For  incompressible  flow  with  V-u  =  0 


V2"u  =  Vxoi  . 


The  solution  of  this  equation  for  u  is  given  by 


if X)  =  _!_  (  —  d3X'  ♦ 

4t  J  f3 

where  r  =  >  X  -  X'  | .  the  integration  is  performed  over  the  domain  of  non-zero  vorticity,  and  <b  is  the  t 
of  an  arbitrary  irrotational  flowfield. 


The  first  obvious  consequence  this  solution  is  that  the  velocity  field  is  specified  by  the  vorticity  field  (  tc 
the  gradient  of  scalar  potential  determined  by  the  boundary  conditions).  The  second  consequence  is  that 
most  important  action  in  any  flowfield  takes  places  where  the  vorticity  resides,  and  important  features  o 
flow  can  often  be  described  by  vorticity  measurement  over  relatively  small  regions  of  the  flowfield. 


Vorticity  is  transported  by  the  flow.  For  inviscid  flow,  the  circulation  about  any  open  material  surface  i 
flow,  which  is  equal  to  the  flux  of  vorticity  through  that  surface,  is  a  conserved  quantity  (Batchelor,  196' 


_d 

dt 


=  0  . 


In  viscous  flow,  the  Navier-Stokes  equations  are  written  (in  the  Eulenan  frame) 


Duij 

7jt" 


yV2u)j 
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i  with  the  usual  summation  convention  on  repeated  indices)  where  ejj  is  the  rate-of-strain  tensor.  Because  oj 
land  e.j)  are  invariant  under  Galilean  transformation,  the  Lagrangian  vorticity  equations  are 


duj 

nr 


2 

e;j  +  »>V  wj 


and  the  intrinsic  rate-of-change  of  the  vorticity  vector  at  a  material  point  is  seen  to  be  due  only  to  the  action  of 
\  iscous  diffusion  or  the  stretching  due  to  e-. 


Finally,  the  enstrophy  per  unit  volume,  \J^\  (the  rotational  analog  of  energy  |u|\  is  directly  proportional  to 
the  local  energy  dissipation  rate  per  unit  mass,  t 

£  *  |  01 1  ^  . 

From  the  discussion  above,  it  is  seen  that  many  of  the  quantities  that  characterize  the  dynamics  of  turbulent  flow 
are  intimately  bound  together  m  the  vorticity  field,  whereas  the  connections  to  the  velocity  field  statistics  are 
often  more  subtle,  and  inherently  less  efficient  descriptors  of  turbulence.  It  is  not  unreasonable  to  expect  that 
vorticity  correlations  will  provide  much  more  powerful  descriptions  of  the  coherent  structures  and  vortex 
interactions  which  are  at  the  leading  edge  of  turbulence  research. 


An  example  of  the  efficiency  of  the  vorticity  field  description  is  seen  in  Figure  1.  Figure  1(a)  shows  velocity 
measurements  in  a  plane  of  a  single-sided  wake  producing  a  2-D  map  of  an  apparently  complex  flow  pattern 
(Perry  and  Chong,  1987).  A  number  of  singularities  are  in  evidence,  including  foci,  nodes  and  saddles 
(F.N.S).  A  full  description  of  the  flowfield  would  ordinarily  require  many  such  planes.  Figure  1(b)  shows  a 
postulated  ’vorticity  skeleton'  of  the  flow,  and  the  resulting  velocity  field  computed  from  Eq.  (2).  The  vorticity 
field  in  this  case  is  a  kind  of  hydrodynamic  shorthand  describing  the  essential  topology  of  the  flowfield. 


In  addition  to  these  very  fundamental  reasons  for  vorticity  studies,  there  are  many  highly  practical  motivations 
as  well.  Propagating  acoustical  signatures  due  to  turbulence  are  closely  related  to  vorticity  dynamics,  and  are 
associated  with  the  time  rate-of-change  of  the  (low)  pressure  field  in  the  cores  of  intense  vortices  (Kambe, 
1986).  Flows  near  boundary  layer  control  devices  and  structure  junctions  are  often  vortex  dominated  flows. 
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(a)  Typical  Instantaneous  (phase-averaged)  Velocity  Vector  Field  for  a 
Single-Sided  Wake  Pattern  Using  a  Hot  Wire  Probe 


<] — 1  Jet  Wake  1 — 


B  7409 


(b)  (top)  Side  View  of  Vortex  Skeleton  for  Single-Sided  Structure. 

(middle)  Oblique  view  of  a  typical  cell.  K  denotes  a  unit  of  circulation. 

(bottom)  Commputed  velocity  field  using  the  Biot-Savart  law  (as  seen  by  an  observer  moving  with  the  e 


Figure  1 

Representation  of  a  Complex  Flow  Pattern  with  a  Vortex  Skeleton  (from  Perry  and  Chong,  1987) 


The  behavior  of  high  helicity  (u<  ■  u)  Hows  (Moffatt,  1969)  and  swirling  flows,  turbomachinery  and  propulsor 
(lows,  is  strongly  affected  by  vorticitv  concentrations  at  inlets.  Periodic  flows  and  wake  development  can  is 
best  described  in  terms  of  shed  vortex  interactions.  Also,  computations  are  frequently  formulated  in  terms  of 
vorticitv  to  take  full  advantage  of  the  economy  of  description  afforded  by  the  vorticity  field  in  many  flows  (e.g., 
Aref  and  Siggia.  1980). 

1.2  Background 

The  development  of  the  VOP  to  date  is  sketched  briefly  in  the  following  sections. 

1.2.1  VOP  History 

The  VOP  was  invented  by  M  B.  Fnsh  at  Cornell  University  as  part  of  his  doctoral  thesis  research  (1981).  The 
first  direct  optical  vorticity  measurements  were  single  component  measurements  in  laminar  Poiseuille  flow  and 
in  a  transitional  boundary  layer  (Frish  and  Webb,  1981).  This  work  was  continued  at  Cornell  by  Ferguson, 
who  developed  the  multi-component  VOP  sensor  based  upon  a  dual-axis  position-sensing  photodiode  detector. 
Ferguson  performed  the  first  direct  optical  measurements  of  two  vorticity  components  simultaneously  (Ferguson 
and  Webb.  1984).  This  version  of  the  VOP  possessed  the  potential  of  three  component  vorticity  measurement. 
The  development  of  the  VOP  was  continued  by  Frish  and  Ferguson  at  Physical  Sciences  Inc.  (PSI)  commencing 
with  an  ONR-sponsored  Phase  I  SBIR  program  in  1987,  described  below. 

1.2.2  Phase  1  Program 

The  Phase  1  SBIR  program  at  PSI  had  two  basic  objectives.  The  first  objective  was  to  demonstrate  that  the 
VOP  technique  could  be  extended  to  three-component  vorticity  measurement  using  a  single  detector  surface  of 
finite  solid-angle.  The  second  objective  was  to  demonstrate  the  volumetric  potential  of  the  VOP,  by  measuring 
vorticity  vectors  at  several  distinct  locations  within  a  large  sampled  volume  simultaneously.  A  method  was 
devised  to  independently  determine  the  positions  of  the  individual  measurement  points  within  the  volume. 
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\ orticity  vectors  at  several  distinct  locations  within  a  large  sampled  volume  simultaneously.  A  method  wa 
devised  to  independently  determine  the  positions  of  the  individual  measurement  points  within  the  volume. 


At  the  successful  conclusion  of  the  Phase  1  program,  the  Phase  II  program,  which  is  the  subject  of  this1  re: 
was  beeun.  The  goals  of  this  program  are  summarized  in  the  next  section. 

1.3  Phase  II  Goals  and  Accomplishments 

The  specific  objectives  of  the  Phase  II  program  were  to:  1)  develop  an  automated  instrument  to  acquire  t 
resolved  vorticity  vector  measurements  at  a  single  well-defined  point  within  a  fiowfield;  2)  demonstrate  th 
VOP  measures  the  correct  vorticity  statistics  in  a  turbulent  fiowfield:  and  3)  make  the  VOP  work  in  watei 

These  three  goais  have  been  met.  An  engineering  prototype  VOP  electro-optical  system  and  a  data  reduci 
and  analysts  algorithm  were  developed  and  tested.  This  VOP  system  was  employed  in  the  study  of  sever: 
precision  Couette  flow  regimes,  including  laminar,  turbulent  and  periodic  flow.  The  measured  vorticity 
statistics  compared  favorably  with  expected  values  from  several  sources.  Finally,  a  new  type  of  vorticity 
particle  -n»^  with  acrylamide  gel  was  invented,  which  permits  the  VOP  technique  to  be  used  with  water 
working  fluid.  In  many  respects,  the  new  water-compatible  particles  are  superior  probes  to  those  used 
previously  m  non-aqueous  working  fluids.  The  use  of  these  particles  for  vorticity  vector  measurements  ir 
laminar  periodic,  and  turbulent  Couette  flow  has  beer  demonstrated. 
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2.  FUNDAMENTALS  OF  VORTICITY  MEASUREMENT 


2. 1  Overview 

From  the  definition  of  the  vorticity  field  as  the  curl  of  the  velocity  field,  u, 

u)  =  Vxu 

it  is  a  straightforward  task  to  devise  an  experiment,  at  least  in  principle,  in  which  the  local  vorticity  vector  may 
be  estimated  from  several  (Eulenan)  velocity  measurements  at  strategically  chosen  positions.  One  would,  as  a 

simple  example,  measure  the  velocity  vector  (u,v,w)  at  four  positions  ,i.e.,  Xg  =  (*o,y0,zg).  (Xg  +  dx,yg,zg), 

(x0.y0  +  dy,Zg).  “d  (Xg.yg.Zg -“dz).  From  the  spatial  Taylor  expansion  of  the  velocity  field  at  Xg,  it  is  evident 

that  good  estimates  of  the  velocity  denvatives  (and  thus  the  vorticity  vector)  at  position  Xg ,  can  always  be 

obtained,  provided  suitable  restrictions  are  placed  upon  the  size  of  the  increments,  dX .  Mathematically,  these 

restrictions  on  the  size  of  dX  need  only  guarantee  the  smallness  of  second  and  higher  order  terms  which 
contribute  the  error  in  gradient  estimates,  while  minimizing  the  differentiation  error  in  the  gradient  calculation 

which  grows  rapidly  with  decreasing  dX  in  the  presence  of  noise.  In  practice,  sophisticated  finite  differencing 
schemes  are  routinely  employed  for  velocity  gradient  estimation  in  experimental  (or  computational)  applications 
(Wallace,  1986.  and  Foss  and  Wallace.  1989).  Unfortunately,  turbulence  m  real  fluids  introduces  difficulties  at 
the  outset. 


For  a  fluid  of  given  viscosity,  v.  the  dissipation  rate,  t,  fixes  the  (Kolmogorov)  microscale  of  the  turbulence, 
1/4 

.  This  is  the  length  scale  at  which  viscous  dissipation  strongly  damps  local  vorticity  or  velocity 


field  fluctuations  and  therefore  is  expected  to  be  a  suitable  order  of  magnitude  of  the  increment,  dX 
(Wynguard,  1969;  Klewicki  and  Falco,  1990).  Selecting  a  finite  dX  in  an  actual  physical  realization  of  such  an 
experiment  usually  imposes  a  restrictive  upper  limit  on  the  Reynolds  number  accessible  to  such  measurements. 

Larger  Reynolds  numbers  result  in  turbulent  length  scales  in  the  flow  which  are  smaller  then  the  probe 
dimensions,  resulting  in  the  washing  out  of  small-scale  structure  and  possible  aliasing  of  the  vorticity 
measurements.  The  Nyquist  sampling  theorem  applies  equally  to  spatial  resolution  issues. 
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A  more  sublie  problem  arises  when  considering  the  time  dependence  of  the  tlowfield.  The  implicit  assum 
of  the  foregoing  discussion  is  that  the  Lagrangian  time  derivatives  ot  the  velocity  components  are  negligib 
comparison  with  the  convective  part  of  the  total  (matenai)  derivative. 

d  ^  d 

dt  dx 

In  other  words,  the  tlowfield  is  assumed  to  pass  through  the  (Eulenan)  probe  region  unchanged,  as  if  "frc 
Taylor’s  frozen  flow  hypothesis  is  then 

J_  - 

ax  ”  “  uc  at 

'Uggesting  that  streamwise  spatial  derivatives  can  estimated  from  time  derivatives  of  measured  quantities  a 
point  where  U .  is  the  convection  velocity.  Clearly  this  is  true  in  turbulent  flow  only  when  Uc  is  sufficier 
large  and  constant.  Great  care  must  be  exercised  in  measuring  low-speed  and  highly  non-parallel  flows  u 
finite  difference  techniques  which  rely  in  part  upon  Taylor’s  hypothesis  (Piomelli.  Balint,  and  Wallace,  1‘ 
The  application  of  hot-wire  anemometry  or  laser-doppler  techniques  (Lang  and  Dimotakis,  1982)  to  finite 
difference  velocity  gradient  measurements  also  suffer  from  velocity  measurement-related  errors  and  probe 
interference  or  interaction,  particularly  near  flow  boundaries.  Limited  spatial  resolution  and  the  problems 
nature  of  high  vorticity  (i.e.,  near-wall)  regions,  obscures  much  of  the  interesting  vorticity  dynamics. 

The  other  major  class  of  vorticity  measurement  techniques  comprises  those  which  attempt  to  directly  obse 
local  deformation  tensor.  By  whatever  means  this  is  accomplished,  e.g.,  grids  optically  ’painted’  in  phot< 
dyes,  etc.  (Fermigter,  Goitre,  Guyon,  and  Jenffer,  1982;  Qiu  and  Falco,  1987),  these  are  essentially  qua 
flow  visualization  techniques.  Local  deformations  of  the  fluid  markers  are  analyzed  to  yield  the  fluid  ang 
velocity  about  at  least  one  axis.  Particle  imaging  velocimetry  (PIV)  (Utami  and  Ueno,  1987;  Agui  and  Ji 
1987;  Adrian,  1991)  is  member  of  this  family.  These  approaches  still  rely  implicitly  upon  the  calculation 
Lagrangian  velocity  field  point  by  point,  and  subsequently,  its  gradients.  These  methods  gain  vorticity  pr 
the  expense  of  spatial  resolution  up  to  a  limit  set  by  the  length  microscale  of  the  flow  and  incur  large  em 
that  for  the  reasons  cited.  The  Eulenan  probes  described  above  make  similar  compromises  depending  on 
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The  VOP  is  a  direct  approach  to  vorticity  vector  measurement.  It  provides  fixed  spatial  resolution  limited  only 
by  particle  size,  but  trades  vorticity  precision  against  temporal  resolution.  No  velocity  measurements  are 
implicit  in  the  method.  The  rotation  of  a  sub-Kolmogorov  scale  fluid  element  containing  a  small  particle  is 
probed  directly  as  it  moves  along  a  streamline.  At  very  high  Reynolds  number,  the  VOP  does  not  alias  the 
data.  Behaving  rather  like  low  pass  vorticity  filter,  the  small  scales  average  out.  The  precision  of  the  vorticity 
measurement  in  the  fluid  element  is  determined  roughly  by  the  length  of  time  for  which  it  can  be  observed. 
This  tune  interval  must  not  exceed  the  Kolmogorov  micro-scale  of  time,  t  =  (We)1''^. 


2.2  The  Vorticity  Optical  Probe  Concept 

The  VOP  is  an  elegantly  simple  concept.  A  small  (20  to  SO  pm)  spherical  particle  suspended  in  a  viscous  flow 
will  quickly  spin  up  under  the  action  of  viscous  stresses  at  the  particle's  surface.  The  torque  exerted  on  a 

sphere  is  proportional  to  difference  between  the  sphere's  angular  velocity  vector,  0,  and  the  angular  velocity  of 
the  surrounding  fluid.  uj/2  (for  local  vorticity,  w).  The  sphere's  final  angular  velocity  (when  the  torque  is 
zero)  is  then 


which  is  approached  exponentially  with  a  relaxation  time  of  r. 


r  =  J_  l!  fP 

r  15  v  Pf 


(4) 


(5) 


where  a  is  the  particle  radius,  and  p p  and  Pf  are  the  densities  of  particle  and  fluid,  respectively.  The  small 
masses  and  moments  of  inertia  of  the  spheres  developed  for  this  purpose  insure  that  they  react  to  the  fastest 
fluctuations  of  velocity  and  vorticity  in  turbulent  flow  (Subsection  2.4.2  below). 


The  symmetric  part  of  the  local  deformation  tensor,  the  rate-of-strain  tensor,  ejj  has  no  influence  upon  the 
orientation  of  the  particle  due  to  the  particle’s  spherical  symmetry.  Thus  a  measurement  of  the  local  vorticity 
vector  in  the  fluid  reduces  to  a  measurement  of  the  probe  particle's  angular  velocity  vector.  The  method  for 
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angular  velocity  measurement  is  to  reflect  a  beam  of  light  from  a  highly  reflective  planar  facet  on  or  in  th 
particle,  and  observe  the  sweeping  reflection. 

To  perform  these  measurements.  Fnsh  (1981)  invented  a  method  of  encasing  microscopic  plane  mirrors  ( 
carbonate  crystal  platelets)  in  20  to  30  |im  clear  plastic  spheres  having  refractive  indices  that  match  the  w 
fluid  in  which  they  are  dispersed  as  probes.  Each  particle  is  mechanically  spherical,  but  optically  planar 
reflections  and  refractions  are  effectively  eliminated  at  the  particle  surface  by  index  matching.  Only  the 
embedded  mirrort  s)  remain  visible.  A  beam  of  laser  light  reflected  from  a  probe  particle  at  a  known  loca 
the  fluid,  is  recorded  as  a  function  of  time.  From  the  reflection's  angular  displacements,  the  particle's  an 
velocity  vector,  and  thus  the  local  vorticity  of  the  fluid,  is  deduced. 

Following  the  discussion  of  VOP  geometry  in  Subsection  2.3  below,  a  complete  discussion  of  vorticity  pr 
panicles,  including  their  properties  in  various  fluids,  manufacturing  methods,  etc.,  will  be  presented  in 
Subsection  2.4.  The  three  component  vorticity  measurement  system  is  described  in  detail  in  Chapter  3. 

2.3  VQP  Geometry  and  Traiectorv  Analysis 

The  basic  VOP  configuration  is  shown  in  Figure  2.  The  laser  is  incident  from  the  positive  z-direction. 
plane  upon  which  the  passage  of  the  reflected  beam  is  observed  is  referred  to  as  the  detector  plane.  The 
of  the  reflected  beam  emanating  from  the  origin,  where  it  intersects  the  detector  plane,  is  henceforward  r 
to  as  a  'reflection  trajectory',  or  simply  'trajectory'.  The  mirror-normal  unit  vector,  n,  has  the  usual  poi 
coordinates,  d  and  9,  defined  with  respect  to  the  positive  z  and  x  axes,  respectively.  The  angular  coordu 
the  reflected  beam  are  d'  and  9'  ■  For  reasons  that  become  clear  in  the  following  sections,  the  primed  ba> 
of  vorticity  vectors  is  preferred  to  that  defined  by  the  coordinate  axes. 
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incident  Beam  and  Detector  Plane 


Primed  Vorticity  Basis  Vectors 

B  6940 


Figure  2 


VOP  Coordinate  System 


2.3.1  Differential  Trajectory  Equations 

In  time  At,  the  vorticity  vector  carries  the  mirror-normal  vector  from  its  orientation  when  first  detected,  rig,  to 
a  new  orientation,  h] ,  by  the  action  of  the  rotation  matrix,  M(co,At) .  The  complete  mirror-normal  trajectory 
consist  of  m  vectors  in  time  mAt  and  has  the  representation; 

The  reflection  trajectory  is  simply  related  to  the  mirror-normal  trajectory  as  described  below.  This  procedure 
was  adopted  for  computing  the  trajectory  simulations  in  Subsection  2.3.2  below,  explicitly  in  terms  of  the 
rotation  matrix  and  a  starting  vector,  qg.  However,  n(t)  is  not  linear  in  the  vorticity  componets.  This 
trajectory  generator  has  been  used  in  a  non-linear  regression  algorithm  for  trajectory  analysis,  but  a  serious 

disadvantage  of  this  approach  is  that  it  is  quite  computationally  intensive.  We  have  therefore  adopted  a 
linearized  method  utilizing  the  rate-of-change  of  the  mirror-normal  angular  coordinates. 
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fhe  raie-ot-change  ot  the  mirror-normal  vector  is  related  to  the  vortictty  vector  bv; 

dn/dt  =  —xh 

2 


Writing  the  mirror-normal  components  as 


nx  =  sin  6  cos  <p 
ny  =  sin  6  sin  <j 
nz  =  cos  6 

differentiating  Eq.  (6)  with  respect  to  time  and  solving  for  d0/dt  and  d<£/dt  gives 

<16'  1, 

—  =  —  iu)vcos<&  -  uiySino  1 

dt  2'  ’ 

d<j>'  wz  cot#.  ,  . 

—  =  _  -  _(a)xcosd>  -  ciiyStnol  ■ 


But  the  reflection  (primed)  coordinates  are  simply  related  to  the  unpnmed  coordinates  as 

<*'  =  <t> 

9'  =  26 

thus  the  relationship  between  thi  rates -of-change  of  the  reflection  coordinates  and  the  vorticity  vector 


—  =  (uiyCosO'  -  ojxsind>') 


do'  wz  1  ,6’.  ,,  •  ,,, 

— =  —  -  —cot — (ujycosd)  +  wvsm <t>  ) 
dt  2  2  2[  *  y 


Expressed  in  the  primed  vorticity  basis; 


— L  =  A oi'7  *  Boi'v  *  Coi'x 

dt  z  y  x 

^  =  Dw'  *  E<i)'v  «•  Fw'x 

Ht  z  y  X 
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where  the  coefficients  A-F  are 


A  =  B  =  —  cos<p',  C  =  -sin<p' 

& 


D 


6' 

cot — simp 


2^2 


E  = 


1  ♦  cot_sin£' 
2 


2\fi 


c  1  J' 

,  F  =  -  —  cot — COS0 
2  2 


Equations  (4)  through  (7)  are  called  the  differential  trajectory  equations  (DTEs).  A  trajectory  consisting  of 
many  (6':. o',)  pairs  at  fixed  time  intervals  generates  many  such  equations  using  finite  difference  estimates  of  the 
4'  and  o'  derivatives.  The  vorticity  components  are  then  extracted  by  a  simple  linear  least-squares  analysis  (see 
Subsection  3.2.8) 


2.3.2  Traicctorv  Simulations 

What  ought  these  trajectories  look  like,  and  what  assurance  can  be  given  of  uniqueness  of  the  vorticity  vectors 
deduced  from  them?  The  initial  orientation  of  mirror-normal  vectors  for  particles  moving  through  the  sampled 
(laser-illuminated)  volume  is  completely  random.  Thus,  for  a  single  fixed  vorticity  vector,  there  exists  a  family 
of  possible  trajectories  on  the  detector  surface  corresponding  to  all  of  the  possible  starting  positions  of  the 
mitTor-normal  vectors.  Figure  3  shows  several  representative  trajectory  families  corresponding  to  the  vorticity 
values  indicated.  Figure  3(a)  shows  the  case  of  pure  w'x,  while  Figure  3(b)  includes  w' r  Clearly  ui'x  and  w'z 
are  more  or  less  directly  related  to  the  rate  and  angle  of  inclination  of  the  trajectories'  passage.  The  addition  of 
ui'y  in  Figure  3(c),  however,  gives  rise  to  a  subtle  change  in  the  trajectory  curvatures.  This  curvature,  and  thus 

ui'y.  becomes  easier  to  measure  as  trajectories  become  longer.  In  fact,  the  geometrical  parameters,  r ,  the  local 

de¬ 
rate  of  passage  of  the  reflection,  the  local  angle  of  inclination  of  the  trajectory,  and  __ ,  the  local  curvature, 

dr 

form  a  convenient  alternate  basis  set  from  which  the  vorticity  components  may  be  estimated.  The  addition  of 
trajectory  curvature  is  what  allows  three  independent  vorticity  components  to  be  determined  from  the  motion  on 
a  (two-dimensional)  surface. 
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(c)  0JX'  =  Wz'  =  wy' 

Figure  3 

Some  Reflection  Trajectories 


b  isia 


Note  that  traiectones  never  cross  on  the  detector  plane:  for  each  initial  mirror  orientation,  distinct  vortic 
vectors  possess  one  and  only  one  trajectory.  The  converse  is  certainly  true:  a  full  trajectory  through  a  giv 
point  in  the  detector  plane  corresponds  to  a  unique  vorticity  vector.  However,  one  must  be  wary  of  the 
singularity  in  the  VOP  geometry  where  trajectory  uniqueness  breaks  down.  Figure  4  shows  a  case  of  sue 
singularity  on  the  detector  surface.  The  fixed  point  about  which  trajectories  appear  to  orbit  corresponds  t 
mirror-normal  which  happens  to  be  parallel  to  the  vorticity  vector.  In  such  a  case,  no  information  can  be 
obtained  about  the  magnitude  of  the  vorticity  vector. 


This  null  vorticity  vector  is  the  special  case  of  a  vorticity  vector  which  is  parallel  to  every  mirror-normal, 
results  in  an  inherent  ’blind  spot’  in  the  VOP  system.  This  does  not  mean  that  a  vorticity  value  of  zero  i 
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be  measured,  merely  that  it  cannot  be  distinguished  from  the  highly  improbable  situation  of  a  finite  vorticity 
vector  parallel  to  the  mirror-normal.  Care  must  be  taken  in  the  interpretation  these  events. 

At  the  center  of  the  detector  (Z  =  0.  X  =  0)  ,  the  coefficients  of  the  w'y  vorticity  component,  B  and  E,  in  the 
DTEs.  vanish.  In  this  case  accurate  measuruments  of  the  ui'x  and  w'z  components  can  be  made  without 
contamination  from  the  uncertainty  in  the  w'y  value.  It  is  reasonable  to  bundle  all  the  uncertainty  into  w'y  where 
it  naturally  resides  (i.e.,  the  primed  vorticity  basis  is  a  diagonal  representation),  rather  than  mix  it  in  the 
unpnmed  frame  between  wv  and  to2- 

2.4  VOP  Particle  Characteristics 

In  evaluating  the  many  VOP  particle  candidate  materials  tested  in  this  program,  mostly  in  the  pursuit  of  water- 
compatibility,  several  basic  guidelines  of  particle  ’design"  have  been  applied  for  all  particle/fluid  systems.  The 
essential  features  of  any  type  of  vorticity  probe  particles,  be  they  solid,  liquid,  or  visco-elastic  material,  are  as 
follows: 

1)  The  optical  properties  of  the  particles  must  meet  the  minimum  experimental  requirements  for  acquiring 
useful  data. 

2)  The  dynamic  response  of  particles  must  be  such  that  they  follow  the  flow  fluctuations  at  all  length  and 


time  scales  of  interest. 


- 1  Intrinsic  panicle  behavior  must  not  compromise  the  simple  interpretation  of  the  measured  angul. 
'.elocitv  vector. 

4)  Particle  interactions  must  not  unduly  influence  the  turbulent  statistics  which  are  to  be  measured. 

Each  ot  these  is  considered  tn  the  following  sections. 


2.4.1  Optical  Properties 

Particles  of  optical  quality  consistent  with  the  requirements  of  the  VOP  detection  system  to  be  discussed 

Chapter  3  have  the  following  features  in  common: 

1 )  Nearly  diffraction  limited  reflection  divergence  from  embedded  mirrors  which  are  free  from  del 
surface  contamination.  These  mirrors  should  be  of  sufficient  size  to  give  acceptably  small  divei 
(Basic  lead  carbonate,  3PbCOj-2Pb(OH>2,  in  the  form  of  hexagonal  crystal  platelets  of  15  /im 
diameter  and  0.06  4m  thickness  have  been  used  almost  exclusively  to  date.) 

2)  Refractive  index  matching  of  the  working  fluid/particle  system  to  a  degree  which  meets  the  abo 
requirement  and  limits  unwanted  particle  surface  reflections. 

3)  Optical  transparency  of  the  particle  material  minimizes  scattered  background  light. 


In  order  to  maximize  the  accuracy  with  which  reflection  trajectories  are  determined  over  the  limited  soli 
subtended  by  the  detector  surface,  the  smallest  possible  angular  spread  of  the  reflections  is  desired.  A  r 
divergence  of  no  more  than  15  percent  of  the  full  cone  angle  subtended  by  the  detector  is  a  sensible  lim 

an  FI. 5  optical  system,  this  corresponds  to  a  divergence  angle  of  $  6  deg.  With  15  nm  mirrors  approx 

5  deg  of  this  is  already  taken  up  by  diffraction.  The  Fraunhofer  diffraction  pattern  of  a  circular  apertur 


KX) 


2Jj(X) 

X 


X  =  1  D  stn0 
X 


The  first  zero  of  Jj(x)  is  at  X  =  3.83,  and 
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e  =  1.22A 

D 


giving  ~5  deg  full  cone  angle  for  green  light  (the  Argon  ion  laser  line  at  S14.S  nm).  The  5  deg  full  cone  angle 
of  the  Airy  disc  contains  84  percent  of  the  total  reflected  power. 

Mismatching  of  refractive  indices  for  fluid  particle  systems  produces  focusing  effects  at  the  curved  particle 
surfaces.  The  particle  surface  near  the  embedded  mirror  shown  in  Figure  S  acts  like  a  pair  of  thin  lenses. 

For  2  percent  matching.  An  -  0.03  with  D  =  15  nm 

1  .  1  1  ,1 

_  —  An  _  +  —  ,  t  mm 

f  rj  xi  3 

The  farfield  defoc using  angle  D  f  =  2.6  deg.  Therefore,  index  matching  in  the  1  to  2  percent  range  is 
necessary  for  holding  total  spot  angular  size  near  the  6  deg  limit.  This  also  leaves  some  room  for  improvement. 
Doubling  the  mirror  diameter  would  nearly  halve  the  cone  angle  and  quadruple  the  total  power  in  the  reflection. 


Ancillary  benefits  of  good  (  5  1  percent)  index  matching  include  minimized  surface  reflections  intensity  into  the 

f  12 


collection  optics,  though  this  is  of  order 


U2  +  nj 


—  1CT*  near  45  deg  incidence,  not  including  geometrical 


REFLECTED  BEAM  „  " 


8-6929 

Figure  5 

Illustration  of  the  Refractive  Contribution  to  Farfield  Reflection  Divergence 
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Npreading  or  reflections  from  the  highly  curved  surfaces.  Such  reflections  would  be  negligible  even  near 
10  percent  index  matching.  The  more  senous  difficulty  of  this  type,  however,  anses  from  low-angle  surf, 
reflections  and  multiple  low-angle  internal  reflections.  In  conjunction  with  beam  extinction  and  diffractior 
to  the  mirrors,  this  can  rapidly  decrease  optical  depth  and  produce  inhomogeneity  in  the  incident  and  refit 
laser  light  at  high  particle  loading  (~  10^  particies/cm^). 

Some  volume  scattering  from  the  particle  material  is  acceptable.  The  fraction  scattered  into  the  detection 
from  the  few  particles  within  the  sampled  volume  at  any  one  time  will  be  orders  of  magnitude  smaller  tha 
reflection  intensities. 


2.4.2  Dynamic  Response 

VOP  particles  must  not  only  track  the  translational  motions  of  the  surrounding  fluid  element,  as  is  usualh 
concern  m  connection  with  LDV  seed  particles,  but  the  angular  velocity  and  accelerations  of  that  fluid  elc 


as  well.  For  the  translational  part,  the  discussion  of  Buchhave,  George  and  Lumley  (1979)  is  recounted  \ 


some  modifications. 


If  the  Lagrangian  velocity  field  is  represented  by  Fourier  components  of  random  phase,  then  the  dynamic 
response  of  the  particle  to  the  component  at  frequency  /,  v  =  vq  sin(2ir/t  *  o).  the  solution  of 

0 ..  ,  - 
_(x  -  v)  =  x 
m 

is  required,  where  x  -  v  is  the  velocity  of  the  particle  relative  to  the  local  Lagrangian  field  and  where  B 
(Stokes)  drag  factor  and  m  is  the  particle  mass. 


The  solution  gives  the  spectrum  of  the  particle/fluid  relative  velocity 
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'a  here  rv  is  the  particle  s  response  time. 


*(/) 


r-v(2irf)2 
1  *  r~(2irf) 


To  follow  velocity  fluctuations  to  £  10  percent,  requires  rv(2ir/)  ~  10'^.  For  $»(/)«  1,  rv  is  obtained  simply 
from  the  Stokes  drag  (Luraley.  19571 


The  Kolmogorov  microscales  tor  two  tlow  conditions  in  water  are  summarized  in  Table  1. 


Table  1.  Comparison  of  Turbulence  Microscales  and  Particle  Characteristics 


Case  1 

Case  2 

for  u  -  1  ms  water  tlow.  f  -  10  cm 

for  6  ms  over  2  cm 

e  =  u2/t  =  10"\  f  =  0.01  cm“/s 

e  =  u3/f  -  108 

n  -  20  4m 

•q  ~  3  pm 

i  -6  cm/s 

v  —  30  cm/s 

;  -  3  t  10"*s 

T  -  10'3S 

/  -  3  x  103  s'1 

/  ~  10^  s'1 

a2  <■ 

Tv  - - 

'l  -10-4s 

3v  p f 

for  40  pm  diameter  neutrally  buoyant  particles  in  water 

The  Kolmogorov  microscales  are  the  local  scales  that  can  be  formed  from  the  local  dissipation  rate  c  ( *u3/f) 
and  the  kinematic  viscosity  v.  These  are: 
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velocity,  u 
length,  r\ 

time,  r 

frequency,  / 


(tv) 


1/4 


1/2 


1/2 


v  and  i  are  estimates  of  typical  velocity  and  length  sacle  of  the  (outer)  flow.  The  response  time  of  a  40 
diameter  neutrally  buoyant  particle  is  included  in  the  table  for  comparison.  Note  that  for  the  most  mten 
turbulence  (Case  2).  some  small  scale  structure  is  likely  to  be  smoothed  out. 


In  the  rotational  response  calculation,  the  treatment  is  fully  analogous.  Following  Chwang  and  Wu  (19' 
torque  on  a  spherical  particle  rotating  at  angular  velocity  (1  ina  fluid  of  vorticity  w  is 


8x  pf  1*^(0  -  2  ui) 


As  stated  in  Subsection  2. 1  above,  the  steady  state  is  achieved  when  fl 
particle  is 


w 

2' 


The  equation  of  motion  t 


15  2  -1 
pp  a2 


n  -  iw 

2 


=  fl 


Again,  using  Founer  analysis,  a  similar  spectrum  is  obtained  for  particle/ fluid  relative  angular  velocities 
the  same  conditions  stated  in  the  translational  discussion  above.  rr,  the  rotational  relaxation  time,  is  obi 
from  the  rotational  analog  of  the  Stokes  drag. 


15  v  pf 


Note  that  a  particle's  rotational  response  time  is  approximately  five  times  faster  than  its  translational  res 
time. 
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The  roie  of  the  local  rate-ol-strain  tensor  is  more  subtle.  By  symmetry,  it  can  generate  no  couple  on  a  spherical 
particle  which  would  tend  to  rotate  the  particle.  Non-sphencal  particles,  on  the  other  hand,  will  tend  to  align 
their  long  axes  with  the  pnncipal  strain  rate  in  the  flow.  In  a  simple  shear,  the  sum  of  a  straining  field  and  a 
rigid  body  rotation,  the  rotation  rate  of  elongated  particles  will  modulate  in  proportion  to  their  degree  anisotropy 
as  thev  tumble.  (Jeffrey.  1922:  see  Subsection  2.4.4).  If  the  strain  rate  becomes  sufficiently  great  to 
appreciably  deform  non-solid  probe  panicles,  the  effect  of  the  strain  rate  on  the  measured  vorticity  is  not 
negligible.  This  is  considered  in  the  next  section. 

2.4.2  lntnnsic  Particle  Behavior 

Rigid  body  rotation  is  by  definition  the  only  behavior  that  solid  particles  in  a  shear  flow  can  exhibit  under 
normal  conditions.  However,  liquid  and  visco-elastic  droplets  interacting  with  the  surrounding  tlowfield  will,  in 
general,  exhibit  internal  circulation  and/or  deformations  in  addition  to  their  translations  and  rotations.  Such 
deformations  can  be  catastrophic,  breaking  particles  or  droplets  into  two  or  more  parts.  Fortunately,  these 
effects  become  problematic  only  at  very  high  strain  rates. 

For  an  immiscible  fluid,  surface  tension  is  able  to  maintain  sphericity  of  small  (  ~25  jim)  droplets  against  very 

high  shear  rates.  For  a  40  utn  diameter  droplet  with  20  dynes/cm  surface  tension  with  the  surrounding  fluid 

dynamics  pressure  fluctuation  must  approach  —  -  lO'*  dynes/cm  to  deforrn  the  droplet  significantly. 

a 

Apparently  then,  liquid  droplets  containing  VOP  mirrors  will  also  function  as  vorticity  probe  particles.  The 
most  significant  effect  on  the  measured  vorticity  within  the  droplet  results  from  the  internal  circulation  induced 
by  relative  velocity  when  the  droplet  phase  is  not  neutrally  buoyant.  A  detailed  discussion  of  liquid  droplet 
VOP  particle  characteristics  is  beyond  the  scope  of  this  report,  but  there  are  many  interesting  implications  for 
two  phase  flow  studies  using  variations  of  the  VOP  technique. 
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In  the  case  of  visco-elastic  probe  particles,  strain-induced  deformations  can  also  lead  to  spunous  modul. 
the  measured  local  vorticity.  For  high  shear  rates,  S,  the  magnitudes  of  both  the  vorticity  and  rate-of-> 
tensor  will  be  large.  So  long  as  the  particle  remains  spherical,  the  instantaneous  rotation  rate  remains  : 


If  panicles  are  deformed,  the  rate-of-stram  tensor  generates  a  torque.  Figure  6  indicates  the  anttcipatec 
deformations  under  extreme  conditions.  In  a  worst  case  scenano,  e.g.,  the  shear  rate  near  the  leading 
flat  plate  in  a  water  tunnel  with  free-stream  velocity  30  ft/s,  shear  rates  can  exceed  lCF*  s’*. 


For  100  ixta  diameter  particles  in  water  with  {S,  a,  v)  =  {10^  s'*.  5  x  lO’’*  cm,  0.01  cm“/s},  the  Res 
number  of  the  flow  in  the  vicinity  of  the  particle  can  be  estimated  as 

a^S 

R  *  _  -  250  . 

»> 

It  is  expected  that  such  a  flow  will  be  nearly  potential  a  short  distance  from  the  particle  and  that  the  eft 
viscosity  will  be  confined  to  a  boundary  layer  near  its  surface  of  thickness  6. 


b  ~ 


-  —  5xl0-4  cm 


The  maximum  viscous  stress  experienced  at  the  particle  surface  is 


t  tS  pjy s^a“  -  104  dynes/cm 


Figure  6 

Deformation  of  Visco-elastic  Particle  at  High  Shear 
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The  maximum  attainable  dynamic  pressure  forces  relative  to  stagnation  points  in  the  particle  frame  are 

Ap  -  p(sa)~  -  10^  dynes/cm“  . 

For  visco-elastic  material,  such  as  acrylamide  gel  (described  in  Chapter  4)  whose  bulk  modulus  is 
5  10°  dynes/cm",  strains  ot  less  than  10  percent  are  the  worst  case.  This  corresponds  to  only  a  few  percent 

vorticity  modulation.  The  yield  strain  of  this  material  is  30  to  40  percent,  i.e.,  larger  deformations  simply 
rupture  the  gel  network. 

2.4.4  Particle  Interactions 

The  Einstein  relation  for  the  bulk  viscosity  of  a  suspension  of  rigid  spherical  particles,  /x*.  is 

*  f ,  5 

=  ji  1  »  -  a 

where  n  is  the  continuous  phase  viscosity,  and  a  is  the  total  volume  fraction  of  the  dispersed  phase.  This 
expression  is  generally  accepted  to  hold  for  a  S  0.02.  For  particles  of  radius  a.  this  condition  on  a  can  be 
expressed  geometrically  with  an  average  'primitive  cell’  of  side  r. 


In  each  vertex  of  the  cube  is  a  single  octant  of  the  sphere,  thus  each  cell  contains  the  volume  of  a  single  sphere. 

4  a3 

The  volume  fraction  is  a  =  —  x — ,  and  when  a  is  0.02,  r/a  =  6.  If  the  mean  distance  between  particle 

3  r3 

centers  is  less  than  6a.  the  Einstein  formula  is  no  longer  an  adequate  description.  This  expression  is  based  upon 
the  approximation  of  single,  independent  particles  perturbing  an  unbounded  flow.  Higher  order  terms  describe 


the  flow  field  perturbations  near  a  particle  due  to  its  nearest  neighbors  (Batchelor  and  Green,  1972b).  as  in 


2 

a  *  Cia  *  ... 

Special  corrections  may  be  required  for  particular  flows.  For  example,  Lin,  Peery,  and  Schowalter  (1970)  have 
shown  that  bulk  viscosity  of  a  shear  flow  with  solid  spherical  particles  require  a  correction  for  inertial  effects  at 
high  shear  Reynold  number,  R,: 
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Following  Batchelor  ( 1967),  the  far-field  correction  due  to  the  presence  of  a  rigid  spherical  particle  emf 


a  pure  straining  motion. 


For  a  plane  strain,  i.e. 

^  =  fl  t  =  -II 
3x  3y 

the  tar-tield  velocity  perturbation  in  the  plane,  z  =  0.  is 

-  Byy]i! 

where  iflx*  -  By?]  is  the  unperturbed  straining  fields  at  the  point  (x.y).  The  r/a  =  6  criterion  is  seen 
equivalent  to  the  requirement  that  the  average  velocity  perturbation  at  neighboring  particle  not  exceed 

-  I  percent  of  the  background  flow. 


For  40  um  diameter  particles,  the  r/a  =  6  condition  gives  a  particle  number  density,  a  S  6* 

10-*  panic  lesion^  and  an  increase  in  bulk  viscosity  of  5  percent.  Particle  centers  have  >  120  /im  avera 
separation;  this  degree  of  VOP  particle  loading  is  never  likely  to  be  exceeded.  We  thus  have  reasonabl 
assurance  that  dilute  suspensions  of  probe  particles  sample  the  background  flowfield  and  not  the  flow  fi 
their  neighbors.  Under  some  conditions  flow  conditions,  however,  panicle  pair  interactions  can  be  imp 

Isolated  strong  panicle  interactions,  which  perturb  the  particles'  rotation  rate,  will  occur  in  shear  flows 
panicles  approach  each  other  closely.  A  panicle  pair  in  a  shear  flow  is  shown  below  in  Figure  7  in  the 
reference  frame  of  particle  B.  Particle  A  approaches  particle  B  with  velocity  «b  where  b  is  an  effectiv 
parameter.  As  b  becomes  small  (<a),  the  particle  interaction  is  pictured  sequentially  in  Figure  8  (e.g. 
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IMPACT  PARAMETER) 


Figure  7 

Particle  A  Approaches  Particle  B  in  the  Uniform  Shear  at 
Impact  Parameter  b  with  Velocity 


S-6939 

Figure  8 


Sequence  of  Pair  Positions  During  Interaction 


2' 


Darabaner  anil  Mason.  1967).  At  closest  approach,  at  or  near  contact,  the  particle  pair  rotate  as  a  rigid 
i.e..  as  a  dumbell.  Rather  than  each  particle  rotating  at  angular  velocity  ui/2,  the  formula  for  the  ngid  F 
rotation  rate  ot  an  ellipsoid  of  revolution  with  a  2: 1  aspect  ratio  due  to  Jeffrey,  applies  approximately  fo 
paired  particles. 


with  bounds. 


1(1  <  0„  <  1(1 

5  P  5 


When  the  particles  touch,  they  quickly  accelerate  through  8/50.  pass  through  0  near  vertical,  and  slow  i 
as  they  approach  alignment  with  the  pnncipal  strain  axis.  Then  they  separate  once  again  and  continue  o 
original  paths.  For  small  b.  the  probability  of  such  an  encounter  is  small  since  the  velocity  of  approach 
off  directly  with  b.  The  mean  time  between  such  interactions  is  then  proportional  to  n/b,  where  n  is  the 
number  density.  For  b  >  d,  the  particle  diameter,  stokes  flow  interactions  drop  off  quickly,  as  (b/d)3. 
the  weakest  attraction  between  particles,  however,  may  result  in  a  meta-stable  pair  which  can  persist  for 
rotation  periods  in  a  laminar  shear.  This  effect  is  observed  experimentally,  and  must  be  invoked  to  exp 
laminar  flow  calibration  data  in  Chapter  4. 


Particles  will  also  interact  with  solid  boundaries.  Some  range  restrictions  (h/a  =  distance  to  wall/radius 
Failing  these,  the  possibility  of  lift  and/or  induction  effects  which  can  induce  body  forces  or  torques  on 
probe  particle  is  introduced.  The  exact  solution  of  the  Stokes  equations  for  the  motion  of  a  neutrally  by 
sphere  in  a  shear  flow  near  a  plane  boundary  has  been  calculated  as  a  function  of  h/a  by  Lin,  Lee.  and  .' 
(1970).  The  calculated  particle  motion  for  local  unperturbed  shear  rate  u0  and  velocity  U0=<ii0h  for  sei 
h/a  values  gives 

Apparently  the  particle  must  come  very  close  to  the  wall  (a  fraction  of  a  micron  for  100  /im  particles)  u 
to  significantly  reduce  the  measured  vorticity.  The  reduced  velocity  of  the  particle  must  also  locally  pei 
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at  h/a  = 


at  h/a  = 


and  at  h/a  = 


10. 


1 


=  0.9997, 


Doh 


1.5. 


Q 

1 

2“o 


0.923, 


U 

o)0h 


1.005,  JL  =  0.436,  — 

1  Woh 


0.9996 


0.921 


=  0.495 


boundary  layer  flow  in  this  case.  Still,  corrections  of  vorticity  statistics  very  near  boundaries  are  not  likely  to 
be  required. 
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3.  THREE  COMPONENT  VORTICITY  MEASUREMENT  SYSTEM 


The  V  OP  apparatus  is  a  hilly  automated  system  which  acquires  reflection  trajectory  data  in  real  time  froi 
-.mall  sampled  volume  at  any  point  in  the  flowfield.  These  data  are  sequences  of  discrete,  random  event 
consisting  of  a  digitized  representation  of  a  reflection  trajectory  traversing  the  detector  surface,  from  wh: 
vorticity  vectors  are  deduced.  The  heart  of  the  three-component  VOP  is  the  dual-axis  position-sensing 
photodiode  detector.  This  is  described  in  Subsection  3.1,  along  with  the  imaging  and  data  acquisition  sy 
A  detailed  account  of  signal  processing  and  data  reduction  procedures  is  provided  in  Subsection  3.2,  foil 
Subsection  3.3  by  a  discussion  of  error  recognition  routines  for  ’perturbed’  or  irregular  trajectories  and/c 
traiectorv  fragments.  Lastly,  a  summary  of  system  characteristics  in  Subsection  3.4  gives  the  dependent 
important  system  performance  measures  on  key  experimental  parameters. 

3. 1  Three  Component  Detection  System 

The  detection  system  of  the  VOP  consists  of  three  features:  the  detector  which  generates  electrical  sigm 
based  on  incident  trajectories;  the  optical  system  which  serves  to  define  the  sample  volume  and  blocks 
extraneous  background  light;  and  the  data  acquisition  system  which  digitizes  and  stores  the  raw  trajector 
for  post  processing  with  a  personal  computer.  Each  of  these  is  discussed  in  turn  below.  The  entire  dati 
processing  program  is  included  as  Appendix  A. 

3.1.1  Position-Sensing  Photodiode 

The  dual  axis  position-sensing  photodiode  (PSD)  employed  in  the  VOP  system,  manufactured  by  United 
Detector  Technology,  Inc.,  generates  analog  signals  from  which  the  instantaneous  position  (X.Z)  and  the 
of  the  reflected  beam  are  derived.  The  PIN  SC-25  PSD  device  is  an  approximately  two  centimeter  squa 
biased  photo-conductive  layer  over  a  uniform  resistivity  substrate.  The  axes  are  defined  by  four  posts  o 
periphery  of  the  active  area,  at  ±Zp  and  at  ±Xp.  Current  drains  through  the  posts  in  proportion  to  tf 
power  of  light  striking  the  surface.  Because  the  substrate  has  uniform  resistivity,  it  behaves  as  a  two 
dimensional  voltage  divider  for  a  light  beam  striking  the  surface  at  a  particular  position.  Thus,  while  th 
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.it  the  currents  draining  through  post-pairs  is  proportional  to  the  power  or  the  incident  light  (X-sum, 

Z-cum  oc  1 1,  the  difference  or  the  currents  is  proportional  to  both  the  position  of  the  centroid  of  the  reflected 
beam  and  its  power,  for  each  axis  (X-difference  at  X-l,  Z-difference  a  Z -I).  The  quotient  of  the  X.Z- 
Jifference  signal  and  the  X.Z-sura  signal  is  therefore  directly  proportional  to  the  X.Z-position  of  the  reflection, 
independent  ot  the  power.  The  quantity  (X-sum  +  Z-sum),  however,  is  proportional  to  the  incident  power 
alone,  independent  of  position. 


The  detector  has  a  nominal  sensitivity  of  0.6  amps/watt  and  an  NEP  figure  of  -  1C*’ 1  watt/y/Hz\  which  vanes 
somewhat  with  wavelength.  The  band  width  of  the  transimpedence  pre-amplifiers  is  approximately  50  kHz  at  a 
gain  ot  106  V  'amp:  thus,  the  total  noise  of  the  detector  is  approximately  1  mV  (at  room  temperature).  A 
reflected  power  of  just  a  few  microwatts  is  sufficient  to  give  a  signal-to-noise  ratio  (SNR)  of  about  100  in  the 
sum  voltage  signals.  The  junction  capacitance  of  the  detector  will  allow  speeds  approaching  1  MHz.  which  can 
be  exploited  in  the  future  with  faster  pre-amplifiers.  The  detector  and  its  sum  and  difference  pre-amps  are 
housed  together  in  a  single  shielded  unit. 

The  reproducibility  of  position  measurements  on  the  detector  is  generally  excellent,  but  is  dependent  upon  the 
signal-to-noise  ratio.  The  physical  position  of  the  centroid,  however,  is  not  linear  in  the  jxisition  deduced  from 
the  output  of  the  detector,  particularly  near  the  edges  of  the  active  area.  A  detailed  d elector  calibration  must  be 
performed  and  implemented  in  the  data  reduction  procedure. 

3.1.2  Optical  Configuration 

A  typical  optical  configuration  is  shown  in  Figure  9.  While  a  number  of  variations  are  possible,  this 
arrangement  used  to  view  the  near  wail  vorticity  structure  of  Taylor-Couette  instabilities,  is  illustrative  of  all  the 
important  features. 

Beginning  from  the  incident  beam  from  within  the  fluid,  the  optical  elements  are  described  in  their  order  along 
the  optical  path.  The  first  optically  active  surface  encountered  in  this  geometry  is  the  outer  flow  boundary. 
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INCIDENT  BEAM 


VOP  Optical  System 


Under  some  conditions,  the  wail  material  should  be  refractive  index  matched  to  the  working  fluid.  This 
desirable  when  working  dose  to  the  wall  at  43  deg,  as  in  the  configuration  of  Figure  9,  since  both  stron 
reflections  and  astigmatism  compromise  the  image  quality.  In  the  case  of  water  as  the  working  fluid,  th 
been  accomplished  by  means  of  a  molded  window  section  of  fluonnated  acrylic  plastic,  which  forms  par 
boundary. 

The  second  optical  element  is  the  coupling  prism.  Together  with  the  index-matched  boundary  window  s- 
they  provide  unrefracted  entry  of  the  incident  beam  and  astigmatism-free  imaging  of  the  sampled  volume 
pnsm  is  custom-made  from  either  Plexiglas  or  acrylamide  gel  (see  Subsection  4.2.2). 
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Next  is  the  plano-convex  output  coupling  lens  —  so-called  because  of  its  function  of  allowing  rays  to  leave  the 
prism  and  emerge  from  the  lens  surface  without  refraction;  the  large  solid-angle  requirements  of  the  three- 
component  VOP  necessitate  this  as  a  means  of  eliminating  both  the  angular  magnification  and  sphencal 
aberration  induced  bv  planar  surfaces.  For  fully  indexed  matched  components,  the  center  of  curvature  of  the 
output  coupler's  convex  surface  coincides  with  the  position  of  the  sampled  volume  in  the  fluid,  which  clearly 
illustrates  the  output  coupler's  function. 

The  imaging  component  consists  of  two  lenses.  The  collimating  (collecting)  first  stage  is  a  Nikon  FI. 4  55  mm 
camera  lens  (inverted).  The  second  stage  is  a  Nikon  F2.8  110  mm  camera  lens  which  provides 
2 X -magnification  of  the  sampled  volume  at  the  image  plane.  Note  that  the  angular  magnification,  however,  for 
rays  leaving  the  sampled  volume,  is  0.5X. 

These  lenses  were  chosen  to  afford  maximum  acceptance  angle  with  minimum  sphencal  aberration.  At  less  than 
F2.  if  is  usually  necessary  to  resort  to  multi-element  optics  to  improve  image  quality  beyond  what  is 
commercially  available  from  precision  achromatic  lenses.  This  was  established  in  the  laboratory  having  found 
aspheric  pairs  at  F1.0  and  achromat  pairs  at  FI. 6  unacceptable  for  imaging.  The  imaging  quality  concerns  are 
made  clear  in  considering  the  next  optical  stage. 

Definition  of  the  sampled  volume  is  accomplished  in  two  dimensions  by  the  incident  beam  waist  and  in  the  third 
dimension  by  the  variable  slit  mask  in  the  image  plane  shown  m  the  figure.  For  a  300  p m  slit  width,  a  point 
spread  function  width  through  the  optical  tram  of  approximately  100  pm,  for  example,  would  result  in  more 
than  50  percent  of  the  imaged  area  being  compromised  by  optical  aberrations.  That  is,  all  or  part  of  the  light  is 
clipped  by  the  slit  edges  m  some  complex  way  as  a  function  of  exit  angle.  Such  difficulties  are  best  avoided 
with  high  quality  imaging  available  with  camera  lenses  used  in  their  proper  mode  (infinite  conjugation  ratio). 
Actual  point  spread  functions  have  been  indirectly  measured  for  this  system  by  observing  the  image  of  a  point 
light  source  ( —5  pm  spot)  passing  through  the  slit  mask  at  successively  smaller  slit  widths.  With  the  camera 
lenses  the  light  continues  to  pass  the  slit  mask  with  the  full  FI. 5  down  to  25  pm. 
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The  final  optical  element  is  an  auto-collimator  also  referred  to  as  the  velocity  compensation  lens.  As  this  r 
implies,  the  function  of  this  lens  is  to  eliminate  the  effect  of  translational  motion  of  the  real  image  of  the 
particle  in  the  image  plane.  A  lens  of  focal  length  f  is  placed  at  the  distance  f  from  image  plane  to  lens 
pnncipal  plane,  and  also  distance  f  from  principal  plane  to  the  detector  surface.  The  position  of  the  light  s 
on  the  detector  is  then  determined  solely  by  the  angle  at  which  the  reflection  exits  the  image  plane.  The  ai 
collimation  principle  is  demonstrated  by  the  illustrative  construction  in  Figure  10.  The  determination  of  an 
velocity  is  not  affected  by  particle  translation  within  the  sampled  volume.  The  auto-collimator  for  this  sysi 
a  40  mm  focal  length.  AR -coated  F2.S  precision  achromat. 


Slit  Mask  in  Auto-collimator  Detector 

image  Plane  Lens  Plane 

greatly  exaggerated) 

vertical  Scale  Flgure  i0  b-sms 

Auto-Collimator  Function  Illustrated.  The  position  of  the  rays  on  the  detector 
plane  depend  on  exit  angle  8  only,  not  on  position  in  the  image  plane 

3.1.3  Data  Acquisition 

The  four  outputs  of  the  PSD.  the  x-sum,  x-difference,  z-sum.  and  z-difference  are  acquired  with  a  LeCroy 
digitizer  module  in  a  Camac  crate.  The  four-channel  6810  has  a  512  K  word  memory  and  a  variable  samf 
rate  of  up  to  1 M  samples/s  simultaneously  on  four  channels.  Each  channel  writes  a  minimum  IK  block  oi 
1 2-bit  data  at  a  uniform  sampling  interval  for  each  trigger  received.  The  trigger  signal  indicates  the  prese: 
a  reflection  on  the  detector  surface  and  the  level  is  set  to  exclude  signals  below  a  desired  signal-to-noise  n 
Wrap-around  buffers  allow  pre-triggering  so  that  data  is  symmetrically  acquired  about  the  trigger  point, 
memory  is  filled  after  writing  only  128  records.  Data  acquisition  must  then  cease  until  this  data  is  transfe 
to  computer  disk  via  the  Camac  GPTB  interface.  This  is  sketched  in  Figure  11. 
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4  •  .0  <SUM>  n  I  (INCIDENT  LIGHT  INTENSITY) 

A  .0  (OIFFI  n  I  •  XCIX  c-  X  CENTROID  OF  INCIDENT  LIGHT) 


'A' 'B  SUM 

e 

Figure  11 


Data  Acquisition  Path  of  the  Three-Component  VOP  System. 
This  system  has  been  used  for  development  and  testing  purposes  only 


This  data  acquisition  scheme  is  far  trom  optimum.  Custom-designed  hardware  in  future  VOP  designs  could 
lead  to  real-time  acquisition  and  processing  capability.  However,  for  the  purposes  of  instrument  development 
the  6810  performed  adequately  and  was  easily  incorporated  into  the  development  prototype  VOP. 


3.2  Data  Reduction 

Data  reduction  involves  many  separate  operations  which  are  listed  in  order  below: 

1 )  Background  subtraction 

2)  Thresholding 

3)  Masking 

4)  Detector  mapping  corrections 

5)  Optical  mapping  corrections 

6)  Signal-to- noise  ratio  statistical  weighing 

7)  Optimum  differencing  interval 

8)  Linear  least-squares  fitting  procedure. 

These  steps  are  briefly  described  in  the  following  paragraphs.  After  these  steps  have  been  performed,  several 
diagnostics  are  used  to  evaluate  the  validity  of  the  measurement;  these  are  discussed  in  Subsection  3.3. 
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All  of  the  dau  reduction  and  analysis  has  been  performed  on  either  a  286  or  386  PC  throughout  the  Phase 
program.  Even  some  of  those  features  which  might  be  considered  'signal  processing"  and  suitable  for  bar 
implementation  have  been  cast  initially  in  software  for  convenience.  For  example,  background  subtraction 
quotients,  and  thresholding  have  all  been  previously  done  with  analog  signal  processors  (Ferguson  and  We 
1983).  These  were  not  used  in  the  Phase  II  program,  however,  in  order  to  accommodate  the  available  dat 
acquisition  system. 

3.2.1  Background  Subtraction 

Unfortunately,  valid  reflections  are  not  the  only  source  of  light  falling  on  the  detector  surface.  The  PSD 
detector  is  not  an  imaging  system  but  rather  an  averaging  system,  measuring  the  centroid  of  the  light  distr 
over  its  surface.  Several  background  sources  and  their  relative  importance  can  be  identified  in  two  catego 
DC  or  slowly  varying  sources,  and  rapidly  fluctuating  sources  (on  the  time  scale  of  the  trajectory  passage; 
the  First  category  is  direct  or  stray  specular  reflections  from  the  incident  beam  which  can  potentially  overv 
the  valid  signal.  Also,  volume  scattering  (Rayleigh,  fine  impurities,  and  fluorescence  of  working  fluids)  f 
the  sampled  volume,  which  is  ordinarily  orders  of  magnitude  smaller  than  a  valid  reflection.  The  latter 
contributions  may  grow  relative  to  the  reflection  intensity  for  larger  sampled  volumes.  In  the  second  catei 
are  multiple  valid  reflections,  i.e.,  more  than  one  reflection  on  the  detector  at  the  same  tune  giving  overla 
trajectories,  and  reflections  or  scattering  from  VOP  particle  surfaces  or  large  particulate  contaminants  or 
bubbles.  The  former  are  of  the  same  order  of  brightness  and  violently  perturb  trajectories,  while  the  latte 
usually  weak  and  affected  trajectories  often  cannot  be  distinguished  from  'clean'  trajectories. 


The  first  category  of  background  disturbance  is  easily  remedied.  Because  it  is  constant  or  slowly  varying, 
(in  time)  background  signals  can  be  measured  and  subtracted  from  the  sum  and  difference  signals,  e.g., 


XDIFF 


XSIS  *  XBIB  ,  XSUM  =  Is 


XS1S  ^  XBIB 


where  Xj  is  the  centroid  of  signal  plus  background,  Xg  is  the  desired  x -position,  and  I§  the  valid  reflectu 
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optical  power,  and  Xg  and  lg  are  associated  with  the  background.  The  background  values  can  be  estimated 
from  the  sum  and  difference  baseline  signals  prior  to  the  event  so  the  correct  values  Xj,  and  1$  can  be  obtained. 
The  actual  background  levels  are  measured  by  seeking  the  most  probable  digitized  value  in  the  individual  1024 
point  records.  This  is  invariably  the  best  baseline  estimate. 

The  second  category  of  disturbance  is  not  fixable.  Multiple  reflection  (coverlapping  trajectory)  events  must  be 
discarded.  Weakly  perturbed  events  must  pass  goodness -of-fit  tests  (see  Subsection  3.3.1  below)  to  be  accepted 
as  good  data. 

3.2.2  Thresholding 

Of  the  1024  point  record  stored  by  the  6810  module,  at  most  only  a  few  hundred  of  the  interior  points  represent 
actual  trajectory  data.  Occasionally,  multiple  events  occur  within  a  single,  recorded  trace.  The  algorithm  must 
decide  where  those  events  are  located  in  the  data  stream  and  use  the  remainder  for  background  corrections. 

This  is  often  a  more  complex  task  than  expected,  particularly  when  the  data  is  characterized  by  extreme 
intensity  fluctuations.  Also,  some  data  sets  are  characterized  by  rapidly  rising  and  falling  edges  whereas  the 
opposite  is  frequently  true  for  equally  valid  events.  The  goal  of  the  thresholding  algorithm  is  to  keep  the 
maximum  amount  of  valid  trajectory  without  introducing  artifacts  or  unnecessarily  clipping  or  chopping 
trajectories. 

When  the  light  power  exceeds  a  preset  threshold,  the  x  and  z-positions  are  calculated  (from  the  quotients  of  the 
digitized  difference  and  sum  signals).  The  thresholding  algorithm  currently  in  use  compares  an  absolute 
threshold  value  with  a  running  n-point  smooth  on  the  sum  of  the  digitized  X-sum  and  Z-sum  signals  (the  total 
incident  power)  with  variable  hysteresis.  The  value  of  n  is  chosen  sufficiently  large  for  each  data  set  to  give 
adequate  smoothing.  This  algorithm  eliminates  rapid  on-off  threshold  flags  due  to  noise  near  the  selected 
threshold  voltage,  and  suppresses  the  possibility  of  splitting  a  single  trajectory  into  multiple  parts  due  to  large 
light  level  excursions.  It  has  the  virtue  of  approximating  a  band-limited  Schmidt  trigger  with  hysteresis  and  thus 
corresponds  to  a  hard ware-implemen table  feature. 
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3.2.3  Maskm 


From  the  discussion  in  Subsection  2.4.4,  one  consequence  of  the  finite  diffraction  angle  from  the  lead  carh 
platelets  is  that  retlections  have  a  finite  dimension  on  the  detector  face.  Frequently,  the  threshold  voltage  1 
exceeded  before  the  (Airy  disc)  reflection  is  fully  on  the  active  area.  This  clearly  gives  a  false  centroid'va 
An  effective  solution  is  simply  to  restrict  the  position  values  that  a  valid  trajectory  can  take  on.  by  ^  < 
software  mask.  A  band  around  the  perimeter  of  the  detector  active  area,  which  is  within  an  average  retlec 
diameter  of  the  detector  edge,  is  likely  to  suffer  the  edge  effect.  All  measured  positions  within  this  band  r. 
are  disallowed  as  possible  trajectory  coordinates.  This  is  implemented  m  software,  just  as  is  the  threshold 
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i  UDT  PIN  SC-25)  used  tn  (his  program,  showing  tines  of  real  physical  position  at  1  mm  separation  on  a 
detector  output  grid  square.  The  calibration  was  done  manually  by  traversing  a  HeNe  laser  beam  across  the 
detector  with  the  precision  translation  stage.  Using  this  map,  the  data  reduction  program  currently  performs  an 
interpolation  on  the  non-orthogonal  grid  to  find  true  position  coordinates. 

3.2.5  Optical  Calibration 

Calibration  of  the  optical  components  (the  imaging  lenses  and  the  auto-collimator)  was  performed  using  a  HeNe 
laser  and  a  60X  microscope  objective  as  a  source  at  the  nominal  sampled  volume  position.  This  reasonably 
approximated  a  point-source  and  filled  the  aperture  of  the  55  mm  collector  lens.  A  fine  reticle  with  a  square 
end  was  placed  perpendicular  to  the  optic  axis  so  that  its  shadow  also  filled  the  aperture.  The  optical  system 
was  properly  focused  and  35  mm  film  was  positioned  at  the  detector  plane  and  exposed.  One  such  photo  is 
included  as  Figure  13.  The  axi-symmetncally  distorted  grid  image  was  manually  digitized  and  fit  to  a  single 


V-501 


Figure  13 

Optical  Calibration.  The  shadow  of  a  square  grid  passing  through  the  optical  element 
provides  a  map  of  optical  distortion 
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parameter  equation  in  r.  the  distance  from  the  optic  axis  m  the  detector  plane.  Based  on  the  assumption  t 
Mmpie  field  curvature  caused  the  distortion,  the  best  fit  was  given  by. 


•actual  =  1-284  sin 


Selector 

1.284 


This  correction  effectively  flattens  the  field.  Figure  14  is  a  typical  digitized  trajectory  for  which  all  of  th 
preceding  corrections  have  been  made  (Subsections  3.2.1  through  3.2.S). 


Figure  14 

Digitized  Representation  of  a  Corrected  Trajectory 


3.2.6  Signal-to-Noise  Statistical  Weights 

As  mentioned  m  Subsection  3.2.3,  light  level  excursions  can  be  quite  extreme.  It  is  not  unusual  to  find 
variation  of  more  than  an  order  of  magnitude  over  a  data  set.  Since  the  noise  power  is  fixed  at  given  gai 
variations  in  signal -to-noise  ratio  (SNR)  are  correspondingly  large.  (The  calculation  of  the  SNR  random 
variable,  and  estimates  of  its  distribution  versus  system  parameters  are  given  in  Appendix  B.)  Clearly,  n 
data  points  should  be  equally  weighted  in  the  fitting  routine.  Their  proper  statistical  weights  are 


and 


Vi 


corresponding  to  the  dependent  variables  0',  and  O',  .  Propagating  the  errors 


detector  signals  through  the  computation  of  6'  and  0'  and  gives 
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D  fl  ♦  x2 
2  SNR-At-yd 


(8) 


D  \J  1  »  z" 
7  SNR-Afyd 


V 


where  At  is  the  sampling  interval.  yd  is  the  y  distance  from  the  imaged  sampled  volume  to  the  detector  plane, 
and  D  is  the  diameter  ot  the  detector  active  area.  Henceforward,  all  distances  on  the  detector  surface  are 
normalized  by  yd,  with  the  system  F-number  yd/D  -  1.  Lengths  are  thus  equivalent  to  angles  on  the  unit 
sphere. 


By  accounting  for  the  true  instrumental  noise,  the  X„2  goodness-of-fit  parameter  described  in  Subsection  3.3. 1 
below,  takes  on  useful  physical  meaning  and  is  used  to  evaluate  how  well  the  trajectory  equations  and  a  single 
vorticitv  vector  represent  a  traiectorv  at  any  instrumental  noise  level.  Uncertainties  in  the  vorticity  components 
properly  reflect  both  the  instrumental  noise  and  the  fit  quality. 


3.2.7  Optimum  Differencing  Interval 

Consider  an  arbitrary  trajectory  traversing  the  detector  with  a  given  constant  vorticity.  In  general,  it  will 


display  some  local  rate  of  passage  on  the  detector,  r  ,  and  some  angle  of  inclination,  \{/,  with  respect  to  the 

detector  plane  coordinate  axes  (x.z),  as  well  as  some  finite  curvature  —  .  Indeed,  recalling  the  discussion  of 

dr 

Subsection  2.3.2,  the  measurement  of  r,  i/-,  and  ££  represent  a  natural  geometrical  basis  for  the  description 

dr 

of  u'x,  w’z,  and  u )'y.  For  a  trajectory  traversing  the  entire  detector  surface  ( -  1  radian),  sampled  at  uniform 

1  1  To 


intervals.  At,  the  number  of  sampled  points  in  the  digitized  record  can  be  estimated  by  N  - 


r  At 


measure  the  mean  value  of  r,  denoted  r  .  over  the  entire  trajectory  with  the  maximum  achievable  precision,  the 


optimum  differencing  interval  Atopt,  on  which  the  time  derivatives  are  based,  must  be  selected. 


The  Taylor  expansion  in  time  of  the  trajectory  coordinate  r  about  rQ  is 


r“  Q 
At 


\ 

2 


?o  At 
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Lut  co  first  order,  the  uncertainties  in  r  and  r  .  denoted  6r  and  6r  .  estimated  from  Eq.  (8),  are  given  tv 


6r  =  a  &r~  -IL 

SNR-At  '  ysr 

where  a  is  a  constant  of  order  umty.  Thus,  with  the  estimate  of  N  above 

&  f  a  1 

/aT 

Clearly  as  At-*0.  the  loss  in  precision  in  r  is  more  rapid  with  decreasing  At  than  the  corresponding  mere 

the  number  of  measurements.  N,  can  offset,  and  5r  grows.  On  the  other  hand,  if  the  intervals  become  i 
great,  errors  due  to  second  order  (acceleration)  terms  in  the  expansion  begin  to  appear  as  estimates  of  the 

derivatives  become  increasingly  poor.  The  optimum  differencing  interval  is  defined  as  that  for  which  the 
sources  ot  error  are  comparable  in  magnitude: 


6r  - 


SNR-At, 


*  If  I  At, 


opt 


opt 


d/  —2  d\l/ 

where  ? ,  a  measure  of  local  acceleration  on  the  trajectory,  is  order  dl-x  .  —I  for  a  typical  trajectory  \ 

d  dr 


order  umty,  so  Eq.  (9)  gives  an  optimum  differencing  interval,  Atopt, 


At  -  1 

alopt  -  , - 

rv/SNR 

or  the  number  of  intervals  for  a  full  trajectory,  nopt 

nopt  ~  ^NR 


and  therefore,  from  Eqs.  (9)  and  (10), 


v/SNRv/N- 


This  analysis  can  be  extended  to  and 


,  and  include  partial  trajectories  as  described  in  Appendix  C 
3r 


4? 


In  terms  of  the  data  reduction  algorithm,  precision  is  considerably  improved  when  derivatives  are  computed 
.'ver  m-point  intervals  where  At0p,  =  mAt,  e.g.. 


j ,  *m*l  “  el  a,  em+2  ~  e2 

tt  i  =  -  .  t)  -t  =  . . 

mAt  “  mAt 

This  approach  will  not  yield  as  great  a  precision  in  w.  or  be  as  statistically  sensitive  to  goodness-of-fit  tests,  as 
non-lmear  regression  directly  on  trajectories,  but  it  is  computationally  many  times  faster. 


3.2.8  Linear  Least-Squares  Fitting  Procedure 

The  least-squares  fit  of  the  data  to  the  differential  trajectory  equations  (DTEs)  is  conventional.  The  three 
vorticity  components  are  obtained  through  the  accumulation  of  appropriate  sums  and  3x3  matrix  inversion,  a 
simple  extension  of  standard  linear  least-squares  methods.  The  linear  DTEs  can  be  written 


d’ 


d>'  =  B-tc 


and  x~  is  therefore 


X 


-  O'fk^k 

*9, 


2 


% 


The  linear  least-squares  technique  determines  the  w  which  minimizes  i.e.. 


duij 


N 


E 


Qt  k"k 


2 

<t>( 


or,  dropping  the  l  index  (=  1....N)  with  the  understanding  that  weighted  summation  is  implied 


(aj'  -  or;*')  =  [a^  *  ,  (i.k*  1,2,3) 
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The  solution  is  provided  by  inverting  the  curvature  matrix.  =  lQ,ak  -  ¥kl 

“Tc  =  [ai«k  *  MkfV'  *  l4iO') 

Uncertainties  are  estimated  from  the  elements  of  the  curvature  matrix. 

3.3  Error  Recognition  and  Rejection 

Some  of  the  trajectories  collected  by  the  VOr  system  will  be  contaminated  or  rendered  useless  by  one  01 
of  the  following  effects: 

1 )  Traiectorv  overlap. 

Two  or  more  reflections  on  a  detector  render  these  events  umnterpretable.  The  probability  of  ti 
overlap  is  data  rate  dependent  (see  Subsection  3.4.3,  below). 

2)  Detector  edge  effects. 

Large  reflection  divergence  angles  often  attributable  to  multiple-mirror  clusters  within  particles, 
spurious  reflection  coordinates  and  speeds  while  the  reflection  is  moving  onto  the  detector  activi 
Detector  masking  accounts  for  most,  but  not  all  of  these  events. 

3)  Non-sphencal  particles. 

Malformed  or  clustered  particles  cause  modulation  of  the  particle  rotation  rate  relative  to  the  Ich 
vorticitv.  (It  is  possible  to  sample  the  rotation  of  such  a  particle  when  df)/9t  is  large.  This  situ 
difficult  to  distinguish  from  a  true,  local  angular  acceleration,  du/d t.  due  to  turbulent  flow 
fluctuations.) 

4)  Traiectorv  fragments. 

Under  some  conditions,  fragments  which  are  too  small  to  provide  reliable  information  on  trajeci 
curvature  and  accelerations,  couple  the  uncertainty  in  u'y  strongly  into  u'x  and  u'  r  This  is  rel 

5)  The  zero  vorticitv  problem. 

A  stationary  reflection  on  the  detector  can  arise  in  the  case  of  zero  vorticity,  a  finite  vorticity  vi 
parallel  to  the  mirror-normal  vector,  or  diffuse  reflection  from  a  contaminant  particle. 
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Trajectories  suffering  from  any  of  the  above  complications  are  referred  to  as  irregular  trajectories.  Some 
strategies  designed  to  identify  such  errors  and  reject  correct  or  flawed  trajectories  are  discussed  below. 


3.3.1  The  Vj-  Statistic 

The  basis  of  least-squares  fitting  algorithms  is  the  minimization  of  the  x2  statistic,  whether  by  linear  or  non¬ 
linear  regression.  The  x2  variable  is  defined  as 

X2  =  Z,2  +  Z22  +  ...  +  Zn2 

where  Z,  are  independent  random  variables  normally  distributed  with  zero  mean  and  unit  variance  ( Bendat  & 
Piersol.  1971).  In  the  content  of  curve  fitting,  such  a  variable  is  usually  constructed  from  the  measured 
dependent  variable.  y  ,  the  known  variance  of  its  Gaussian  noise,  Oj,  and  the  assumed  functional  form  for 
ylx....{p}),  evaluated  at  the  independent  control  variablefs),  x .  We  obtain 

,  (y;  -  y(*i,...))2 

L,  =  5 

or 


2  £  (y\  -  y<*i....>) 

*  '  A  2 

1=1  <7 


The  best  fit  is  defined  as  the  set  of  parameters  {p}  in  y(Xj  ,{p})  which  result  in  the  minimum  attainable  x2- 
The  number  of  degrees  of  freedom  for  a  fit  with  p  independent  parameters  is  v  =  n-p  and  the  expected  value  of 


2L  -  l  . 


when  the  variances  <Tj2  have  been  correctly  estimated,  and  when  the  deviation  from  y(xj  _,{p})  is  attributable 
entirely  to  noise.  If  the  former  condition  is  met,  increases  in  x,,2  are  assumed  to  arise  from  systematic 
deviation  of  the  measured  values  from  the  fitted  expression.  This  is  the  basis  of  goodness-of-fit  tests.  A  large 
X,2  value,  t.e.,  x2  >  >  1,  is  a  strong  indication  of  a  poor  fit  to  the  data  because  the  function  y(Xj  {p}) 
cannot  describe  the  data  for  any  parameter  values. 
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Unfortunately.  \J~  is  not  an  absolute  test  of  the  adequacy  of  the  fitting  function.  The  distribution  of  sig 
noise  ratios  of  trajectories  in  a  data  set  can  be  broad.  Identical  trajectories,  each  containing  the  same  oh 
detect,  are  represented  in  Figure  15.  Figure  15(a)  possesses  a  large  SNR  whereas  Figure  15(b)  has  a  io' 
and  instrumental  position  noise  is  comparable  in  variance  to  the  deviations  from  the  true  trajectory.  (Sue. 
error  might  be  due  to  a  second  reflection  striking  the  detector  while  the  first  is  traversing  the  detector.  ! 
eases,  the  centroid  of  the  pattern  is  strongly  deflected  from  a  legitimate  or  regular  trajectory  shape.)  In 
case,  the  best  fit  curve  shown  as  the  solid  line,  is  likely  to  be  very  similar.  On  the  other  hand,  if  is 
partitioned  into  the  contributions  from  systematic  deviations,  3 ,  from  the  fitted  curve  and  those  due  to  r 

noise  o.  a  1  and  on  ot  — ! —  in  Figures  15(a)  and  15(b),  respectively,  then 

A  snra  snrb 


2 

Xj, 


i 

<r 
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But,  from  the  discussion  above,  a  ^  <  3  and  ag 


3.  therefore 


=  1 


2 

X*B  = 


*  2 


\J-  can  be  quite  insensitive  as  a  goodness -of-fit  parameter  when  the  precision  of  the  measurements  is  poor  and 
parameter  uncertainties  are  large. 


This  situation  is  exacerbated  when  fitting  the  differential  trajectory  equations.  Noise  magnification  of  the  aA 
and  oB  values  as  a  result  of  differentiation  often  eliminates  the  goodness-of-fit  utility  of  the  x~v  test  entirely.  In 
this  case,  an  auxiliary  approach  is  needed  for  trajectory  evaluation. 

3.3.2  Orthogonal  Polynomial  Decomposition 

The  differential  trajectory  equations  (DTEs)  describing  a  single  reflection  and  constant  vorticity  can  be  fit  to  any 
measured  trajectory.  Not  all  measured  trajectories,  however,  result  from  single  reflections  at  constant  vorticity. 
Irregular  trajectories,  those  which  could  not  have  been  produced  by  the  equations  or  which  are  pathological 
(near  the  VOP  singularity),  must  be  identified  and  rejected.  A  fast  algorithm  is  needed  to  quantify  trajectory 
shapes  and  decide  whether  they  belong  to  the  subset  of  legitimate  or  regular  trajectories  which  the  DTEs  can 
produce.  Trajectories  are  described  by  (fl'(tj),  <p'(tj)).  the  set  of  N  angular  coordinates  of  the  reflection  sampled 
at  fixed  time  intervals.  At.  The  O'  and  <t>'  projections  of  each  trajectory  are  easily  fit  to  generalized  orthogonal 
polynomials  using  a  very  efficient  linear  algorithm  (Bevington,  1969). 

Orthogonal  polynomials  are  defined  as  a  basis  set  of  functions  PQ(t)  on  an  interval  [tj,  t2]  such  that  an  arbitrary 
function  0(t)  can  be  expressed  as 

anPn(t)  with  orn  -  f  Pn(t)0(t)dt  and 
*1 


m » E 

n 


j  P„(t)Pm(t)dt»0  for  n*m 
1 


47 


In  other  words,  the  an  are  independently  specified  and  do  not  alfect  the  other  am,  m;*n.  The  polynoi 

S'(tj)  =  a0  +  ajt^-tj)  -r  Orit^-t-,)2  f  ... 

<?•(!,)  =  i)0  +  hjitj-t,)  +  rt2(trt2)2  +  ... 

where  the  tn  depend  only  on  the  number  of  points,  N,  oq  represents  the  mean  8'  value,  aj  gives  the  av 
slope,  a-,  the  average  curvature  and  so  on.  Based  upon  the  symmetries  of  the  VOP  geometry  and  trajei 
equations,  polynomials  up  to  the  third  order  adequately  span  nearly  all  of  the  allowed  trajectory  space, 
an  arbitrary  trajectory  shape  can  be  expressed  as  a  vector  in  8-dimensional  vector  space 
aQ-aj, a2.a3.Ao.fii, 62,83]  =  |a,Bj.  \v  -  are  also  generated.  Some  representative  orthogc 

polynomial  decompositions  (OPDs)  are  reproduced  in  Figure  16. 


The  figure  gives  the  impression  that  the  OPDs  span  a  la, uj  space  that  is  much  larger  than  that  spanned 
DTEs:  Figure  16b  is  not  a  legitimate  (regular)  trajectory.  To  reject  this  trajectory,  a  (7-dimensional)  f 
could  be  found  tn  !a.u]  within  which  trajectories  are  considered  to  be  regular.  However,  with  the  assu 


B-7 


(a)  Regular  Trajectory  (b)  Irregular  Trajectory 

(i.e.,  not  producible  with  the  DTE 


Figure  16 

Some  Representative  Trajectories  with  their  Orthogonal  Polynomial  Representations 
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that  all  regular  trajectories  are  contained  in  [a.fl],  a  trajectory  that  cannot  be  fit  by  the  third  order  OPD  is 

unlikely  to  be  regular.  It  follows  that  the  \ ?  goodness-of-fit  statistics  for  the  OPD  can  be  immediately 
employed  to  reject  many  irregular  trajectories.  These  statistics,  computed  from  trajectories  directly  rather  than 
their  time  derivatives,  are  more  sensitive  tests  of  goodness-of-fit,  as  noted  in  the  preceding  section.  The  regular 

trajectory  boundaries  in  [a.hj  have  not  yet  been  fully  explored. 

A  useful  extension  of  the  OPD  xy~  test  for  edge  effects  is  used  routinely.  For  OPD  values  above  a 
threshold,  the  trajectory  is  re-analyzed  with  a  smaller  detector  area.  An  additional  swath  of  10  percent  of  the 
detector  diameter  around  its  perimeter  is  added  to  the  software  mask.  If  the  OPD  x„~  value  drops  below  yet 
another  threshold,  the  error  is  assumed  to  have  been  edge  effect-related  and  this  trajectory  is  analyzed  over  the 
new  smaller  detector  area.  Otherwise  it  is  rejected. 

3.3.3  Proximity  Tests 

The  last  two  classes  of  problems  noted  above,  trajectory  fragments  and  zero  vorticity,  are  closely  related.  The 
proximity  tests  are  so-called  since  they  monitor  the  proximity  of  the  measured  vorticity  to  the  VOP  singularity: 
w  parallel  to  the  mirror-normal  vector.  Both  zero  vorticity  and  the  shortest  trajectory  fragments  may  suffer  the 
effects  of  the  VOP  singularity,  the  former  by  virtue  of  actually  being  parallel  to  the  mirror-normal,  and  the 
latter  by  spunous  curvature  estimates  at  finite  SNR.  The  discussion  of  the  actual  implementation  of  these  tests 
has  been  deferred  to  Appendix  D.  The  parameters  upon  which  these  tests  are  based  must  first  be  introduced. 

3  4  System  Characteristics 

The  three-component  VOP  measurement,  though  simple  in  concept,  involves  a  number  of  subtleties  in  analysis 
and  interpretation  which  complicate  the  experimental  technique  in  practice.  In  this  respect  the  VOP  technique  is 
rather  like  Laser-Doppler  Velocimetry  (LDV)  which  suffers  the  effects  of  finite  size  sampled  volumes,  noise, 
sampling  bias  due  to  velocity  fluctuation,  etc.  The  VOP  system  measures  a  reflection  from  mirrors,  log- 
normally  distributed  in  diameter,  moving  in  a  fluctuating  turbulent  velocity  field,  passing  through  a  gaussian 
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beam  protile  at  random  position.  The  reflection  (and  sometimes  tluctuating  background  light)  strikes  the 
detector  at  a  random  initial  <X,Z)  position,  for  a  time  period  (and  trajectory  size  )  dependent  upon  the  v 
vector,  the  velocity  vector,  the  position  of  the  particle  in  the  sampled  volume,  and  the  mirror  diameter, 
ot  the  independent  random  variables  (less  the  number  of  relations  correlating  some  of  them)  in  the  afeov 
description  is  a  convincing  argument  for  care  in  analyzing  and  interpreting  data.  On  the  other  hand,  tht 
also  quite  rich  in  information  about  the  vorticity  and  velocity  fields,  and  yields  simple  parameters  which 
used  to  monitor  instrumental  precision. 

The  following  subsections  outline  the  necessary  correction  to  raw  vorticity  component  PDFs,  as  well  as 
guidelines  for  interpreting  the  statistics  obtained  from  these  distributions.  First,  an  analysis  of  the  expec 
distribution  at  constant  u  is  presented.  This  is  the  resolution  function.  Next,  a  correction  for  vorticity  s. 
non-uni formitv  due  to  both  instrumental  and  flowfield  parameters,  called  the  sampling  function,  is  descr 
This  function  accounts  for  vorticity  variation  at  fixed  instrumental  resolution.  A  method  for  correcting 
vorticity  PDFs  using  both  the  resolution  and  sampling  functions  is  outlined.  Lastly,  issues  relating  to  sj 
and  temporal  resolution  are  considered. 

3.4.1  Vorticity  Resolution  Function 

The  vorticity  resolution  function  is  defined  as  the  distribution  that  would  be  obtained  by  measuring  a  vo 
{-function.  The  width  of  this  distribution  is  due  entirely  to  instrumental  noise,  and  it  is  this  function  th. 
convolved  with  measurement  of  fluctuating  vorticity  which  must  be  accounted  for: 

R|w0)  =  £  R(w)  -  u0)d3<J 

< 4 ) 

and  the  measured  PDFs: 

Pmeasured<">  =  J_  P<«')  R(“  ~  «')  d3(i/ 

0) 
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to  general,  the  form  of  R(u)  will  vary  with  u,  but  it  is  not  possible  to  measure  a  vorticity  6-function  at  every 
point  in  vorticity  space  in  order  to  catalogue  the  behavior  of  R(w).  On  the  other  hand,  if  the  uncertainties  in 
vorticity  measurement  and  their  distributions  can  be  described  as  a  function  of  u  and  relevant  experimental 
parameters,  then  R(  ui)  can  be  directly  calculated.  In  the  following,  separability  of 
R(u>)  *  Rx(wx)  -  Ry(uiy)  -  R^fuij)  has  been  assumed. 


Approximate  geometrical  relationships  elicit  the  essential  features  of  the  VOP  trajectory  equations.  The 
uncertainty  of  any  particular  vorticity  measurement  can  be  estimated  from  relations  of  the  vorticity  components 

to  r,  and  (the  alternate  vorticity  basis  described  in  preceding  sections).  These  values,  in  turn,  are 
dr 

simply  related  to  the  basic  trajectory  data,  where  errors  are  derived  from  known  noise  and  trajectory 
characteristics.  When  the  noise  and  trajectory  characteristics  (i.e.,  angular  size,  etc.)  are  distributed  in  known 
ways,  depending  upon  experimental  parameters,  then  all  of  the  information  needed  to  generate  a  resolution 

function  is  at  hand.  The  propagation  of  measurement  error  through  r,  \ p,  and  ^  is  described  in  Appendix  C. 

dr 

Since  the  uncertainty  of  given  vorticity  measurement  is  based  upon  these  parameters,  it  is  important  to  known 
with  what  probability  a  given  SNR  or  rjsj  (and  thus  vorticity  component  uncertainty,  a;)  will  occur.  The 
distributions  of  signal-to-noise  ratio  and/or  transit  times,  and  trajectory  angular  size  are  derived  in  Appendices  B 
and  E,  respectively. 


Figure  17  shows  the  geometry  of  an  arbitrary  trajectory.  The  figure  is  a  representation  of  the  unit  solid  angle 
chord  (detector)  surface  on  the  unit  sphere,  where  a  is  the  distance  of  the  centroid  of  the  trajectory  from 

detector  center  (x=0,  z=0),  y  is  the  distance  from  the  centroid  to  the  effective  center  of  curvature  of  the 
T  33"]  -1 

trajectory  with  7  —  ,  and  c  is  the  distance  from  detector  center  to  the  apparent  center  of  curvature, 

[drj 
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Figure  17 

Trajectory  Geometry 


For  modest  |  c  |  values  (  $  1),  the  approximate  rotation  rate  of  the  reflection  about  the  apparent  center 
curvature  is  denoted  wc  =  f— ,  and  the  vorticity  components  can  be  estimated  by 


In  terms  of  the  alternate  basis  r,  i£,  and  Eq.  (12)  are 
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•a here  (he  i . ,  notation  in  Eq.  (14)  indicates  the  square  root  of  the  sum  of  the  squares.  The  approximate  total 
length,  r^-,  ot  the  N -point  traiectorv.  is  the  fraction  of  the  (1  radian)  angular  size  of  the  detector,  and  r.  is  the 
minimum  angular  spacing  between  successive  trajectory  points.  As  expected,  the  curvature-related  contribution 
to  the  uncertainty  in  m'x  and  ui'z  vanishes  at  detector  center  (az  =  0,  ax  =  0).  In  this  case,  the  second 
(smaller)  term  dominates  in  the  expression  for  6oj'x  and  6oi'z. 


In  the  large  |  c  |  limit  in  which  vorticity  vectors  are  nearly  perpendicular  to  the  mirror  normal,  the  u>c 
approximation  tails.  In  this  case.  :uj  *  r,  and  c  *  y  >  a.  giving 


The  large  |c]  limit  gives  the  same  uncertainties  in  the  vorticity  components.  This  error  analysis  is  therefore 
expected  to  be  valid  for  any  arbitrary  vorticity  vector. 
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Having  described  the  dependence  of  the  errors  on  trajectory  characteristics,  the  error  distributions  can  b< 
derived  in  £q.  ( 14)  for  the  component  uncertainties,  the  dominant  terms  m  the  av  a,  <x  (y'SNR  Tfjj  1  . 
(he  determination  of  the  distribution  p(Oj)  since  the  distributions  of  p<SNR)  and  pfr^)  are  known.  With  t 
power  threshold  set  sufficiently  high  (Subsection  3.2.2),  the  variation  in  SNR  produces  a  modest  unifdrrr 
hroadening  ot  the  a,  distribution.  The  variation  in  then  dominates  the  error  distribution.  This  distnb 
as  a  function  of  l  (the  total  angular  size  of  trajectory  on  the  unit  sphere)  is  derived  in  Appendix  E.  In 
Urge  i  approximation  the  distribution  of  trajectory  lengths  is  given  by 


0<rN<R.  f^R 


where  R  is  the  angular  measure  of  the  detector  radius.  (The  random  variation  of  ax,  and  »z  is  not  indue 
this  error  analysis,  and  is  yet  another  source  of  broadening.)  Since  a  oc  (i/sNR"t^)  * ,  we  obtain 


thus  for  large  a[t  Eq.  (IS)  is  approximately 


ffi,min  <  a\  <  °l,max 


The  broadening  due  to  the  SNR  distribution  (and  the  contribution  of  ax  and  az),  suffices  to  remove  the 
(integrable)  singularity  at  <Xj  min  so  that  Eq.  (16)  is  a  good  approximation  over  the  full  range.  The  uppe: 
limit  <r]  max  is  either  a  function  of  the  cutoffs  imposed  in  the  proximity  tests  (Appendix  D)  for  a  preset  \ 
This  error  distribution  estimate  is  borne  out  by  the  calibration  data  (single  vorticity  vector)  presented  in 
Chapter  4  where  the  error  distribution  for  each  cu  component  in  a  laminar  flow  are  presented.  The  know 
distribution  of  errors  leads  naturally  to  an  approximate  resolution  function  described  below. 
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i'nr  a  subset  or  measurements  or  the  quantity  x.  derived  from  a  normal  parent  distribution  characterized  by  a  , 
the  contribution  to  (he  total  distribution  tunction  in  the  range  Act.  is 


pxtx:o.) 


cJ. 


A  a 
— 


<  o,  <  a. 


Aa 


•shere  a,  is  the  amplitude  of  the  |th  contnbution  and  Wj  is  the  statistical  weight.  The  expected  form  of  the 
resolution  tunction  in  the  limit  or  a  continuous  a  distribution  is  then 
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These  resolution  function!  s)  are  denved  from  calibration  data  with  ffmin  and  <xmax  estimated  from  the 
experimental  error  distributions.  The  calculation  of  the  moments  of  this  distribution  is  straightforward. 


55 


The  moments  ot  the  normalized  Gaussian  distribution  are  x  =  0.  \~  =  <r .  -  0,x"^  =  3a4,  etc.  Then 
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The  vanance  and  tlatness  relative  to  a  normal  distnbution 
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Note  that  for  unbounded  amax,  the  resolution  function  width  ap,  is  v’T  time  greater  than  the  gaussian  val 
a_,n  =  Equations  (18)  and  (19)  are  plotted  versus  ma^  in  Figure  18. 
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Figure  18 


RMS  Distribution  Widths  and  Flatness  Factors  for  Resolution  Function  Versus 
Relative  to  the  Normal  Distribution  Values 


,Jmax 

arain 


It  is  usually  impossible  to  set  (the  mean  total  trajectory  angular  size)  to  an  arbitrary  high  value,  since  it  is 
determined  in  part  by  the  vorticity/velocity  ratio: 


=  |  U)|  Tjr  = 


(20) 


where  f  is  the  sampled  volume  dimension  parallel  to  the  flow  and  rtr  is  the  transit  time.  As  described  in 
\ppendix  B.  in  geometries  where  the  flow  is  inclined  to  the  beam,  the  variation  of  the  traversed  sampled 
volume  dimension  and  therefore  transit  time,  results  in  a  range  of  l  at  constant  o>.  Turbulent  velocity 
tluctuations  contnbute  to  the  f  ,  broadening  m  a  similar  fashion. 

For  modest  i  ^  in  a  turbulent  flow  a  better  estimate  for  p(ctj)  is 


W  *  T 
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Then  Eq.  (18)  becomes 
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These  resolution  functions  are  broader  and  their  moments  can  be  computed  as  before.  The  second  mome 
this  case  (Eq.  (211)  gives 
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The  resolution  function  width  diverges  weakly  with  unbounded  trmax.  though  is  still  less  than  3amm  a 


-  -  30.  It  is  convenient,  in  anv  case,  to  truncate  the  error  distributions  bevond  a  present  «r  „  levi 

ft  '  LU<lA 

umm 


order  to  produce  reliable  bounds  on  <rR  and  minimize 


In  general,  some  intermediate  form  of  p(ffj)  is  expected,  and  under  differing  experimental  conditions  the 
inclusion  of  other  factors  may  be  necessary,  e.g.,  SRN  and  t  broadening.  If  p(fw)  is  known,  then 


P(gi:r) 


P(f  J  drN 


where  Psnr  is  the  SNR  distribution  function  (Appendix  B), 


ar2 

tr2 


is  substituted  for  SNR  from  the  expi 


for  ar  and  pfr^:  f^)  is  computed  in  Appendix  E.  In  experiments  where  turbulent  flow  fluctuations  domi: 
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nstrumental  broadenim:.  tnc  detailed  torm  ot  pttr  )  is  obviously  of  less  importance  to  the  shape  ot  the  PDFs  and 
..dues  ot  the  moments. 

! t  must  be  emphasized  that  the  a  ,.)m  and  a.  m  x  limits  the  expression  tor  the  normalized  resolution  function  R. 
-cale  directlv  with  the  vorticitv  rhus  the  resolution  functions  vary  systematically  with  vorticity.  complicating 
the  deconvolution  of  instrumental  broadening  from  true  vorticity  fluctuations.  This  inhomogeneous  broadening, 
however,  preserves  the  integral  probability  of  vorticity  component  PDFs.  When  the  probability  of  sampling  a 
particular  vorticity  is  influenced  by  the  vorticity.  corrections  are  required  to  eliminate  this  bias.  In  order  to 
preserve  the  simple  form  ot  the  resolution  functions  above,  turbulent  vorticity  PDFs  are  corrected  for  vorticity 
dependence  via  sampling  tunctions. 

;  4  2  Vorticitv  Sampling  Function 

What  is  the  probability  per  unit  time  of  measuring  vorticitv  ui,  under  fixed  experimental  conditions,  and  how 
-hould  that  value  be  weighted  to  obtain  a  true,  unbiased  vorticity  distribution?  The  answer  to  this  question  con¬ 
stitutes  a  vorticitv  sampling  function.  The  true  vorticity  PDF,  p(  w)  is  related  to  the  measured  PDF.  pm(w,)  by 

Pmiwxl  3  j  S(w)  P(ux.wy,wzjdwydujz 

•  >r  example,  where  some  normalization  is  imposed  over  the  range  of  integration,  and 

p(w)  =  S  "lM  pm(w) 

where  S( uij  is  the  vorticitv  sampling  function.  Note  the  probability  of  measuring  uij  is  influenced  in  general  by 
the  values  of  all  components.  The  data  rate  at  any  point  in  vorticity  space  is  given  by 

R  =  nAu  S(w)  <22) 

where  n  is  the  probe  panicle  number  density  and  A  is  the  cross-sectional  area  of  the  sampled  volume  normal  to 
the  flow  velocity  u. 
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Proper  statistical  weights  are  given  to  vorticitv  components  when  generating  PDFs  and  statistics.  To  the  e 
that  the  statistical  weighing  factor,  !  ia~.  vanes  systematically  with  vorticitv.  this  must  be  corrected  for. 
sampling  corrections  to  the  PDFs  for  data  reduced  with  the  optimal  differencing  scheme  can  be  deduced  tr 
the  form  ot  a,. 

a,  <* 

The  explicit  appearance  of  r  in  tr,  can  be  corrected  directly  by  the  sampling  factor 

Sr(w)  «  (r)"2  *  \u'l  *•  (2u/x)2  \1\2  cc'""1  . 

A  more  subtle  contnbution  with  increasing  r  is  the  tendency  of  the  ryj  distnbution  to  favor  longer  trajectc 
As  seen  in  Appendix  E.  the  mean  trajectory  length  is  a  function  of  tv.  Vanations  m  r ,  and  therefore  I 
affect  the  PDFs  through  ryj.  (Also,  functional  the  form  of  the  resolution  function  is  weakly  dependent 
which  is  not  stnctly  a  sampling  function  type  correction  as  defined.)  In  the  large  approximation  discus 
ihe  preceding  section,  however,  this  effect  is  negligible  and  will  not  be  examined  here. 


i/SNR 


"rM 


An  important  sampling  function  factor  arises  from  the  increased  probability  per  unit  time  at  higher  vortici 
that  reflections  will  be  swept  into  the  detector  solid  angle.  Again,  for  sampled  volume  transit  time  rtr.  wc 
r  ~  i  which  leads  to 


Sq(w) 


r  x  ♦  4Fr  Ttr 

16tF2 


where  F  is  F-number  of  the  optical  system  (which  has  been  approximated  by  1  in  order  to  simplify  all  of 
foregoing  discussions).  The  first  terra  m  the  brackets  is  the  fraction  of  the  solid  angle  4t  subtended  by  tf 
detector  surface  and  the  second  term  is  the  additional  solid  angle  swept  out  by  the  detector  surface  in  time 
in  a  frame  rotating  at  r  *  *iJTtr  In  the  general  case  where  individual  rtr  values  are  unknown,  this  appr 
is  inappropriate.  Fortunately ,  an  equivalent  correction  exists,  based  on  detector  transit  tune,  rjy.  Detector 
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:  ran  .si!  ume  weighing  is  needed  in  the  ease  of  the  VOP  for  preeiselv  the  reasons  sampled  volume  transit  time 
weighing  is  required  in  Laser-Doppler  Velocimetry  (LDV).  A  proper  time  average  is  obtained  for  discrete 
r.mdom  events  ot  finite  duration,  like  VOP  data  and  'burst'-tvpe  LDV  data,  only  when  these  events  are 
weighted  bv  the  time  interval  over  which  they  are  observed  (Buchhave.  George,  and  Lutnlev.  1976) 

“true  =  Eu)iTDi  ’  ETD.  (25) 

i  i 


To  ventv  that  7q  is  the  proper  weighing  factor,  its  mean  value  can  be  estimated  using  the  same  arguments 
which  led  to  Eq.  (24) 


rD  = 


*  u 

luT 


4fF  t, 


tr 


(26) 


where  the  second  term  in  the  denominator  corrects  for  the  probability  that  the  reflection  leaves  the  detector  area 
in  time  r.r  At  zero  voracity,  this  expression  correctly  reduces  to  tq  =  t’u/  |  u  |  =  rtr.  The  mean  detector 
transit  time  weighted  data  rate,  combining  Eqs.  (22).  (24),  and  (26),  is  thus  dependent  of  u  and  ui. 


nA| u|  SqQ  rD  = 

16FZ 

where  V  is  the  sampled  volume.  Detector  transit  time  weighing  simultaneously  eliminates  both  vorticitv  and 
velocity  sampling  bias. 


This  correction  simultaneously  accounts  for  the  enhanced  probability  of  measurement  at  high  vorticity  and  the 
reduced  contribution  of  these  short  events  to  the  time  average.  This  is  illustrated  in  Figure  19.  In  the  sampled 
volume  transit  time  limited  case.  TQ  =  r,r,  and  detector  transit  time  is  identical  to  the  conventional  LDV  form. 
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Figure  19 

Transit  Time  Corrections.  Correct  averages  are  obtained  by  weighing  vorticity  values 
in  proportion  to  the  time  interval  during  which  the  particle  is  observed 

Potentially  the  most  difficult  sampling  corrections  are  those  near  zero  vorticity,  or  the  VOP  singularity  wf 

Sg(ai)  -»  0  and  pm(w)  -*  0  ,  as  |  oij  -»  0  . 

This  is  related  to  the  interpretation  of  near  zero  vorticity  by  the  proximity  tests  (Appendix  D)  and  has  not 
been  fully  clarified,  nor  has  any  explicit  form  of  SQ  been  derived.  In  such  cases  there  is  little  available 
information.  Assumptions  about  the  form  of  pftu'y)  are  probably  required  (i.e.,  that  the  vorticity  field  is  r 
dimensional  and  that  the  wy'  contribution  can  be  ignored)  in  order  to  fill  in  p(aj'x)  and  p(oj'z)  near  |  coj  = 
an  unbiased  fashion.  Statistically  significant  null  vorticity  caused  intermittent  irrotational  flow  (oi=0)  can 
be  distinguished  from  the  normal  background  errors.  This  is  only  expected  to  be  troublesome  for  turbulei 
flows  with  little  or  no  mean  shear  or  significant  intermittency,  and  will  be  pursued  in  future  work. 

Another  problematic  form  of  sampling  bias  is  the  direct  interaction  of  the  vorticity/velocity  ratio  with  the 
proximity  tests.  Trajectory  fragments  are  typically  obtained  when  the  convection  velocity  is  high  relative 
local  vorticity.  Since  proximity  tests  limit  the  allowable  contamination  of  u'z  and  w'x  measurement  by  t 

large  uncertainty  in  w'y,  two-component  fits  are  generally  acceptable  only  near  the  detector  center.  This 

effectively  smaller  detector  area  creates  a  sampling  bias  whose  effects  were  observed  experimentally  and  . 
noted  in  Chapter  5. 
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“he  complete  vonicitv  sampling  function  is 


Site)  =  Sq(ui) Sr(a;)  7^3 


(28) 


which  -an  be  incorporated  directly  into  the  statistical  weights  dunng  dau  analysis,  t  At  present.  S0  is  set  to 
unity.  i 


The  sampling  function  assures  that  no  vorttcuy  vector  is  systematically  over-  or  under-counted  m  the  measured 
PDFs.  However,  the  implicit  r  dependence  in  the  vorticity  resolution  function  results  in  inhomogeneous 
broadening  of  the  PDFs.  Formally,  ihe  deconvolution  of  the  PDFs  alluded  to  in  the  previous  section  should  be 
performed  on  the  three-dimensional  lomt  PDF,  p(aix,i*)v,u;z)  with  the  resolution  function  R(mv  wv,  tcz).  Since 
:ne  hreauth  ot  the  resolution  functions  scale  approximately  with  |  w  j ,  i.e..  Pui/w  *  constant,  the  deconvolution 
.an  be  linearized  by  transforming  to  logarithmic  coordinates  (scaled  by  Iniuii ).  The  final  PDFs,  p(tOj),  are  the 
proiections  ot  the  corrected  PDF  onto  the  coordinate  axes.  Such  a  computation  is  likely  to  be  impractical, 
however,  and  approximate  method  will  need  to  be  employed.  Empirical  corrections  to  the  moments,  for 
example,  as  employed  in  Chapter  5,  may  be  a  generally  satisfactory  solution. 


3.4.3  System  Spatial  and  Temporal  Resolution 

The  intrinsic  spatial  resolution  ot  the  VOP  (particle  diameter)  has  been  discussed  in  Chapter  2.  No  discussion 
has  yet  been  offered  of  the  effects  of  finite  sampled  volume  size.  Apart  from  the  direct  effects  on  the  vorticity 
statistics  due  to  spatial  gradients  of  the  mean  vorticity  field,  which  are  readily  correctable,  there  is  the  more 
important  issue  of  the  vorticity/ velocity  ratio,  -^fi.  If  the  velocity,  u,  becomes  large  at  given  vorticity  and 


i“ 


sampled  volume  dimension  (rJ  (parallel  to  u),  then  the  transit  time.  rtf,  drops  off  along  with  trajectory  angular 
size  r^t.  Since  a  much  smaller  angular  displacement  is  observed,  the  precision  of  vorticity  measurement  suffers. 
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Conversely,  it  we  insist  upon  maintaining  a  minimum  mean  r^j  =  rvj  ■  ,  then 

.  >  rN,min 

u  ~i^r 

T^r 

which  is  a  condition  upon  the  sampled  volume  dimension  parallel  to  the  mean  velocity.  If  rN  mjn  =  0. 1  ra 

and  -  1  cm'1,  then  i.,  must  be  >  1  mm. 

ui  J 


Clearly  the  particle  s  position  is  not  known  within  this  volume  without  auxiliary  data  (i.e.,  simultaneous  v. 
measurement).  This  is  of  little  concern  for  vorticity  statistics,  when  the  mean  flow  parameters  change  onl 
slightly  in  the  streamwise  direction  and  the  Lagrangian  vorticity  fluctuations  are  slow.  This  is  problematic 
however,  in  transforming  from  the  Lagrangian  to  the  Eulenan  viewpoint  where  the  spatial  vorticity  structu 
convecting  through  the  sampled  volume  is  sought  based  upon  the  random  time  senes.  The  minimum  com 
distance  is  typically  several  (6  to  8)  particle  diameters  at  high  particle  loading  (Subsection  2.4). 

The  temporal  resolution  of  the  VOP  is  a  more  difficult  issue.  Independent  of  flow  fluctuations,  the  instru: 
uncertainty  in  the  vorticity  components,  5coj.  can  be  expressed  as  an  uncertainty  relation 

5«i  td  =  Cj 

where  tq,  detector  transit  time,  is  the  time  for  which  the  particle  is  observed  and  C,  is  a  constant  (or  func 
which  depends  upon  expenmental  parameters.  The  longer  the  parttcle  reflection  can  be  observed  the  bette 
(he  precision  of  the  vorticity  measurement.  A  changing  vorticity,  however,  imposes  an  upper  limit  on  rjy 
In  a  turbulent  flow,  an  eddy  is  not  usually  expected  to  persist  for  more  than  a  single  rotation,  i.e.,  the  vor 

—  —  rN 

vector,  w,  changes  by  Ofui)  m  the  time  ~2xlw.  The  time  for  which  the  particle  is  observed  is  rjy—  — 


5w  _ 
ui  2t 
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independent  ot  instrumental  noise.  The  error  associated  with  being  always  arbitrarily  close  to  short-lived  fine 
scales  ot  turbulence  (the  'support'  ot  the  dissipation)  produces  an  additional  broadening  in  turbulent  PDFs 
which  does  not  exist  in  laminar  distributions.  Interestingly  therefore,  even  \  ~  distributions  should  have  strong 
turbulence  signatures,  since  non-constant  vorticity  (vortex  stretching)  degrades  fit  quality  (see  Subsection  5.2.2). 
Fhe  implications  ot  Eq.  (29)  deserve  further  study. 


The  other  temporal  resolution  issue  is  data  rate.  Specifically,  what  is  the  peak  mean  data  rate?  Since 
overlapping  traiectones  can  be  recognized  and  rejected,  events  of  mean  duration  tq,  can  be  crowded  together 

with  mean  time  between  events  rg,  until  overlap  eliminates  more  events  per  unit  time  than  are  added  by 
decreasing  rg  For  random  data,  a  Poisson  distribution  describes  the  probability  that  an  event  will  occur  at  time 
t .  given  an  event  at  time  t  =  u.  and  a  mean  time  between  events,  rg,  or  event  rate  Rg  ( =  l/r^), 

-r/Tp 

PE<t:rE>  =  e 

For  events  of  mean  duration  tq.  the  (average)  total  probability  of  overlap  is  the  probability  that  the  next  event 
begins  within  time  tq. 

7D 

C  -t/rp  ,  (  -  m/rpi 

^overlapt^-T)!  =  Ti  I  e  dt  =  1,1 -e  / 

But  tor  each  overlap,  two  events  are  lost.  Thus  the  net  data  rate  is  the  mean  event  rate.  Rg,  less  the  overlap 
rate 

RDATA  =  RE  -  2RE(l  -  e  r°/TE)  (valid  for 


The  peak  mean  data  rate  is  that  for  which  Rqata  is  a  maximum: 


dRDATA 

drE 


at  RET)  =  TD,TE  *  0  3 1 
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The  mean  time  between  events  must  be  about  three  times  the  mean  event  duration  to  achieve  the  highest 
rate  uhich  Rpg^K  -  0.46RE.  The  peak  mean  data  rate  for  -  0.5  ms  is  -300  Hz.  At  rate. 

half  of  the  events  must  be  discarded  due  to  trajectory  overlap.  Under  these  circumstances,  the  importan 
reliable  error  recognition  algorithms  is  evident.  The  peak  instantaneous  data  rate  is  considerably  higher, 
only  by  the  minimum  spacing  of  the  shortest  events.  These  closely  space  events,  however,  are  extreme! 
and  very  short  time  autocorrelation  data  is  sparse  for  most  applications. 
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4  water  compatible  vop  particles 


fhe  description  or  the  VOP  to  this  point  has  been  quite  general  in  regard  to  fluid/particle  systems,  stressing  the 
need  for  refractive  index  matching  and  the  advantages  of  neutral  buoyancy,  etc.  For  water  as  the  working  fluid, 
however,  these  simple  prescriptions  constituted  a  real  technical  challenge:  the  list  of  possible  particle  materials 
satisfying  these  needs  is  very  brief  (see  Table  2  below). 

In  practical  terms,  water-compatibility  is  essential  for  the  development  of  the  VOP  into  a  widely  used  research 
instrument.  Most  fluid  dynamics  laboratories  around  the  world  have  either  water  tunnels  or  wind  tunnels. 
Oreanic  fluids  cannot  be  senouslv  considered  tor  large-scale  applications  (with  the  notable  exception  of  the 
Pennsylvania  State  Universitv  glycerine  tunnel  in  which  the  PMMA  VOP  particles  could  work). 

To  solve  this  problem,  a  new  class  of  vorticity  probe  particle  made  of  Acrylamide  gel  which  provides  excellent 
index  and  density  matching  was  invented.  The  quality  of  the  match  is  due  primarily  to  the  fact  that  gel  is 
comprised  of  90  to  95  percent  water.  After  a  brief  review  of  particle  manufacturing  methods  and  materials,  the 
new  Acrylamide  process  will  be  described  and  VOP  water  calibration  data  in  laminar  Couette  flow  will  be 
presented. 


Table  2.  Candidate  Materials 


Material 

nD 

Status 

Fluorocarbon  solids 

Pentadecafluorooctvl 
Methacrylate  (PDFOMA) 

-  1.355* 

1.8* 

Mirrors  incorporated,  low  yield: 

VOP  tested 

Fluorocarbon  liquids 

Perflonnated  polyether 
Perfluorophenanthrene 

-1.31* 

1.334 

2.0* 

2.0 

Mirrors  not  yet  incorporated 

Gels 

Agar  (polysachamdes) 
Acrylamide  (5  to  10%) 

7 

— 1.34* 

7 

1.04* 

Not  yet  tested 

Best  probe  particles  to  date 

"Values  are  concentration,  or  molecular  weight  dependent. 
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4  I  Review  or  Particle  .Manufacturing  Techniques 


For  a  description  or  the  onginal  VOP  particle  process,  refer  to  Fnsh  and  Webb  ( 1981).  The  technique  i 
hnetly  recounted  below  to  introduce  the  concepts  of  suspension  polymerization. 

4.1.1  Suspension  Polymerization 

All  types  of  VOP  particles  are  manufactured  bv  suspension  polymerization  techniques  involving  four  bas 
features: 

1 )  Liquid  monomer,  from  which  the  solid  polymer  is  made  in  the  presence  of  a  polymerization  inn 
dispersed  as  small  droplets  in  a  medium  in  which  it  is  insoluble. 

2)  The  medium,  either  aqueous  or  non-aqueous  depending  upon  the  type  of  monomer,  actively  or 
passively  inhibits  coalescence  of  the  dispersed  monomer  droplets. 

3)  The  surface  tension  between  the  phases  maintains  sphericity  of  the  droplets  during  the  penod  in 
the  material  cures. 

4)  Lead  carbonate  mirrors  dispersed  in  the  monomer  pnor  to  creating  the  suspension  remain  withn 
droplets. 

The  pearl  pigmentation  additive  ZTX-B  (formerly  a  product  of  the  Mearl  Corporation,  now  available  on 
foreign  suppliers),  containing  basic  lead  carbonate  mirrors  was  designed  by  the  manufacturer  for  use  in  . 
plastics.  It  contains  a  plasticizer,  dibutyl-pthalate  (DBT),  and  ketones  as  solvents  or  thinners.  This  slur 
readily  miscible  with  methyl  methacrylate,  the  monomer  of  polyimethyl  methacrylate)  (PMMA)  from  wi 
onginal  panicles  were  made.  The  target  concentration  of  mirrors  in  the  distilled  monomer  will  produce 
average  of  one  to  two  mirrors  per  30  jim  diameter  droplet.  It  is  evidently  important  for  the  mirrors  to  t 
well  dispersed  in  the  monomer  to  achieve  the  most  uniform  distribution  of  mirrors  in  the  suspension.  A 
quantity  of  benzoyl  peroxide  is  added  to  this  solution  (1  g/liter)  as  an  initiator.  A  suspension  of  small  n 
droplets  is  made  in  a  buffer  solution  by  vigorous  agitation  (a  blender  works  well). 
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The  bulfer  solution  consists  of  low  molecular  weight,  80  percent  hydrolyzed  poly(vmyi  alcohol)  (PVA)  at  10  to 
20  g  liter  in  distilled  water.  The  PVA  acts  as  an  active  dispersant  which  inhibits  thinning  of  the  buffer  solution 
layer  between  droplets  when  they  interact.  The  surfactant  action  of  the  PVA  also  aids  by  inhibiting  the 
formation  of  new  surface  by  creating  a  skin -like  layer  on  the  droplet.  Small  quantities  of  salts  are  added;  NaCl 
modifies  surface  tension,  and  reduces  the  tendency  for  electrostatic  forces  to  promote  coalescence,  and 
ammonium  thiocyanate.  NH4SCN.  inhibits  the  formation  of  emulsion  polymer.  This  suspension  is  heated  and 
stirred,  and  the  droplets  cure.  The  product  is  a  fine  powder  of  hard,  clear,  spherical  PMMA  particles  and  the 
mirrors  remain  inside  of  spheres  of  sufficient  size  to  accommodate  them.  The  desired  size  distribution  is 
obtained  by  repeated  gravitational  sedimentation. 

This  is  the  model  of  suspension  polymerization  that  has  guided  and  informed  the  development  of  water 
compatible  VOP  panicles  using  new  matenals.  The  first  approach  to  water  compatibility  used  the  fluonnated 
acrylic  plastic  poiy(pentadccatluorooctylmethacrylate)  or  PDFOMA.  Several  crucial  differences,  due  to  the 
troublesome  nature  of  fluorocarbon  chemistry,  limit  the  usefulness  of  the  PMMA  model  in  this  case.  The 
second  and  most  fruitful  avenue  employed  acrylamide  gel.  This  water  soluble  polymer  network  effectively  turns 
the  PMMA  model  inside-out,  with  a  non-aqueous  buffer  and  an  aqueous  polymer  solution  as  the  dispersed 
droplet  phase.  After  summarizing  the  fluid/panicle  systems  that  have  been  used  or  are  still  in  development,  a 
description  of  the  water  compatible  particles  is  given  in  Subsection  4.2. 

4.1.2  Summary  of  Fluid/Particle  Systems 

Table  3  lists  the  fluid/particle  systems  conceived  to  date,  along  with  the  relevant  physical  properties  of 
the  matenals  and  the  particle  dynamic  response.  For  PMMA  probe  particles,  both  p-cymene  and  DBT  are 
routinely  used  as  working  fluids  in  small  scale  experiments  and  in  our  precision  Couette  flow  apparatus 
described  in  Chapter  5.  P-cymene  is  a  light,  strong-smelling,  aromatic  hydrocarbon.  Its  toxicity  is  low,  but  it 
is  quite  flammable.  OBT,  as  mentioned  above,  is  a  commercial  plasticizer  with  30  times  the  kinematic  viscosity 
of  p-cymene.  These  two  fluids  provide  a  convenient  means  of  spanning  a  large  Reynolds  number  range  with  a 
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Table  3.  Comparison  of  Particle  Response  Characteristics  tor 
Particle/Fluid  Combinations 


Densities 

Wf 

Refractive 

Indices 

w 

Translational 
Response  Tu 
(ms) 

Rotational 
Response 
Tr  (ms) 

F 

Vis 

*Y 

30  um 

Diameter 

Particles 

t^2 

t*2 

PMMA 

P.Cvmeme 

1. 2,0.9 

1.49,1.49 

0.18 

0.04 

0.0t 

pmma 

Dibutylpthalate 

1.2,1.04 

1.49,1.49 

0.006 

0.001 

0.1 

PMMA 

60  wt  %  Z  1  in  water 

1. 2.2.0 

1.49.1.49 

0.06 

0.025 

0.0 

PDFOMA 

Water 

1.75.1.0 

1.355.1.33 

0.1 

0.026 

0.0 

Acrylamide 
gel.  8  wtS 

Water 

1.04.1.0 

1.34,1.33 

0.08 

0.016 

0.0 

LLNL 

Aerogel* 

Air 

0.005,0.00 

1 

1.00,1.00 

0.03 

0.008 

0.1 

Polv(methlmethacrylate)  =  PMMA 

Polv(peniadecafluoroctylmethacrylate) 

•Not  yet  attempted 

=  PDFOMA 

single  flow  apparatus.  Aqueous  ZnI->  solutions  are  not  currently  used  because  they  are  quite  dense  and 
corrosive. 


New  silica  aerogel  technology  from  Lawrence  Livermore  National  Laboratories  has  made  anair-compatih 
a  possibility.  These  "frozen  smoke’  aerogels  have  been  produced  with  densities  of  -0.004  g/cc  and  rel 
index  near  unity.  A  VOP  probe  particle  of  this  type  could  be  used  in  wind  tunnels.  The  potential  of 
manufacturing  such  aerogel  spheres  is  under  study. 

4.2  Water  Compatible  Particles 

Two  successful  approaches  to  water-compatibility  ,  the  fluorinated  acrylic  plastic  and  acrylamide  gel  met 


are  described  below. 


4  2. 1  PDFOMA  Particles 


PDFOMA  is  a  perfluoroalkvl  ester  of  methacrylic  acid  which  is  used  in  the  contact  lens  industry,  co- 
polymenzed  with  other  materials.  It  is  quite  expensive,  about  $3  per  gram.  Its  polymerization  chemistry  is 
essentially  the  same  as  other  acrylics,  except  for  its  extreme  sensitivity  to  the  presence  of  oxygen  which  inhibits 
polymerization.  Azobis-isobutyronitrtle  (AiBN)  is  the  preferred  initiator.  The  first  diffirulty  encountered  with 
PDFOMA  is  that  basic  lead  carbonate  mirrors  do  not  disperse  readily  the  monomer  The  suspending  agents 
m  ZTX-B  (or  P80-H)  are  not  compatible  with  the  monomer.  Careful  multiple  washings  in  alcohol  are  required 
to  remove  the  suspending  materials  (e.g.,  even  a  mono-layer  of  DBT  on  mirror  surfaces  is  apparently  enough  to 
promoted  clumping  of  mirrors  in  PDFOMA).  Final  dilution  must  be  one  part  in  -  10®.  Even  so,  the  strongly 
hydrophobic  perfluonnated  chains  are  not  compatible  with  the  basic  lead  carbonate  mirrors.  In  order  to  render 
the  mirrors  equally  hydrophobic  they  are  'teflon-coated'  by  treatment  with  a  dilute  perfluoroalkvl  tnethoxystlane 
volution.  The  adhesion  of  the  siianol  groups  to  the  mirrors  is  still  rather  poor,  but  evidently  enough  to  give 
good  dispersion  in  the  monomer.  This  monomer  solution  must  be  de-gassed  to  eliminate  oxygen  and  excess 
alcohol.  The  buffer  solution  must  also  be  modified.  The  hydrophobic  droplets  und  to  cluster  quickly.  Increased 
concentrations  of  higher  molecular  weight  PVA  are  required.  NH4SCN  is  no  longer  needed  but  higher  NaCI 
concentrations  are  used.  Pure  nitrogen  is  bubbled  through  the  buffer  while  heating  and  stirring  prior  to  and 
Junng  polymerization  to  prevent  oxygen  poisoning. 

A  transmission/ re  flection  photomicrograph  of  the  resulting  product  is  shown  in  Figure  20.  The  boundaries  of  the 
spherical  poly(PDFOMA)  particles  are  clearly  visible.  These  particles  are  from  10  to  50  jim  in  diameter  and 
only  25  percent  contain  mirrors  which  are  barely  visible  m  the  photograph,  home  mirrors  can  be  seen  nearly 
edge-on.  while  others  show  clear  reflections.  Some  mirrors  exhibit  a  poorly  understood  surface  contamination 
which  reduces  their  reflectivity. 

PoIy(PDFOMA)  particles  have  been  tested  in  water.  The  reflections  are  adequate,  but  the  data  rejection  rate  is 
high.  Further,  these  particles  are  too  dense  to  be  used  in  many  applications,  particularly  when  larger 
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Transmission/Reflectton  Photomicrograph  of  PDFOMA  VOP  Particles 


particles/ mirrors  are  needed.  The  success  of  the  acrylamide  gel  process,  however,  obviates  the  need  to 
PDFOMA  process  further. 

-I  2.2  Acrylamide  Gel  Particles 

Acrylamide  gel  is  a  visco-elastic  material  that  is  commonly  used  as  an  electrophoresis  medium  in  bioch> 
research  and  biotechnology  industries.  The  gel  is  made  by  polymerizing  acrylamide  with  small  quantit; 
acrylamide  in  a  dilute  aqueous  solution,  resulting  in  cross-linked  polymer  network  which  forms  a  rubbe 
that  is  90  to  95  percent  water  by  weight.  The  polymerization  reaction  requires  the  addition  of  an  initia 
ammonium  persulfate,  and  proceeds  rapidly  when  a  catalyst,  tetramethylethylenediamine  (TEMED)  is  a 
the  solution.  The  mechanical  and  optical  properties  of  the  resulting  gel  material  are  controllable  via  pol 
concentration  and  degree  of  cross- linking,  but  typically  has  an  index  of  refraction  of  —  1.34  and  densir 
1.04  g/cc.  Satisfactory  VOP  particle  characteristics  have  been  achieved  using  an  8  percent  by  weight  a 
solution  ot  acrylamide/bis-acrylatmde  at  a  50:1  ratio. 


To  the  author's  knowledge,  suspension  polymerization  (in  a  non-aqueous  medium)  has  never  before  been  used  to 
manufacture  acrylamide  gel  spheres.  This  procedure  is  not  a  straightforward  extension  of  textbook  methods.  In 
particular,  difficulties  are  encountered  because  the  polymer  forms  rapidly  (often  within  seconds)  after  addition 
of  the  TEMED  catalyst  to  the  unpolymenzed  acrylamide  solution,  thereby  preventing  the  break-up  into  droplets. 
Furthermore,  the  mirrors  tend  to  clump  in  the  aqueous  solution  and  tend  to  migrate  to  the  boundary  between  the 
solution  and  the  suspending  medium. 

The  first  of  these  difficulties  is  solved  by  selecting  a  suspending  medium  with  special  properties.  "Thoroil",  a 
product  of  the  Kurt  J.  Lesker  Company,  is  chemically  inert,  is  immiscible  with  the  gel  ingredients,  has  low 
toxicity  and  is  inexpensive.  Most  importantly,  it  is  sufficiently  polar  to  dissolve  the  TEMED  catalyst.  This 
latter  property  is  critical  to  the  success  of  the  gel  particle  manufacturing  process  since  by  dissolving  the  catalyst 
directly  in  the  oil.  the  catalytic  action  is  inhibited  until  the  acrylamide  solution  is  dispersed  in  the  oil.  The 
catalyst  then  enters  the  spherical  droplets  by  diffusion  and  curing  begins.  The  problem  of  mirror  clumping  and 

migration  is  remedied  by  the  addition  of  high  molecular  weight  polyethylene  oxide  ( tS  0. 1  percent).  This  has 

the  effect  of  producing  a  highly  non-Newtonian  fluid  which  effectively  immobilizes  the  mirrors.  Acrylamide 
polymerization  is  also  inhibited  by  oxygen.  All  of  these  fluids,  including  the  suspension  oil.  must  be  thoroughly 
degassed.  A  more  detailed  discussion  of  acrylamide  particle  manufacture  is  included  as  Appendix  F. 

A  phase  contrast  photomicrograph  of  an  early  batch  of  acrylamide  gel  particles  is  shown  in  Figure  21.  The 
hexagonal  mirror  shape  can  clearly  be  seen.  The  index  match  is  so  close  that  without  phase  contrast,  the  gel 
sphere  boundaries  would  be  invisible.  Many  of  these  spheres  are  more  than  100  nm  in  diameter,  and  some 
mirror  clumping  is  visible.  Quality  control  and  continuous  processing  can  improve  these  particles  and  better 
control  their  size  distribution,  though  they  have  already  provided  high  quality  vorticity  data  in  water  without  the 
problems  associated  with  the  PDFOMA  process. 
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Figure  21 


V-503 


Phase  Contrast  Photomicrograph  of  Acrylamide  Gel  VOP  Particles 

4.2.3  Water  Calibration  Data 

The  acrylamide  gel  particles  were  first  tested  in  water  using  a  precision  Couette  flow  apparatus.  This  dt 
described  in  Chapter  5.  Suffice  it  here  to  say  that  the  particles  were  uniformly  distributed  in  a  laminar 
with  a  spatially  and  temporally  constant  vorticity  field  with  expected  value 

(w'x, u/y,u!’z)  =  (0,  0,  1.35  s  *) 

The  measured  vorticity  distributions  an-  presented  in  Figure  22.  The  wz  (spanwise)  distribution  Figure 
exhibits  a  strong  peak  at  1.25s  with  a  FWHM  -  10  percent.  It  should  be  emphasized  that  this  is  a  call! 
Small  variations  in  the  positions  of  optical  elements  result  in  a  net  departure  from  the  nominal  angular 
magnification.  The  calibration  factor,  i.e.,  the  ratio  of  the  expected  and  measured  peak  values  quoted  a 
should  be  used  to  scale  all  data  acquired  in  the  same  optical  configuration  thereafter.  Figure  22(b)  is  th 
(streamwise)  whose  mean  value  is  near  zero  as  expected,  with  —5  percent  FWHM.  The  ny  (normal) 
distribution  is  many  times  broader  than  the  others.  This  is  the  anticipated  consequence  of  the  VOP  geo 
and  finite  instrumental  precision  examined  in  Chapter  3. 
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Figure  22 

1  Jiininar  Flow  Histograms  (R/Rc  =  0.5)  Water  Calibration  Data 

Several  aspects  of  these  calibration  data  require  some  discussion.  In  particular,  the  existence  of  a  broad 
secondary  peak  at  -w'z/2  must  be  accounted  for.  Examination  of  the  individual  events  comprising  the 
secondary  peak  appear  to  be  valid  data.  Indeed  they  are.  This  peak  is  the  result  of  stable  particle  pairing  at 
low  shear  rate  (u'z)  due  to  weak  surface  interaction.  Also  the  widths  of  these  distribution  should  be  derivable 
from  the  component  error  distributions  and  the  resolution  function  arguments  of  Chapter  3. 


From  the  expression  of  Jeffrey  (1922),  the  angular  velocity  Q,  of  an  ellipsoidal  particle  in  a  uniform  shear,  u0, 
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1  *  «,  sm-7nu.<0t 


■'.here  a„.  B„.  anil  7^  are  coefficients  depending  upon  the  aspect  ratio,  17,  of  the  particle.  Paired  sphenc 
particles  ot  radii  a  , .  and  ai  are  taken  to  behave  as  ellipsoids  of  aspect  ratio 

aj  *  a-> 

V  =  - —  ■ 

max(a|,a->) 


For  a^  =  a->.  rj  =  2.  The  expression  for  p^fO),  the  probability  density  function  (PDF)  of  observing  the 
rate  0.  is 


PtjIO)  =  c 


an 

IT 


Differentiating  Eq.  (33)  and  eliminating  t,  gives 


“ 

P,(0)  = 

<v„«2 

Lsl  -I]  - 

1  BO  1)  ] 

a  1  |2 

BO  Bj 

in  the  interval,  <  Q  <  .  The  PDF  for  a  single  sphencal  particle  reduces  to  po(0)  =  5(0  -  0r 

the  total  PDF  is 


p(0)  =  ( 1  -0po(tl)  *  tp^(0) 


where  f  is  the  fraction  of  paired  particles  observed.  For  the  calibration  data  shown  here,  f  *  0.4. 

If  the  steady  slate  is  taken  to  be  paired  fraction  f.  at  the  shear  u >g,  then  a  calculation  of  pairing  (i.e..  col 
rate  should  allow  an  estimate  of  pair  lifetimes  in  the  shear.  If  collisions  are  assumed  to  take  place  for  a 

of  impact  parameter  b  $  d  ( =  2a),  where  d  is  the  particle  diameter  (Subsection  2.3.4)  then  the  total  co; 
frequency  is 
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d  _ 

h  =  j^8nVd^  -  u>0b  db 

=  ino)0d3 

ao)0 

r 

w here  a  is  the  volume  traction  of  particles  in  the  flow.  If  the  pairs  have  a  mean  lifetime  rp;ur,  then  the 
equilibrium  pair  fraction  is 

f  =  Qfai°rPair  (34) 

T 

For  the  low  shear  calibration  data,  a  =  10'3,  wQ  -  1  s'*,  and  f  -  0.4,  so  rpajr  ~  lO^s,  suggesting  that  the 
pairs  survive  indefinitely  at  this  shear  rate. 


This  result  is  somewhat  surprising.  Figure  23  is  a  map  of  the  relative  motion  of  sphere  centers  in  the  plane  of  a 
uniform  shear,  taken  from  Batcheior  and  Green  (1972a).  Open  sphere  trajectories  are  those  for  which 
(b/a)“  >  0  in  the  preceding  analysis  where  a  is  the  particle  radius,  but  for  (b/a)“  <  0  (imaginary  impact 


Figure  23 

Interaction  of  Two  Freely-Moving  Spheres  in  a  Linear  Shear  Flow  (from  Batchelor  &  Green,  1972a) 
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parameter,  b!)  sphere  traiectones  are  closed,  i.e.,  the  pairs  are  stably  bound.  The  imaginary  impact  par. 
'imply  interpreted  as  the  inaccessability  of  the  (b/ap  <  0  region  to  unperturbed  trajectories.  Apparent: 
brownian  motion,  triplet  interaction,  or  flow  field  fluctuations  (perhaps  also  sphere  imperfections,  comp 
'urtace  interactions)  can  transport  spheres  in  relative  motion  across  the  (b/a)“  =  0  boundary  (in  both 
directions).  Pair  lifetimes  can  be  quite  large  tn  laminar  shear  flows.  After  a  sufficient  period  of  time 
'hearing  tlow.  starting  from  some  initial  distribution,  it  appears  that  a  steady  supply  of  single  particles  i 
available  only  because  triplets,  etc.,  are  unstable  in  the  flow  and  the  odd  particles  is  always  stripped  lo< 
unstable  pairs  (e.g.,  in  a  turbulent  shear  flow)  the  lifetime  approaches  one  turnaround  time  l/ui0,  and  th 
equilibrium  pair  concentration  reduces  to  —  air.  Evidence  for  pair  disruption  by  turbulence  is  seen  in 
Chapter  5. 

The  component  error  (ffj)  distributions  are  shown  in  Figure  24.  As  was  argued  in  Subsection  3.4.1.  the 
seen  to  be  distributed  approximately  as  l/ffj“.  From  estimates  of  <rmm,  and  <rmax  the  expected  resolutii 
functions  can  computed  and  compared  to  the  measured  calibration  PDFs.  In  the  special  case  of  wz,  the 
resolution  function  is  convolved  with  Eq.  (33)  to  include  the  pair  contribution.  The  mean  aspect  ratio  i 
to  be  rj  =  1.6  and  f  =  0.4.  These  PDFs  are  shown  in  Figure  25.  The  essential  features  of  these  data 
apparently  accounted  for  in  this  analysis.  There  appears  to  be  an  excess  noise  not  seen  in  the  cr-  distnb 
however,  which  is  particularly  apparent  in  the  wings  of  the  wv  PDF.  Further  study  is  required. 
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Figure  24  Figure  25 


Measured  Error  Distributions  Corresponding  Resolution  Functions  Calculated  from  the  Measured 

to  the  Histograms  of  Figure  22.  Error  Distributions  (Compare  with  Figure  22) 
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5.  VORTICITY  MEASUREMENTS  IN  COLETTE  FLOW 


"jvlor-Couette  tlow.  the  flow  between  concentric  rotating  cylindrical  surfaces,  continues  to  be  an  impor 
iaboraiorv  tor  the  study  of  the  transition  to  turbulence  and  of  large  scale  coherent  vortical  structures.  V< 
measurements  have  been  performed  in  six  distinct  Couette  tlow  regimes,  which  are  enumerated  in 
Subsection  5.1  1.  Calibration  data  is  acquired  in  the  laminar  tlow  regimes.  Evidence  of  coherent  vortu 
structure  is  observed  in  the  outer  boundary  rotation  turbulence  data.  The  statistics  of  these  turbulent  tlo 
Mtnilanties  to  standard  turbulent  boundary  layers.  Inner  boundary  rotation  flows  show  the  characteristic 
patterns  of  the  ‘spectral  evolution"  to  turbulence  (Coles,  1965)  as  verified  by  tlow  visualization  at  Revr 
number  corresponding  to  the  measurements.  Spatial  and  temporal  periodicity  of  the  components  of  the 
ire  measured,  which  are  compared  with  numencal  simulations. 

5.  1  Couette  Flow 

In  selecting  a  tlow  field  in  which  to  verify  the  operation  of  the  VOP  both  in  laminar  and  turbulent  flow 
conditions,  the  following  criteria  were  applied:  1)  the  flow  apparatus  must  be  useable  with  all  of  the  VC 
working  fluids  used  previously,  culminating  with  water;  2)  it  should  be  a  device  of  small  to  moderate  sc 
3)  it  must  possess  well  defined  laminar  regimes  where  the  vorticity  distribution  is  well  known;  4)  it  mu 
range  ot  accessible  Reynolds  number  sufficient  to  reach  fully  developed  turbulence  with  statistics  comp; 
other  standard  flows;  and  5)  the  potential  to  display  coherent  vortical  structures.  It  is  difficult  to  concei 
tlow  which  meets  these  criteria  more  precisely  than  Couette  tlow.  A  benchtop  precision  Couette  flow  d; 
be  made  to  function  with  virtually  any  fluid  in  relatively  small  volumes.  Sub-cntical  flows  posses  a  cor 
vorticity  vector  everywhere  in  the  flow.  Depending  upon  which  boundary  surface  is  rotating,  Couette  f 
vorticity  field  exhibit  astonishing  variety  and  richness. 

5.1.1  Couette  Flow  Regimes 

As  mentioned  above,  there  is  great  variety  in  Couette  flow  behavior  which  can  be  generally  divided  inti 
categones:  flows  dominated  v  >  the  rotation  of  the  inner  (cylindrical)  boundary  or  those  dominated  by  i 
boundary  rotation. 
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Let  the  radius  and  angular  velocity  of  the  inner  and  outer  cylindrical  boundaries  be  denoted  r- and  rQ,fl0 
respectively.  The  boundary  conditions  on  the  velocity  field  are  Ug(r  =  rj)  =  0^  and  Uj(r  =  rQ)  =  Q0r0,  the 
time-independent  < sub-crnical)  solution  ot  the  Navier-Stokes  equations  gives  (Landau  and  Lifshit.  1959) 


Ur  =  Uz  =  0 


Ug  =  ar*_ 


where  a  = 


ro°o  '  r,'n. 


and  b  = 


!°i  '  VVo 


From  the  equation  for  the  curl  operator  in  cylmdncal  coordmates,  the  only  non-zero  component  or  the  vorticity 

is 

a(rUg) 
dr 

=  2a 

Thus  the  vorticity  vector  has  a  constant  magnitude  and  direction  everywhere  in  flow  field.  Note  that  when 

*)  *> 

fl0r“  -  0 , r “  =0.  then  a  =  0  and  the  flow  is  lrrotational.  This  is  a  special  point.  The  necessary  and 
sufficient  condition  tor  the  stability  of  the  (tnvisctd)  flow  field  is  the  Rayleigh  criterion 

^{rUg(r))2  >  0 

w  hich  reduces  to 

Vo  >  v2 


For  the  case  of  outer  boundary  rotation  (Qo=0)  the  flow  is  stable.  For  inner  boundary  rotation  (0o=0)  the 
flow  is  unstable.  For  finite  viscosity,  there  is  a  critical  Reynolds  number  for  inner  cylindrical  rotation  marking 
the  onset  of  more  complex  flows. 


iiiich  cylinder  can  be  thought  ol  as  possessing  its  own  Reynolds  number  based  upon  boundary  surface  spec. 


and  gap.  g  =  r  ,-r| 


R0  = 


^orolro  ril 


.  R,  = 


Vi<ro  -r.) 


the  heicht  of  the  fluid  annulus  is  usually  much  greater  than  g.,  i.e.,  h/g  £  30.  A  How  is  generally  desert! 
hv  particular  type  of  coordinates  (R0,R,)  on  a  map  of  versus  RQ.  One  interesting  aspect  of  Couette  flo\ 
that  the  tlow  at  (Rq.R,)  is  not  necessanly  unique.  Several  different  stable  states  may  be  reached  depending 
upon  how  (  R0.R.)  is  approached. 


\  map  ot  R ,  \ersus  R.  from  Busse.  Gollub,  Maslowe.  and  Swinney  ( 19851  is  reproduced  below  m 
Future  26(a).  The  outer  boundary  rotation  dominated  flow  map  is  indicated  in  Figure  26(b)  by  greatly 
extending  the  Rn  axis.  Only  the  cases  of  R0  =  0,  R^  >  0  and  R0  *  0.  1^=0  are  considered.  Moving  up 
R  axis  from  0.  at  R,  =  Rc  the  Couette  flow  becomes  unstable  and  secondary  flows  known  a  Taylor  vortu 
appear.  These  vortices  are  stationary  counter-rotating  vortex  pairs  sucked  up  and  down  the  span  (see  Figi 
in  Section  5.3.1).  At  R,/Rc  -  1.2  the  secondary  tlow  itself  becomes  unsuble  to  azimuthal  disturbance,  x 
flow  selects  one  of  several  possible  periodic  sutes  known  as  wavy-vortex  tlow  (Coles,  1965).  This  sute 
characterized  hv  traveling  azimuthal  waves  superposed  (non-linearlv)  on  a  Taylor  vortices.  At  R^/RQ  -  1 
incommensurate  frequencies  appear  in  the  wavy-vortex  spectra  ( Fenstermacher.  Swinney.  and  Gollub.  197 
giving  rise  to  quasi-  periodic  flows.  With  increasing  Rj,  broadband  background  fluctuations  begin  to  appe 
finally  overwhelming  the  periodicity  which  is  no  longer  detecuble  at  R,/Rc  23.  The  basic  Taylor  vorte 
pattern  still  persists  with  fine  turbulent  structure  everywhere  tn  the  flow.  This  process  is  referred  to  as  tlv 
'spectral  evolution"  to  turbulence  (Coles.  1965). 

The  behavior  along  the  R0  axis  is  distinctly  different.  For  R^  =  0,  is  suble  and  persists  in  the  giving  the 
expected  Couette  value 
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Stt  ABOVE 


n  Huk 


(b)  Note  That  Onset  of  Fully  Developed  Turbulence  is  Dependent  on  Flow  Conditioning 

Figure  26 

Couette  Flow  Map.  VOP  Measurements  at  • 


20n 


1  - 


(35) 


until  wall-bounded  shear  flow  instabilities  akin  the  ordinary  boundary  lay  type  initiates  a  catastrophic'  transition 
to  turbulence.  This  turbulence  is  initiated  at  isolated  spots  in  the  boundary  lay  which  grow,  dominate  the  flow, 
and  then  disappear,  with  a  pulsating  sort  of  regularity.  In  this  respect,  turbulent  Couette  flow  is  similar  to  pipe 
flow,  exhibiting  tntenmttency  with  increasing  Reynolds  number,  until  the  intermittency  factor,  ,  approaches 
unity  when  the  flow  is  fully  turbulent. 


83 


5.1.2  Ciuelte  How  Apparatus 


The  Couette  tlow  device.  designed  and  constructed  for  this  program,  is  shown  below  in  Figure  27.  The 
relevant  parameters  are  specific  in  Table  4.  This  design  is  unusual  n  several  respects.  First,  it  is  actuall 
Couette  flows  rather  than  one.  as  seen  in  the  figure.  The  middle  cantilevered  cylinder  rotates  freely  abou 
central  cylindrical  spindle.  A  fixed  outer  cylinder  was  required  on  which  flats  could  be  cut  or  molded  for 
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Figure  27 

Precision  Couette  Flow  Apparatus  (timing  belt  drive) 


Table  4.  Couette  Flow  Apparatus  Parameters 
Inner  Fluid  Annulus 

rQ  =  4.477  cm,  r^  =  3.944  cm,  rj/rQ  =  0.881 
G  =  0.559  cm,  h/g  =  39 

Outer  Fluid  Annulus 

rQ  =  4.738  cm,  r,  =  5.127  cm,  rj/rQ  =  0.894 
G  =  0.611  cm,  h/g  =  36 
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mounting  of  the  output  coupling  optics  tor  the  VOP  described  in  Subsection  3. 1.2.  The  rationale  tor  this 
approach  denves  from  the  need  to  accommodate  the  fixed  outer  cylinder  constraint,  and  the  desire  to  study  both 
t  the  basic  flow  tvpes  described  in  the  previous  section  with  a  single  apparatus.  As  only  the  central  cylinder 
:  -tales,  the  outer  fluid  annulus  constitutes  an  inner  boundary  rotating  flow,  and  the  inner  fluid  annulus  gives  an 
uter  boundarv  rotating  tlow 

The  above  design  provides  this  flexibility  at  low  cost  since  no  elaborate  seals  or  bearings  need  be  immersed  in 
ihe  working  fluid.  This  last  feature  contributes  much  to  the  overall  advantage  ot  this  design  since  three 
different  working  fluids  are  used  which  no  single  set  of  seals  could  accommodate. 

The  fixed  central  spindle  and  base  plate  were  machined  from  stainless  steel  as  is  the  routing  chuck  which  ndes 
■n  sets  of  thrust  and  roller  bearings  pressed  onto  the  spindle.  The  plexiglas  cylinders  were  precision  machined 
and  polished.  The  outer  fixed  cylinder  is  mounted  on  a  sutnless  steel  base  plate  and  the  central  routing 
cylinder  is  mounted  in  the  chuck  with  adjusting  screws  for  alignment.  The  routing  cylinder  is  aligned  and 
balanced  in  place.  Baffles  cut  from  teflon  suppresses  irregular  end  disturbances.  Toul  runout  due  to  machining 
and  alignment  flaws  is  less  than  two  thousandths  of  an  inch  up  and  down  the  full  span  ot  the  annuli.  The 
central  cylinder  is  turned  with  a  variable  speed  synchronous  AC  motor  with  a  frequency  variable  synthesized 
three-phase  controller.  Timing  belts  and  pulleys  are  used.  Cylinder  roution  rates  are  stable  against  long  term 
Unit  to  -0  01  percent  < i.e. .  drifting  by  no  more  than  -2tt  in  10^  routions).  Periodic  speed  fluctuations  are 
observed  during  a  single  roution  of  -0.1  percent.  The  Reynolds  number  ranges  accessible  to  this  apparatus 
using  three  working  fluids  is  summarized  below  in  Table  5. 
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translation  stage,  The  surface  at  the  stainless  steel  inner  cylinder  was  black -passivated  to  eliminate  specular 
reflections  form  its  surface. 

The  vorticitv  components  measured  in  this  geometry  and  their  relationship  to  the  flowtield  are  indicated  in  the 
figure.  Explicitly,  the  standard  VOP  primed  basis  is  related  to  the  cylindrical  coordinate  vorticity  components 

-'z-u'0'wr>  by 

uiz'  =  wz 

1  , 

2  '  0  r  s  (36) 

1  , 

w‘v  =  -=  (“fl  -  <*M  = 

t/2 

where  ois  and  <cn  are  the  components  in  and  normal  to  the  plane  of  the  principal  axis  of  the  rate-of-strain  tensor, 
respectively.  The  component  in  this  plane  should  thus  experience  the  most  rapii'  stretching  and  consequent 
intensification. 

5.2.1  Transition  to  Turbulence  in  P-cvmene 

Probability  density  functions  (PDFs)  of  two  vorticitv  components  for  the  outer  boundary  rotating  flow 
are  presented  in  Figure  29(a)  through  (d).  Starting  at  1  Hz  rotation  rate  (a)  (cylinder  surface  speed  of  28  cm/s) 
and  increasing  by  factors  of  2  up  to  case  (d)  at  8  Hz  (surface  speed  2.25  m/s)  where  the  boundary  layer  appear 
to  be  fully  turbulent.  (Some  residual  intermittency  may  persist  at  this  Reynolds  number  -25,000  as  indicated  in 
Figure  26(b)  but  a  laminar  peak  is  difficult  to  distinguish  against  the  background  turbulence.)  Figure  30(a) 
through  (c)  are  flow  visualization  photographs  corresponding  approximately  to  the  1.  4,  and  8  Hz  cases.  The 
How  visualization  particles  have  been  introduced  into  the  inner  fluid  annulus  only.  Figure  30(a)  is  the  purely 
laminar  case.  Figure  30(b)  show  the  formation  of  fine-grained  turbulent  zones  which  migrate  through  the 
flowtield.  In  Figure  30(c)  the  turbulence  has  essentially  filled  the  inner  fluid  annulus. 
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Outer  Cylinder  Rotatom  f  low  Visualization 


In  Figure  29.  the  PDF  vorticity  axes  have  been  scaled  exactly  with  the  increasing  rotation  rate,  in  order  t 
cisual  comparison  will  convey  the  constancy  of  do >/u  up  to  the  strong  transition  to  turbulence  at  8  Hz.  T 
j.'  ,  (  =u'n)  distributions  are  not  included  since  the  instrumental  breadth  of  these  distributions  dominates  th< 
turbulent  fluctuations.  This  is  not  to  say,  however,  that  the  mean  values  and  second  and  higher  moments 
he  obtained  from  these  distributions.  The  question  will  be  examined  in  the  next  section  (Subsection  5.2.1 
helowi. 

There  is  some  evidence  at  low  speed  (1  and  2  Hz)  of  particle  pairing  as  was  noted  in  Chapter  4.  Howev 
4  Hz.  the  characteristic  side  lobes  have  been  reasonably  suppressed.  This  can  be  interpreted  as  good  evi 
that  weak  intermittent  turbulence  is  an  adequate  mechanism  for  disrupting  stable  pairs,  and  it  follows  tha: 
developed  (lows  are  unlikely  to  have  significant  pair  populations. 

5.2.2  Vorticity  Statistics  and  Profiles 

Table  6  contains  the  statistical  parameters  which  can  be  extracted  from  these  data  with  some  confidence, 
each  rotation  rate,  the  vorticity  scale,  the  Reynolds  number,  the  distance  from  the  inner  (stationary)  cylir 
which  the  measurements  were  made,  and  the  corresponding  laminar  vorticity  (from  Eq.  (35))  are  provide 
Note  that  the  number  of  measurements  in  each  PDF  is  not  large.  The  error  in  the  second  moments  is  of 

Aw./y  N  -  or  -5  to  10  percent  for  these  data.  Because  of  the  inefficient  nature  of  the  current  data  acqu 

and  reduction  procedure,  several  hours  and  30  to  40  data  sets  are  required  to  generate  a  few  histograms 
those  in  Figure  29,  though  the  data  is  collected  for  each  set  in  a  matter  of  seconds.  (  A  new  system,  cun 
the  design  stage,  will  acquire  data  and  process  vorticity  data  in  real  time.)  For  the  conditions  specified  1 
table  the  component  mean  vorticity  uj ,  rms  fluctuation  Aojj  and  an  estimate  of  the  rms  instrumental  broa 

based  on  the  error  distributions,  bwv  are  given.  Instrumental  broadening  scales  with  rotation  rate  as  expe 
The  widths  of  the  z,  s,  and  n  components  of  vorticity  can  be  seen  to  be  accounted  for  by  instrumental  m 
until  the  onset  of  intermittent  turbulence  at  4  Hz.  Turbulent  fluctuations  then  apparently  dominate  the  w 
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Table  0.  Vorticity  Statistics  tor  the  PDFs  of  Figure  29. 
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distribution  widths.  This  is  not  so  in  the  case  of  o;n,  however,  where  instrumental  noise  continues  to  dominate. 
\i  this  icsel  ot  statistics  iN).  estimates  ot  the  turbulent  contribution  to  the  sec  ind  moment 

l-Hif  -  ($“>nf 

are  unreliable  and  very  sensitive  to  the  estimate  of  owin.  The  same  may  be  said  for  higher  moments  ,i.e., 
skewness  and  flatness  factors,  for  the  PDFs  of  oiz  and  u  .  To  obtain  the  flatness  of  the  ui7  distribution  to 
<  10  percent  accuracy  would  require  two  orders  of  magnitude  increase  in  the  volume  of  data. 


Mean  and  fluctuating  vorticity  profile  data,  normal  to  the  inner  cylinder  surface,  was  also  obtained  only  at  a  few 
points  at  or  below  8  Hz.  Hie  normal  v  coordinate  will  always  be  referred  to  the  stationary  boundary.  The 
mean  vorticity  profile  and  rms  fluctuations  are  shown  in  Figure  31.  No  significant  variations  in  the  vorticity 
across  the  boundary  layer  were  observed  for  laminar  cases  (2  Hz  data  is  shown)  which  maintained  their  Couette 
flow  values  at  all  stations.  For  the  turbulence  case  (8  Hz)  the  rapid  increase  of  mean  vorticity  u> z  near  the  wall 

is  clear  indication  of  the  development  of  boundary  layer  type  turbulence.  Very  little  data  on  this  flow  is 
available  in  the  literature  (e.g.,  Taylor,  1936,  Coles,  1965  and  Reichardt,  1956).  Near  the  inner  wall,  however, 
the  flow  is  quasi-parallel  and  centripetal  accelerations  vanish.  The  flow  in  this  region  should  behave  much  like 
a  boundary  layer  in  turbulent  channel  flow. 
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"  Corrected  for  instrumental  and  gradient  broadening 
^From  Kim.  Moin.  and  Moser  (1987) 

♦  wx  fmm  Kastrinakis  and  Eckelmann  (1983) 

Figure  31 

Comparison  of  Measured  Couette  Row  Vorticity  Statistics  with 
Turbulent  Boundary  Layer  Computations  and  Data 

On  this  assumption,  the  law-of-ihe-wail  vorticity  profile  was  superposed  upon  a  uniform  background  von 
Assuming  a  symmetnc  vorticity  distribution  across  the  flow  (which  is  not  strictly  true,  e.g.,  Taylor,  1 9 3< 
discussion  of  Reynolds  number  dependence  below)  the  wz  profile  data  was  fit  to  the  proposed  vorticity 

distribution  with  the  constraints  that:  1)  the  estimated  (y)  integral  of  mean  z  vorticity  across  the  flow  wa.- 
to  the  known  outer  boundary  velocity  (2.25  m/s);  and  2)  that  y  was  non-dimensionalized  in  the  usual  wa; 


u 

V 

V 

where  uiq  is  the  mean  span  wise  vorticity  at  the  inner  wall  (y=0). 


The  best  fit  was  obtained  for  estimated  ui0  —  1000  s*^  and  the  mean  background  vorticity  ~360  s'*,  anc 
plotted  in  the  figure. 
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Hccau.se  or  the  unite  width  ot  the  laser  beam  (  150  «tm  or  about  6».  units  for  the  above  data),  regions  of 
unsiantial  spatial  vorticitv  gradients  produce  an  additional  broadening  (of  Auj2  in  this  case).  From  the  estimated 
dope  oi  ihe  mean  u;_  protile.  and  the  known  width  of  the  gaussian  beam  profile,  a  simple  subtraction  of  the 
.-timated  gradient  broauenini!  sutfices  to  correct  doz.  Based  upon  the  arguments  of  Chapter  3  and  the  measured 
.  rror  distributions  ot  the  vorticitv  components,  estimates  of  the  widths  of  the  vorticitv  resolution  functions  bo; 
•sere  derived,  ignoring  the  sanation  with  vorticitv.  The  true  rms  fluctuations  are  obtained  from  the  difference  of 
:he  squares  ot  the  measured  rms  values  and  the  instrumental  noise,  &u)j. 


fhe  rms  vorticitv  fluctuations,  when  corrected  tor  both  instrumental  and  vorticitv  gradient  broadening,  show 
.uantative  agreement  witn  turbulent  boundary  layer  data.  The  direct  numencal  simulations  of  Kim.  Moin.  and 

'  loser  1  1^87)  for  example,  trom  wmch  the  rms  vorticitv  curves  are  taken,  agree  with  the  measured  _ ,  at 

u’o 

least  to  edge  ot  the  butter  laser  iv  _  -40).  These  data  are  also  consistent  with  the  hot-wire  vorticitv  of 
dreamwise  vorticitv  ot  Kastnnakis  &  eckelmann.  and  the  9-wire  hot-wire  measurements  of  Balint.  el  al.Bevond 
this  point  the  Cuuette  values  are  likely  to  be  greater  due  to  the  large  mean  shear  away  from  the  wall.  The 
values  are  more  interesting.  Tiev  appear  to  be  larger  than  both  the  Au >ziui0  and  the  curves  predicted  by  the 


-imulalion.  This  should  not  be  surprising,  however,  since  o.\ 


=  _L(UJX 

t/2 


Uly) 


lies  in  the  plane  of  the 


principal  axis  oi  the  raie-ot -strain  tensor,  where  the  vorticitv  fluctuations  due  to  vortex  stretching  are  anticipated 
to  be  most  intense. 


These  comparisons  are  quite  sensitive  to  the  estimate  of  tug.  Further  the  argument  employed  to  tit  the  data  to  the 
Uw-ot-the-wali  plus  background  vorticity  is  tenuous.  With  increasing  Reynolds  number  the  Couette  flow 
vorticitv  protile  appears  to  become  increasingly  asymmetric  and  the  integrated  vorticitv  near  the  inner  wall  will 
bear  no  simple  relationship  to  the  outer  boundary  velocity.  A  more  detailed  protile  is  required.  More  extensive 
measurements  at  RQ  =32.000  are  presented  below. 
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3  2.3 


rnree-Component  Vorticity  Fluctuations  near  the  Wall 


\t  10  Hz  outer  boundary  rotation  rate,  corresponding  to  R0  =  32,000,  two  profiles  were  obtained:  the  I 
u.ith  counterclockwise  rotation  so  that  the  most  accurately  measured  components  are  oiz  and  ois  in  the  no 
,>t  Fieure  28.  and  the  second  with  the  rotation  in  the  opposite  sense,  so  that  in  the  same  vop  geometry  tl 
accurately  measured  components  were  w2  and  wn.  Thus  the  mean  and  fluctating  values  of  three  orthogo. 
vorticitv  components  could  be  accurately  obtained  without  greatly  increasing  the  number  of  data  points  r 


uz  Au>z  Accs  Ati)n 

The  profile  data  for  the  mean  spanwise  vorticity, - and  the  rms  fluctuations  -  ,  _  ,  _  in 

u>o  t*t0  ui0  ui0 

and  n  components  are  shown  in  Figure  32.  These  values  have  seen  normalized  by  the  measured  wall  vt 
=  v)75  s'1  The  beam  waist  for  this  data  was  -  100  fim  corresponding  to  4.6  y+  unit.  The  rms  fluL 


Figure  32 


Profiles  of  Mean  Spanwise  Vorticity  —  and  rms  Fluctuations 

wo 

A/oT  Acj*  Aojn 

_ _  ,  _  ,  _  for  10  Hz  Outer  Boundary  Rotation  (Rq  =  32,000) 

wo 


94 


i 


"jve  been  corrected  lor  irwrumenui  noise  and  gradient  broadening  as  described  above  !Tus  Iigure  combines 
:ne  z  data  tor  (he  clockwise  and  counterclockwise  rotation  which  were  round  lo  be  in  agreement  within 
•xpenmenia!  error  as  expected,  except  tor  the  change  of  algebraic  sign.  The  mean  values  of  the  other 
mponents  are  zero  within  experimental  error. 

rhe  mean  spanwise  vorticitv  profile  is  clearly  boundary  layer-iike.  The  scaling  is  no  longer  in  doubt  since 
measurements  of  mean  wall  vorticity  were  performed  deep  within  the  viscous  sublayer  (y  ^  -  2.3).  At  v  _ 

<  10.  the  tluctuating  vorticitv  statistics  show  good  agreement  with  the  simulatins  of  Kim.  et.al,  though  it  must 
oe  remembered  ur  and  <*■  represent  a  45  dea  rotation  of  the  at.,  uv  basis  in  the  simulations.  Neither  Awi  „ 

>  r  w  \  v  i.rms 

-.or  rms  -hould  tend  to  zero  at  the  wall,  as  is  required  ofAui.,  rms  since  all  x  and  z  spatial  derivatives  are 
Jentically  zero  at  the  wall  The  expected  persistance  of  strong  vorticity  fluctuations  beyond  y  _  -  40  due  to  the 
non-zero  background  shear,  is  observed.  The  vorticity  fluctuations  in  the  plane  of  the  principal  axis  of  the  rate 
t  strain  tensor.  Aw,.  rms.  are  profoundly  amplified  well  beyond  y_  -  40.  which  in  a  boundarv  layer  would  be 
(he  logarithmic  region.  The  could  be  taken  as  confirmation  of  the  numerical  observations  of  Rogers  and  Moin 
1987),  that  all  turbulent  shear  flows  should  exhibit  strong  vortical  concentrations  inclined  ai  -45°  to  the  mean 
flow  direction,  in  the  direction  of  the  principal  strain.  The  near  equality  of  spanwise  and  normal  fluctuations  is 
ilso  consistent  with  their  observation  in  direct  numerical  simulations,  that  vortex  filaments  (those  distinguishable 
:rom  the  background)  under  compression  tend  to  "buckle"  and  turn,  rather  than  weaken.  This  level  of  detail 
•vithin  the  sublaver  has  never  before  been  observed  due  to  serious  question  regarding  accuracy  (Taylor’s 
h  vpoihesis)  and  spatial  resolution. 

These  values  are  generally  larger  than  those  measured  by  in  channel  flow  by  Balint  et  al.  ( 1987)  using  the  9- 
wire  hot-wire  probe.  While  it  is  difficult  to  compare  these  flows  directly,  it  seems  clear  that  the  vorticity 
optical  probe  performs  admirably  under  the  worst  conditions,  i.e.,  high  Reynolds  number,  near  wall.  etc.  The 
questionable  nature  of  Taylor's  hypothesis  in  turbulent  shear  flows,  particularly  in  the  wail  region  (Piomelli  et 
ai.,  1989;  and  Champagne,  1978)  further  enhances  confidence  in  the  VOP  which  does  not  depend  upon  Taylor's 
hypothesis  in  any  way. 
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rhe  PDFs  up  to  now  have  been  treated  as  one-dimensional.  In  tact,  they  are  projections  of  the  measure 
.  ’rticitv  '.ector'  .  .to  selected  vorticity  basis  vectors.  Since  the  components  are  measured  simultaneouM 
possihl  •  .iew  the  projections  ot  vorticity  vectors  onto  selected  planes,  forming  two-dimensional  joint 
l"hese  then  represent  the  probability  of  finding  a  vortex  filament  or  tube  .or  just  random  background,- w 
particular  anitle  ot  inclination  in  the  plane.  In  a  manner  analogous  to  the  quadrant  analysis  of  the  Reyn* 
••tress.  W illmanh  and  Lu  ( 1972).  such  vorticity  maps  could  potentially  lead  to  a  more  quantitative  under 
of  coherent  structure,  in  particular,  "coherent  vorticity"  (Hussain,  1986).  For  the  turbulent  case  at  8  H. 
iwo-dimensionai  PDFs  are  presented  as  scatterplots  in  the  wz.u>s  plane  at  two  locations  in  the  flow. 
Figure  33(al  is  the  u iz.ws  scatterplot  at  v ^  =  12.  and  Figure  33(b)  at  y^  =  28.  Although  the  statistics 
poor,  the  distinctlv  different  character  of  these  planes  is  clear.  At  v  ^  =  12  the  transverse  fluctuation  u 
'trong  as  expected  and  large  vorticity  excursions  appear  to  be  forming  a  pattern,  which  is  not  quite  disc- 
At  y  ,  =28.  however,  even  at  this  level  of  statistics,  a  pattern  is  emerging  which  has  been  emphasized 
Jotted  triangle. 
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Figure  33 

u>v  uix  Scatterplots 
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‘  or  (he  ease  01  the  10  Hz  rotation  rate,  the  |oint  PDFs  ot  the  protile  data  are  shown  in  Figure  34(a)  and  34(b) 
as  seatterplots  in  the  plane  and  tc..ucn  plane  respectively.  Again,  a  pattern  emerges  which  is  not  evident 

trom  the  profiles  tn  Figure  3  1  Two  significant  features  stand  out  tn  these  joint  PDFs:  1 )  Aside  from  the  mean 
.'low  core  which  is  becoming  more  symmetrical,  there  is  a  distinct  tendency  for  large  streamwise  ( ojs)  excursion 
to  have  small  spanwise  excursion,  and  vice-versa.  This  suggests  the  symmetrical  bending  and  amplification  of 
’.  orticity  for  reasonable  traction  ot  the  sample  data.  This  is  consistent  with  the  counter-rotating  vortices  of  the 
hairpin  vortex  model,  particularly  in  the  <cz.ojs  plane  though  it  is  not  discemable  whether  the  vortices  occur  in 
>patially  contiguous  pairs:  and  21  There  is  also  the  smattering  of  very  intense  activity  scattered  about  the  uj  ,ojs 
plane  near  y,.  that  is  not  as  pronounced  elsewhere.  Little  can  be  made  of  this  at  this  stage  (e.g.,  in  terms  of 
filament  painng  as  in  Pumir  and  Sicgia.  19871  without  vastly  improved  statistics  over  a  much  larger  domain.  It 
appears,  however,  that  a  larce  contribution  to  the  vorticitv  fluctuations  is  associated  with  isolated,  intense 
'.orticity  concentrations.  Higher  moments  are  needed  at  higher  Reynolds  number.  The  VOP  may  be  uniquely 
'Uited  to  this  task. 


Other  auxiliary  data  types  might  prove  interesting  in  some  applications.  From  the  discussion  in  Subsection  3.4. 

•> 

it  will  be  recalled  that  the  goodness-of-fit  parameter  is  sensitive  to  angular  acceleration,  i.e..  __ . 

interpreting  this  as  an  error  in  the  vorticitv  measurement.  But  since  instrumental  errors  can  be  largely 
accounted  for  by  the  procedures  outlined  in  previous  chapters,  this  fluid  dynamically  driven  feature  should 
appear  as  a  residual  error  which  can  be  characterized.  Recall  that  under  the  action  of  a  straining  field  eij  an 
aligned  vortex  filament  will  intensify  as 


dui: 

-  OC  u  ■  ejj 

dt  1  11 


(37) 


This  should  appear  as  an  "acceleration  tail"  in  the  distributions,  recalling  the  discussion  pertaining  to 

Eq.  (26)  in  Subsection  3.4.3.  Figure  35(a)  and  (b)  are  the  \ *  distribution  for  the  2  Hz  laminar  data,  and  the 

8  Hz,  y^  =28  data  respectively.  Some  information  on  both  the  strength  of  the  filament  and  the  local  strain 
rate  mey  be  accessible  in  this  way. 
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Figure  35 


X~  Distributions  for  Laminar  and  Turbulent  Flow. 


The  elucidation  of  spatial  and  temporal  vorticity  structure  is  considered  in  the  following  sections  describing 
measurements  in  the  inner  boundary  rotation  regimes  of  Couetle  flow. 

5.3  Inner  Cylinder  Rotation  Data  in  Water 

We  return  now  to  the  Rj  axis  of  Figure  26(a).  The  flowfield  with  inner  boundary  rotation  is  of  fundamentally 
different  character  as  described  in  Subsection  5.1.  It  is  inherently  unstable  and  viscosity  inhibits  the  breakdown 
of  flow  until  a  critical  Reynolds  number  Rc,  is  exceeded.  The  subcritical  constant  vorticity  flow  data  has 
already  been  presented  in  Section  4.  Those  measurements,  as  well  as  the  remaining  measurements  described  in 
the  following  sections  were  performed  in  water,  using  the  new  acrylamide  probe  particles. 
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f  or  ihese  measurements  performed  near  the  wail  of  the  outermost  (stationary)  boundary,  we  return  to  th 
.ontigurauon  ot  Figure  9  in  Chapter  3.  Thus  the  relationship  between  the  pnmed  VOP  vorttcity  basis  v 
.ind  the  local  tlow  tield  cylindrical  coordinates  is  u >zt  =  uiz,  u'v.  =  -o)r,  and  cox.  =  wg,  which  is  usually  i 
in  which  most  compulations  are  performed  in  Couette  flow  simulations  (e.g.,  Marcus  (1984)).  Because  i 
angles  of  optical  entry  and  exit  in  this  configuration,  all  of  the  index -matching  techniques  discussed  in  Si 
were  employed  at  the  optical  interfaces,  including  an  index-matched  molded  window  section  (PDFOMA 
an  acrylamide  gel  prism.  The  entire  optical  tram,  including  the  incident  beam-steering  optics,  was  mour 
micrometer  driven  translation  stages,  so  that  any  point  in  the  flow  field  was  reproducibly  selectable,  pan 
up  and  down  the  span  (z-direction).  The  gel  pnsm  simply  slides  up  and  down  the  span  with  the  optics. 

5.3.1  Wavv  Vortex  Flow 

Flow  visualizations  at  three  supercntical  Reynolds  numbers  are  shown  in  Figure  36.  Figure  36(a)  shows 
standard  stationary  Taylor-Couette  rolls  at  Rj/Rc  *  1.1.  These  rolls  have  stationary  inflow  and  outflow 
boundaries  corresponding  to  the  positions  where  the  streamwi.se  vorticity  (cog)  changes  sign.  The  photogi 
Figure  36(b)  is  R/R.  *  15.  This  is  near  the  onset  of  the  wavv-vortax  regime,  where  the  Taylor  couet 
exhibit  periodic  azimuthal  traveling  waves.  Figure  36(c).  at  Rj/Rc  *5.0,  still  in  the  wavy-vortex  regime 
a  disturbance  in  the  flow  due  to  change  w  azimuthal  wavenumber  propagating  up  the  span  as  the  flow  s< 
into  a  stable  state. 

The  value  of  Rc  for  this  system  is  s  127.  For  water  at  25°C.  R^/R,.  =  5  corresponds  to  an  inner  boun 
rotation  rate  of  0.35s  Hz  which  lies  squarely  in  the  wavy  vortex  regime.  This  is  a  convenient  starting  pc 
senes  of  vorticity  measurements  since  the  amplitude  of  the  vorticity  fluctuation  are  strong  relative  to  the 
and  yet  this  is  sufficiently  far  from  the  quasi-penodic  regime  to  minimize  possible  flow  irregularities.  A 
penod  of  time  the  flow  settled  into  one  of  several  preferred  states  (Coles,  1965)  where  it  remained  stabl 
thereafter.  If  the  flow  was  stopped  during  the  course  of  the  experiments,  great  care  was  taken  to  insure 
upon  restart  it  reached  the  same  final  state.  This  proved  to  be  difficult  since  the  azimuthal  wave  states  a 
number  of  span  wise  rolls  (which  are  rather  degenerate  at  large  height/gap)  must  be  reproduced.  Figure 
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Flow  Visualization  Photographs  for  Inner  Boundary  Rotation 


the  (low  visualization  photo  at  R,'RC  =  5.  has  not  yet  reached  a  stable  state,  showing  a  pronounced  irreg 
When  the  tlow  had  settled  in  to  7  azimuthal  traveling  waves  and  15.5  waves  (31  Taylor  vortices)  on  the 
vomcity  measurements  were  performed  at  8  Hz  station  at  a  depth  of  y  =  0. 1  g  (one-tenth  of  the  gap)  fr< 
outer  boundary.  The  wavy-vortex  penod  for  this  state  was  measured  to  be  1.357s,  and  the  span  wise  ■ 
wavelength.  \  =  1.4  cm.  Figure  37  is  the  measured  mean  vorticity  components  versus  z/X.  Again,  for 
hundred  points,  the  uncertainties  in  the  values  of  wz  and  ux(=u)q)  are  a  few  percent,  but  the  uncertainty 

iv(=wr)  is  comparable  to  its  value. 

Since  the  time  at  which  the  individual  data  points  are  recorded  is  also  logged,  sparsely  and  randomly  san 
penodic  dau  can  be  used  to  reproduce  the  temporal  signature  of  the  (vorticity)  wave.  Figure  38  is  singl 
of  -  50  points  acquired  in  just  over  90  seconds.  Thus  the  mean  sampling  rate  is  less  than  the  period.  B 
of  the  accurately  measured,  stable  wave  period,  however,  the  time  behavior  of  the  vorticity  during  a  sing 
cycle  can  be  reconstructed  by  folding  the  record,  modulo-T  into  a  single  cycle.  (  It  is  preferable  to  avoid 
situation  through,  by  sampling  adequately  in  the  first  instance.)  Figure  39  is  a  synthetic'  time  history  o: 
wavy  vortex.  The  trend  is  clear  though  a  bit  noisy.  There  is  a  distinct  phase  relationship  between  the 
components  which  cannot  be  easily  interpreted  with  extensive  fiowfield  mapping. 


Wavy  Vortex  Spatial  Vorticity  Structure  (R/Rc  =  5,7  wave  azimuthally  and  15.5  on  the  span) 
Span  wise  variation  of  mean  vorticity  component  values  near  the  wall, 
over  one  spatial  wavelength  of  wavy  vortex  flow  y/g  =  0.1,  X  =  1.4  cm 
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Figure  38 

Wivy  Vortex  Flow  (water)  Vorticity  Versus  Time  Sample  Record 
•  In  data  collected  over  2  min  interval,  no  periodicity  is  visible 
-  Wavy  vortex  period,  r,  measured  1.347s  (viewing  image  plane) 
-  Mean  sampling  period  —  2s 


This  method  is  reminiscent  of  the  phase  averaging  procedures  used  in  the  eduction  of  coherent  structures  by 
conditional  sampling.  In  this  case,  the  problem  is  a  rather  trivial  one  since  the  data  is  periodic  and  the  phase 
information  can  thus  be  stored  along  with  the  data.  Some  phase  jitter  may  be  present  m  the  data  due  to 
incipient  instabilities.  The  7  azimuthal  wave  flow  at  Rj/Rc  =  5  appears  to  be  only  barely  stable,  owing  perhaps 
to  vibration  or  other  flow  irregularities.  A  Reynolds  number,  the  7  azimuthal  wave  flow,  is  not  expected  to  be 
stable.  Since  the  flow  is  essentially  laminar,  pair  interactions  have  not  been  ruled  out  as  additional  sources  of 
noise.  Finally,  at  the  outflow  boundaries  (near  reduced  time=0)  where  the  vorticity  is  large,  the  velocity  may 

be  larger  still,  resulting  in  reduced  j  and  thus  reduced  precision  in  vorticity  measurement.  The  size  of  the 

error  bars  (nns  distribution  widths,  not  shown)  associated  with  the  points  in  Figure  39  vary  systematically 
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Figure  39 


Wavy  Vortex  Periodogram  (R/Rc  =  5,  7  azimuthal  waves,  t  =  1.347) 


across  the  figure  with  the  minimum  near  the  center.  This  sampling  bias,  which  is  also  periodic,  is  like! 
affect  the  mean  vorticity  profiles,  particularly  for  a i'v.  Correction  strategies  have  been  discussed  in  Ch. 


The  temporal  autocorrelation  is  defined  as 

/  UiWtojft+Tjdt 

^ - 

/"  t** 

Autocorrelations  for  discrete  random  data  have  been  described  by  Frish  (1981)  and  others  and  will  not  f 
recounted  here,  except  to  note  that  the  time  axis  is  simply  divided  uniformly  into  bins  of  sufficient  widti 
obtain  reasonable  statistics  in  each  bin  (properly  normalized  including  statistical  weights)  and  still  insure 
adequate  time  resolution.  A  necessary  refinement  of  discrete  auto-correlations  for  finite  time  events  is  i 
time  weighting  (Buchhave,  et  al.,  1979),  which  is  the  auto-correlation  anaolog  of  transit  time  weighting 
uniformly  spaced  bins  are  then  suitable  for  Fourier  transformation  in  order  to  estimate  the  power  spectn 
the  random  signal.  This  process  has  been  carried  out  with  the  discrete  time  approximation  of  Eq.  (38)  i 
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ensemble  of  records  like  those  in  Figure  39.  Figure  40  is  an  autocorrelation  of  coz  for  8  records  (a  total  of 
-400  points  distributed  over  -  80  bins).  The  reflection  symmetry  is  an  artifact  of  the  wrap-around  algorithm 
employed.  While  the  autocorrelation  must  be  unity  at  r=0,  the  reduced  amplitude  everywhere  else  is  due  to  the 
partial  loss  of  coherence  attributable  to  uncorelated  noise.  A  fitted  cosine  function  shows  that  65  percent  of 
4>z2,  or  ~  80  percent  of  the  rms  Atuz  is  in  the  periodic  structure  of  the  signal.  (Approximately  40  percent  of 
or  64  percent  of  Auig.  is  due  to  the  periodicity.)  Figure  41  is  the  FFT  of  the  autocorrelation  which 
contains  no  additional  information,  (t  should  be  noted,  however,  that  with  sufficiently  long  records  and 
adequate  sampling,  these  power  spectra  could  potentially  show  as  much  temporal  detail  as  the  velocity  records 
in  the  famous  study  of  Fenstermacher  et  al.  (1979).  The  difference  that  vorticity  vector  measurement  could 
make  is  that  amplitudes  and  phases  of  frequency  components  can  be  directly  associated  with  the  dynamical 
variables  that  drive  the  flow  physics  and  thus  the  instabilities. 

5.3.2  Turbulent  Tavlor  Vortex  Flow 

By  increasing  R1/Rc  an  order  of  magnitude  to  50,  all  vestiges  of  periodicity  in  the  flow  are  lost,  and 
one  is  well  into  the  turbulent  Taylor  vortex  flow  regime.  A  striking  feature  of  this  flow,  however,  is  the 
spanwise  wave  structure,  which  is  clearly  visible  in  the  flow  visualization  photo  in  Figure  42.  At  much  higher 
values,  this  'coherent  structure'  still  persists,  though  less  visibly,  in  the  presence  of  intense  turbulence. 


<i»z  Autocorrelation  for  Wavy  Vortex  Flow  (R/Rp  =  5,  7  azimuthal  waves,  r  =  1.34/s) 
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Linear  Power  Spectrum  of  Span  wise  Vorticity 
(from  FFT  of  autocorrelation  in  arbitrary  units) 


Flow  Visualization  Photograph  of  Turbulent  Taylor  Vortices,  R^/Rc  =  50 
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The  R,  R.  =  50  spanwise  structure  is  again  measured  across  one  span  wise  wavelength  at  y  =  O.lg.  The 
histograms  of  these  measurements,  as  well  as  the  wz,tog  scatterplots,  are  included  as  a  function  of  position  on 
the  span  in  Figure  43. 

The  scatterplots  show  an  alternating  pattern  of  left  and  right  inclined  vortex  lines  near  the  plane  of  the  wall, 
which  is  to  be  expected  from  counter-rotating  streamwise  vortices  in  the  mean  flow.  The  fluctuations,  however, 
seen  clearly  in  the  scatterplots.  exhibit  several  interesting  features  .  The  spanwise  vorticity  fluctuations  are  most 
intense  near  inflow  and  outflow  boundaries,  (where  uq  =0).  Figure  44(a)  and  (b)  are  the  mean  and  rms 
vorticity  values  for  these  histograms.  The  near  vanishing  of  the  spanwise  vorticity  at  z/X  —  0.7  is  apparently 
the  inflow  boundary  position,  where  low  speed  fluid  accumulates  at  the  outer  wall.  The  small  values  of  ojx  and 
-l'z  result  m  enormous  uncertainties  for  uv  (noted  in  the  Figure  44(b)  as  the  ‘two-component  range  ).  The  rms 
fluctuations  normalized  by  the  magnitude  of  the  mean  vorticity  |u|  are  fairly  uniform  along  z.  suggestive  of 
relatively  homogeneous  through  highly  anisotropic  turbulence.  Again,  detailed  flowfield  mapping  is  needed  for 
more  definitive  descriptions. 

At  still  higher  Reynolds  numbers,  the  fluctuations  increase  dramatically  in  intensity.  A  few  histogams  are 
vhown  in  Figure  45  at  R^/Rc  =  150  for  the  condition  indicated  in  the  figure.  The  vorticity  fluctuations  in  all  of 
the  turbulent  Taylor-vortex  flow  data  show  strong  interactions  with  the  local  mean  flow.  The  eduction  of  the 
prevalent  structures  as  the  turbulence  accomodates  itself  to  the  local  flow  would  be  an  interesting  study. 

5.4  Summary  of  Couette  Flow  Measurements 

Three-component  vorticity  measurements  have  been  performed  in  the  both  outer  and  inner  boundary  rotation 
type  Couette  flow.  Summaries  of  observations  about  the  structure  of  the  vorticity  fields  in  each  of  these  flows 
are  given  below. 
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Turbulent  Taylor  Vortices  R^/Rc  =  50:  Spanwise  Profile  (y/g 
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Figure  43  (Continued) 

Turbulent  Taylor  Vortices  =  50:  Spanwise  Profile  (y/g  =  0.1) 


Figure  44 

Spatial  Variation  of  Normalized  RMS  Vorticity  Fluctuations  Near  the  Wall 
Corrected  for  Instrumental  Noise 
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Figure  45 

Turbulent  Taylor  Vortex  Flow  R/R,,  =  150 
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5.4.1  Outer  Boundary  Rotation 

For  outer  boundary  rotation  Couette  flow,  laminar,  intermittent,  and  turbulent  vorticity  distributions  were 
measured.  It  was  found  the  laminar  Couette  flow  exhibited  uniform  vorticity  throughout  the  flowfield,  at  the 
values  expected  from  the  exact  solution  of  the  equation  of  motion.  Intermittent  turbulence  appeared  and  was 
detected  at  the  anticipated  Reynolds  number.  Fully  developed  turbulence  has  mean  and  rms  fluctuating  vorticity 
profiles  in  y,  the  coordinate  normal  to  the  inner  stationary  bounday,  which  resemble  standard  plane-boundary 
layer  turbulence  very  close  to  the  wall.  Farther  from  the  wall,  the  vorticity  fluctuations  appear  to  be  dominated 
by  the  fluctuating  component  of  vorticity  which  lies  parallel  to  the  principal  axis  of  the  mean  flow  rate-of-strain 
tensor.  Joint  probability  density  functions  of  the  turbulent  vorticity  show  evidence  of  structures.  The  z  and  s 
component  of  the  vorticity  are  only  weakly  correlated,  though  apprently  not  independent  (Tennekes  & 

Lumlev,  1972).  It  is  difficult  to  distinguish  structures  from  background  vorticity,  however,  and  these 
observations  await  improved  statistics  from  greater  quantities  of  data  for  clarification. 

5.4.2  Inner  Boundary  Rotation 

The  inner  boundary  rotation  experiments  served  as  both  a  test  of  the  new  acrylamide  vorticity  probe  particles  in 
water,  and  a  test  of  the  VOP's  ability  to  distinguish  temporal  and  spatial  vorticity  field  structure  in  the  presence 
of  noise  or  turbulent  flow  fluctuations.  The  laminar  inner  boundary  rotation  Couette  flow  demonstated  that  the 
acrylamide  particles  behave  as  expected,  and  also  were  a  clear  indication  of  the  particle  pairing  behavior 
anticipated  in  laminar  shear  flow.  This  pairing  is  observed  to  be  significantly  disrupted  by  the  presence  of 
turbulent  flow  fluctuations.  The  wavy-vortex  flow  experiments  demonstrate  the  elucidation  of  spatially  and 
temporally  periodic  vorticity.  These  measurements  are  limited  by  the  currently  available  maximum  record  lenth 
that  can  be  acquired.  Unfortunately,  this  limitation  neccessitated  skipping  over  most  of  the  interesting  range  of 
Reynolds  number  where  the  spectral  evolution  to  turbulence  transpires,  since  several  incommensurate 
frequencies  are  difficult  to  resolve  without  long  temporally  continous  records.  Lastly,  turbulent  Taylor  vortex 
flow  exhibits  fascinating  interactions  of  the  vorticity  fluctuations  with  the  mean  flow,  and  deserves  further 
study. 
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6  CONCLUSION 


o.  1  Summary 

A  system  tor  simultaneous  measurement  of  the  three  components  of  the  local,  instantaneous  vorticity  v'ect 
laminar  and  turbulent  liquid  flow  has  been  assembled  and  tested.  This  system  is  an  extension  of  the  sing! 
component  Vorticity  Optical  Probe  system  developed  by  Frish  and  Webb.  The  most  powerful  features  ot 
approach  to  vorticity  measurement  are  its  superior  spatial  resolution,  limited  only  by  the  size  of  the  probe 
particles,  and  its  outstanding  performance  in  high  vorticity  and  low  velocity  regions  without  the  need  to  t: 
Taylor's  hypothesis.  In  this  respect,  the  VOP  is  very  complementary  to  other  velocity-based  vorticity 
measurement  methods.  Not  only  are  the  three  vorticity  components  measured,  any  two  of  them  more  act 
ihan  the  third,  but  uncertainties  in  the  vorticity  components  as  well.  Other  useful  statistical  parameters  ai 
derived  during  data  reduction.  This  auxiliary  information  allows  vorticity  resolution  and  sampling  functu 
he  derived  as  function  of  experimental  parameters.  The  three-component  VOP  is  'self-characterizing',  in 
sense  that  it  measures  its  own  precision  as  it  acquires  data. 

The  VOP  is  now  water-compatible.  The  invention  of  a  new  class  of  probe  particles  made  of  acrylamide 
offers  the  possibility  of  new,  widespread  applicability  of  the  VOP,  particularly  in  larger  scale,  high  Reyn 
(low  facilities.  The  new  acrylamide  particles  are  ideally  suited  to  large  scale  application,  since  they  are 
-  90  percent  water  by  weight  and  readily  mass- producible  at  low  cost. 

The  three-component  system  has  been  tested  in  six  distinct  Couette  flow  regimes.  Laminar,  transitional  ar 
turbulent  flows  were  studied  for  inner  and  outer  boundary  rotation.  The  laminar  flows  exhibit  the  expecu 
spatially  and  temporally  constant  vorticity  anticipated  from  the  analytic  Couette  flow  solution  of  the  equal 
motion.  For  inner  boundary  rotation,  the  spatial  and  temporal  behavior  of  the  vorticity  components  the  n 
outer  boundary  in  wavy-vortex  flow  at  R^/Rc  =  5.0  has  been  measured.  Comparisons  with  recent  numer 
simulations  of  this  flow  compare  favorably  (Coughlin,  1991).  Turbulent  taylor  vortex  flow  shows  strong 
interactions  between  the  fluctuating  vorticity  and  the  local  mean  flow.  Outer  boundary  rotation  data  at 
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R0  =32,000  (based  on  gap  and  surface  speed)  shows  very  large  spanwise  vorticity  fluctuations  at  the  wail, 

>  40  percent,  after  corrections  for  instrumental  broadening.  These  measurements  also  compare  well,  close  to 
the  inner  wall,  with  direct  numerical  simulations  of  turbulent  boundary  layers.  Farther  from  the  wall  the 
difference  in  these  flow  types  becomes  evident.  The  mean  background  shear  in  the  Couette  flow  dramatically 
amplifies  the  vorticity  fluctuations  in  the  s-direction,  inclined  at  45  deg  to  the  flow  direction  which  is  also  the 
principal  strain  direction.  This  observation  agrees  with  recent  simulations  of  turbulent  shear  flow,  which 
suggest  that  hairpin  vortices  are  an  important  feature  of  all  turbulent  shear  flows.  Scatterplot  representations  of 
joint  aiz,u>s  PDFs  show  some  evidence  of  intense  vortex  filaments  rotating  toward  the  principal  strain  axis, 
though  this  must  be  confirmed  with  further  measurements.  These  measurements  show  the  potential  power  and 
versatility  of  the  VOP  technique,  and  invite  new  approaches  to  old  problems. 

6.2  Limitations  of  the  Three-component  VOP 

As  with  any  fluid  flow  measurement  system,  there  are  a  number  of  fundamental  limitations  which  the  single 
point  VOP,  by  the  nature  of  its  operation,  can  never  overcome. 

The  sense  in  which  vorticity  measurements  are  time-resolved  is  somewhat  different  from  other  types  of 
measurement.  The  VOP  trades  off  some  temporal  resolution  to  ’see'  the  rotation  of  the  particle-containing  fluid 
element,  in  order  to  avoid  the  loss  of  spatial  resolution  which  is  invariably  the  price  of  finite-difference  methods 
of  velocity  gradient  measurement.  The  vorticity  is  measured  at  a  probe  particle’s  location  which  is  convecting 
through  the  sampled  volume.  It  may  move  a  considerable  distance  during  the  measurement,  depending  on  its 
vorticity  and  the  allowed  sampled  volume  dimensions.  This  is  the  sense  in  which  the  VOP  is  locally 
Lagrangian',  and  the  Lagrangian  microscale  of  time  is  the  limiting  period  for  accurate  vorticity  measurement. 

Up  to  this  point,  this  does  not  really  qualify  as  a  limitation.  Longer  time  of  observation  is  required  for  vorticity 
of  lesser  magnitude,  but  the  convection  velocity  through  the  sampled  volume  does  not  necessarily  cooperate. 

The  vorticity/velocity  ratio  problem  and  the  VOP  'singularity’  (w |  [^)  are  closely  related.  Short  trajectories  at 

finite  signal-to-noise  are  almost  invariably  interpreted  as  having  large  (spurious)  curvature,  i.e,  orbiting  about  a 
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nearby  point,  reflecting  the  indeterminate  nature  ot  the  a>v'  component  at  a  point.  For  larger  trajectorie.c 
-ituauon  improves  rapidly  giving  a  useful  measure  of  this  component,  albeit  with  greater  uncertainty  than 
■ihers.  Fhe  vorticity/ velocity  ratio  often  drops  below  the  lower  limit  of  usefulness  (typically  from  0. 1  to 
1  ^m  ‘  depending  on  the  sampled  volume  dimensions).  Here  one  must  either  make  assumptions  about  th 
vorticity  field  in  order  to  measure  low  vorticitv  values  (e.g.,  that  it  is  essentially  two-dimensional,  which 
usually  a  good  approximation  when  fitting  only  two  components  near  detector  center),  or  put  a  count  in  a 
special  bin  for  near  zero  vorticity  and  move  on.  At  the  end  of  the  data  reduction,  if  there  is  a  sufficient 
of  counts  in  this  but.  then  it  may  be  statistically  and  fluid  dynamically  significant,  e.g.,  an  irrotational  (Z 
spike  in  the  PDF  due  to  mtermittency. 

The  VOP  has  sampling  limitations.  Overlap  of  finite  duration  trajectories  trajectory  are  not  allowed  in  thi 
currently  used  detection  scheme.  The  mean  time  between  trajectories  must  be  at  least  three  times  the  me 
detector  transit  time.  This  implies  an  average  spacing  of  the  probe  particles  in  the  flow,  and  thus  typical 
the  structures  that  will  be  resolved.  This  statement  makes  the  coi  lection  between  the  local  Langrangian 
vorticity  measurements  and  the  Eulerian  structure  of  the  flow.  But  here  is  Taylor's  hypothesis  is  creepin 
in:  the  notion  the  spatial  structure  can  be  derived  from  a  time  sequence  of  values  of  a  given  quantity  at  a 
point  in  the  flowfield  is  still  a  subject  of  debate,  particularly  in  turbulent  shear  flows.  Though  the  VOP 
eliminates  this  controversy  for  individual  vorticity  measurements,  these  questions  pertaining  to  extended  ! 
'tructures  persist,  as  for  all  other  single  point  measurements  of  flow  variables. 

6.3  New  Measurements,  New  Methods 

As  an  expansion  of  single  point  measurement  capability,  a  second  detector  could  be  added,  as  indicated  i. 
Figure  9  in  Chapter  3,  in  order  to  perform  two-component  velocity  measurements  simultaneously  with  th 
component  vorticity  measurements  for  the  same  particle  in  the  sampled  volume.  Instantaneous  vorticity-' 
correlations  would  thus  be  accessible,  including  significant  u  and  ut u  contributions  and  components.  T 
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-ccond  detector  cun  tv  iu>t  <ls  easily  configured  as  a  second  VOP.  Two-point  spatial  vorticity  correlation  is 
-till  a  measurement  ot  tundamentat  importance  in  the  study  of  coherent  structures  in  turbulence. 

Finalise  it  appears  that  the  V  OP  must  one  day  go  the  way  of  many  other  measurements— toward  multiple  point 
measurements,  conditional  sampling,  \ phase  average;  methods,  and  ultimately  imaging  techniques.  The 
photograph  in  Figure  4o  below  is  of  VOP  probe  particles  in  a  shear  flow  illuminated  bv  a  sheet  of  laser  light. 
Sequences  ot  bnght  flashes  along  the  apparent  streamlines  indicate  the  presence  of  vorticity.  As  the  particles 
route,  reflections  from  some  enter  the  camera  lens  producing  a  bnght  spot.  Only  a  time  base  and  a  length 
scale  are  needed  to  extract  \ortictties.  The  logical  extrapolation  from  these  concepts  is  direct  vorticity 
.  isuaiization  techniques. 

0.4  Closing 

Another  tool,  the  three-component  water-compatible  Vorticity  Optical  Probe,  has  now  been  added  to  the 
impressive  collection  already  marshalled  in  the  pursuit  of  a  fundamenul  understanding  of  fluid  flow.  The  new 
VOP  is  uniquely  suited  to  the  problems  of  high  shear,  high  Reynolds  number,  non-parallel  and  near-wall  flows 
that  have  plagued  the  expenmenul  fluid  dynamics  community  for  decades.  It  is  to  be  hoped  that  this  new 
building  block  will  enlarge  our  knowledge  of  this  most  common  and  most  complex  phenomenon  —  turbulence. 


Figure  46.  Vorticity  Visualization? 
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Appendix  A 


VOP  Data  Reduction  and  Hardware  Simulation  Program 


C 

VERS.  6  DEC  90 

C 

C 

SLARGE 

C  ARRAY  AND  COMMOM  DECLARATION 

c 

REAL *4  RDATAX(1024) 

REAL *4  DDATAX(1024) 

REAL*4  SDATAX(1024) 

REAL *4  RDATAY(1024) 

REAL *4  DDATAY(1024) 

REALM  SDATAY(1024) 

REALM  CHIOPD(IO) 

INTEGERS  SQRFLG 
COMMON/FLG/SQRFLG 

COM  MON/WD  AT /W(  1024),  Y(  1024),  SINTX!  1024), SINTY(  1024),  NCOUNT 

common/cdat/period.scale.ydet.snoise.btx,bty,cdhi,cdlo,cutmax,o 

COMMON/OMEG/TM(  1024),  WX(  1024),  WY(  1024).  WZ(  1024).NU(  1024) 

COM  MON/ST  AT/D  WX(  1024),  DWY(  1 024) ,  DWZ(  1 024),CH12(  1 024) 
COMMON/MAP1/VXMAP(20.20).VZMAP(20.20) 

COMMON/M  AP2/XP(2Q) ,  ZP(20) 

COM  MON/TIM  R/TOTAL,  LST.TSC  A.  NCSEG 
COMMON/BLIN/BERMAX.ISCRAP.1BENAB 
COMMON/PROX/PROXC.PROXW.PROXP.CLIPR 
CHARACTER*20  FXS.FXD.FYS.FYD.TSMP.ZNAME.FMAP 
INTEGER*2  OFF6,OFF7,OFF8.OFF9,OFF10 
LOGICAL  LPOS . LTRS , LROUND , LTRS 1 . LTEST 
DATA  1DFILE/5/ 

DATA  NSM/7/ 

DATA  NWORDS/1024/ 

FMAP=  C:\FORT50\MAP.DAT' 

C 

C 

26  WRTTE(*,’(A\)')’  ENTER  X  SUM  FILENAME;  ' 

READ  (V(A)\ERR  =  26)  FXS 


WRJTE(*,’(A\)')’  ENTER  X  DIFF  FILENAME; 
READ  (V(A)’.ERR  =  27)  FXD 


l  WR1TE(*, '(A\)')'  ENTER  Y  SUM  FILENAME;  ' 
READ  (*,'(A)’,ERR  =  28)  FYS 


29  WRITE(*,’(A\)y  ENTER  Y  DIFF  FILENAME;  ’ 
READ  (*,’(A)’,ERR  =  29)  FYD 
C 
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WRITE)*. ’(A  )•)•  ENTER  TIMESTAMP  FILENAME;  ’ 

READ  (*.’(A)’.ERR  =  29i  TSMP 

C 

i:  READ  DETECTOR  MAP  DATA  FILE 

C 

OPENUDFILE.  FILE  =  FMAP.  STATUS  =  OLD’,  FORM  =  FORMATTED  ) 

DO  13  L  =  1.20 

READ(IDFILE.*’PZ 

ZP(L)  =  PZ-.863 

XP(L)  =  . 909-0. 1  *(L-  h 

DO  12  M=l.20 

READ(IDFILE.*)V\MAP(M.L).VZMAP(M.D 

12  CONTINUE 

13  CONTINUE 
CLOSE(IDFILE) 

OPEN  DATA  FILES 

OPENI6.  FILE  =  FXS.  STATUS  =  OLD’.  FORM  = ’UNFORMATTED’) 
CALL  CAT(6) 

OPEN(7,  FILE  =  FXD.  STATUS  =  OLD’.  FORM  = ’UNFORMATTED’) 
CALL  CAT(7) 

OPEN18.  FILE  =  FVS.  STATUS  =  OLD’,  FORM  = ’UNFORMATTED’) 
CALL  CAT(8) 

OPEN(9.  FILE  =  FYD,  STATUS  =  OLD’.  FORM  = ’UNFORMATTED’) 
CALL  CAT(9) 

OPENdO,  FILE  =  TSMP.  STATUS  =  OLD’.  ACCESS=  SEQUENTIAL’, 

I  FORM  =  BINARY’) 

CALL  CAT(IO) 


SET  INITIAL  CODES  AND  COUNTERS 


NEV  =  I 
ISEG  =  0 
IBDSG1  =999 
IBDSG2  =  999 
IAUTO=0 
ISM  =0 
ITURN  =0 
NEND=0 
NOPT  =  2 
IFLAG=0 
IFLAG2  =  0 
ISCRAP=0 
IBENAB=0 


WRITE(*,’(A\)’)’  ANY  BAD  SEGMENTS? 
READ(*,’(I4)’)IBDSG 
IF(1BDSG.EQ.  I)THEN 
WRITE(*,’(A\)’)’  BAD  SEGMENT  1 
READ(*,  *(I4)’)IBDSG  1 


123 


nnnnon  non  o  n 


WR1TE(V(A  >')•  BAD  SEGMENT  2 

READ(*.-(I4)-)IBDSG2 

ENDIF 

C 

WR1TE(*.’(A\)T  CHANGE  DEFAULT  SUM  NOISE?(l  =YES);- 
READ)*, '( 14 ) '  }ICH 
IF(ICH.EQ.  DTHEN 

3?  WRITEfVfADT  ENTER  SUM  NOISE  LEVEL  (VOLTS);- 
R£AD(*,.-(F12.7)',ERR  =  37)SNOISE 
IF(SNOISE.LT.0.001)GO  TO  37 
ENDIF 


YDET  =  -I  9 

WRITE(V(A\)-)-  DETECTOR  PLANE  POSITION?  (DFLT  =  -1 ,9);- 
READ(V(U)*)NDET 

IFfNDET.EQ.  1)READ(*.'(FI2.7)’,ERR  =  38)YDET 
WRrTE(*,'(A\)’)’  Y  =  ’ 

WRITE(*,  1200IYDET 

SET  DEFAULT  PARAMETERS 


NMIN  =  25 
NMAX  =650 
TRS  =  0.S 
TRS1  =TRS 
TOTAL =0.0 
HYS  =  10. 

NSM  =  17 
CDHI  =  20. 

CDLO  =  2.0 
CUTMAX  =5.0 
BERM  AX=  15.0 
CLIPR=0.06 
RCUT  =  2. 

PROXC  =  .990 
PROXW=  998 
PROXP=  994 
QDECR  =0.065 
Q=0.65 
QDEF=Q 

START  BIG  LOOP 

DO  11  1CSEG=  l.NCSEG 
ISEG  =  ISEG  + 1 

BASELINE  SUBTRACTION 

CALL  BASE(6,OFF6,SDATAX,NEND,TTME) 
IF(NEND.EQ.  DIFLAG  =1 
CALL  BASE(7,OFF7,DDATAX,NEND,TIME) 
IF(NEND.EQ.  DIFLAG  =1 
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CALL  BASEfS.OFFS.SDATAY. .VEND. TIME) 

IF(NEND.EQ.  DIFLAG  =  1 

CALL  BASE(9.0FF9.DDATAY. SEND. TIME) 

!F!NEND.EQ.  I IIFLAG  =  I 

CALL  BASEdO.OFFIO.DDATAY.NEND. TIME) 

TOTAL  =  TOTAL  *  TIME 
IFIIFLAG.EQ.  DTHEN 
NSEG  =  ISEG-I 

WRTTE(*.'(A')')'  END  OF  FILE  ENCOUNTERED  AT  SEG#' 

WRITE!*.  120DNSEG 
GO  TO  107 
ENDIF 

C 

IFHSCRAP.EQ.  UTHEN 

WRITE!*. '(A\)')'  BASELINE  PROBLEM  IN  SEG  #' 

WRITE!*.i:01)!CSEG 
ISCRAP=0 
GO  TO  1 1 
ENDIF 

IF!!ICSEG.EQ.IBDSG1>.0R.!ICSEG.EQ.IBDSG2))THEN 

WRITE!*, '(AN)*)'  SKIPPING  BAD  SEG  #’ 

WRITE!*.  120DICSEG 

GOTO  It 

ENDIF 

C 

1200  FORMAT(Fl2.7) 

1201  FORMATII3) 

C 

IFIIAUTO.NE.  DTHEN 
DO  1  I  =  I.NWORDS 
wn)  =  (I-D*PERIOD 
1  CONTINUE 

ENDIF 

8  IRET  =  0 

IFdAUTO.NE.  DTHEN 

C 

SWITCH  TO  AUTOMATIC  (NON-INTERACTIVE)  REDUCTION  MODE? 

47  WRITE!  V(A))T  INTERACTIVE  OR  BATCH?  (1  =BATCH);’ 

READ(*.,(I1)'.ERR  =  47)IAUT0 

C 

C  VIEW  RAW  DATA? 

C 

251  WRITE!*. ‘(A))’)'  DISPLAY  X  SUM?  (1  =YES);‘ 

READ(*,’(II)'.ERR  =  25DITD 
IF(ITD.E01)CALL  FTDISPfSDATAX.I.O) 

C 

252  WRITE! V(A\)’)'  DISPLAY  X  DIFFERENCE?  (I  =  YES);' 

READ!  V(I1)\ ERR  =  252)ITD 
IF(ITD.EQ.i)CALL  FTDISP(DDATAX.l.O) 

C 
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253  WRITE!*. (AW)’  DISPLAY  Y  SUM.'  !l=YES);’ 

READ(*.’(1 1)‘.ERR  =  253)ITD 

IFdTD.EQ.  DCALL  FTDISP(SDATAY.  1 .0) 

c 

254  WRITE(*.'(AW)’  DISPLAY  DIFFERENCE?  !1  =YES);’ 
READ(*.'dl)’.ERR  =  254)ITD 

I FIITD . EQ.  1  )C ALL  FTDISP(DDATAY,  1 ,0) 

C 

C  CHANGE  DEFAULT  PARAMETER? 

C 

56  WRITE!*.  1 300VTRS 

256  WRITE!*. ’(AW)’  SET  NEW  THRESHOLD?  (1  =  YES);' 

RE  ADI  •.’(11)’,ERR  =  56)ICH 

1F(ICH.EQ.1)THEN 
R£AD(*.’(FI2.5)’,ERR  =  56)TRS 
IF(TRS.LT.O.OI)GO  TO  56 
ENDIF 

57  WRITE!*.  130DHYS.NSM 

257  WRITE!*. '(A\)')’  SET  NEW  VALUES?  (1=YES);’ 

READ!*,‘(I1)',ERR  =  57)ICH 

IFdCH.EQ.  1)THEN 

WRITE(*.'(A\)’)'  SET  HYSTERESIS! +  /-  %  OF  THRESHOLD);  ’ 
READ!  •  '(F 1 2.5)’  .ERR  =  57)HYS 
WRrrEC.'(A\)')'  SET  SMOOTHING  PAD  LENGTH;  ' 
READ!*.‘(I3)’,ERR  =  57)NSM 
ENDIF 
C 

IF(IBENAB.EQ.O)THEN 

WRITE!*. '(AW)’  BASELINE  CHECKER  IS  OFF.  TURN  ON?(l=YES)’ 
READ(*.'(II)')IBENAB 
ENDIF 
C 

58  WRITE!*.  1302)BERM  AX 

258  WRITE!*. '(A\)’)'  RESET  MAX  BASELINE  ERROR?  (1=YES);' 
READ!*.’(II)’,ERR  =  58)ICH 

IFdCH.EQ.  1)THEN 
READ(*,'(F12.5)'.ERR  =  58)BERMAX 
ENDIF 
C 

WRITE!*.2I2)QDECR 

WRITE!*, ’(AW)’  NEW  EDGE  DECREMENT?  (1=YES);’ 
READ(*,’(I1)’)ICH 

IFdCH.EQ.  1)READ(*.’(F12.5)’)QDECR 
C 

WRITE(*.212I)CLIPR 

WRITE!*, ’(A\)’)’  NEW  END  CLIPPING  FRACTION  (1  =  YES);* 
READ(*,’(I1)’)ICH 

IFdCH.EQ.  1)READ(*,’(F12.5)’)CLIPR 
C 

WRITE(*.21I)RCUT 

WRITE(*.’(A\)’)’  NEW  FIT  IMPROVEMENT  RATIO  ?  (1  =  YES);’ 
READ(*,’(H)’)ICH 

IFdCH.EQ.  i)READ(*,’(F12.5)’)RCUT 
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WRITE!*. ’(A  )V  CHANGE  PROXIMTY  COSINESS  1  =  YES);1 

READ(*.'(I4)')ICH 

IFdCH.EQ.  DTHEN 

WRITE!*, ‘(ADT  ENTER  Wlln  COSINE  (DEF  .998);1 
READ(V(F12.7)’)PROXW 

WRITE!*, ‘(AO)1  ENTER  n  il  C  COSINE  (DEF  .990);1 
READ(*.'(FI2.7)')PROXC 

WRITE!*. '(ADT  ENTER  W  II  C  COSINE  (DEF  .994);1 

READ(*.’(FI2.7)HPROXP 

ENDIF 

C 

c 

WRITE  (*,  1303)NMIN 

WHITE!*. '(AD1)1  NEW  MIN  TRAJ  LENGTH  ?  (l=YES);’ 

READ(*.'(I1)')ICH 

IFdCH.EO- 1  )READ(  *.  '(I3)')NMIN 

WRITE  (*.  1 304)NMAX 

WRITE!*,1!.*))1)1  NEW  MAX  LENGTH  ?  ( I  =  YES);’ 

READ(*.'(I1)’)ICH 

IFdCH.EO-  l)READ(*.'d3)')NMAX 

C 

ENDIF 

1 300  FORMATC  DEFAULT  THRESHOLD  IS  CURRENTLY1. F8.4.) 

1301  FORMATC  HYSTERESIS  CURRENTLY1, F8.4.1,  WITH1. 13. 1  PT  SMOOTHING1) 

1 302  FORMATC  DEFAULT  BASELINE  ERROR  (NOISE  LEVELS)1, F8. 4.) 

1 303  FORMATC  DEFAULT  MIN  TRAJ  LENGTH  IS  CURRENTLY  M3.) 

1304  FORMATC  DEFAULT  MAX  TRAJ  LENGTH  IS  CURRENTLY  M3.) 

C 

IFdAUTO.NE.  DTHEN 

C 

DO  10  J  =  I.NWORDS 

LTRS  =  (SDATAX(J)  -t-  SDATA  Y(J)).GT.TRS 

LPOS  =  !SDATAX(J).GT.0.0).AND.!SDATAY(J).GT.0.0) 

!F(LPOS.AND.LTRS)THEN 

RDATAX(J)  =  DDATAXfJVSD  ATAX(J) 

RDATAY(J)  =  DDATAY(J)/SDATA  Y(J) 

ELSE 

RDATAX(F)  =0.0 
RDATAY(J)=0.0 
ENDIF 

10  CONTINUE 

DISPLAY  RAW  QUOTIENTS? 

WRITE!*. '(A))1)1  DISPLAY  X  QUOTIENT?  (1=YES);’ 

READ(V(I1)’)ITD 

IFdTD.EQ.  1)CALL  FTDISP(RDATAX.l.O) 


WRITEC.TAV)1)1  DISPLAY  Y  QUOTIENT?  (I=YES);’ 
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READ(*.'(I1)')ITD 

IFflTD.EQ.  1)CALL  FTDISP(RDATAY.l.O) 

C 

ENDIF 

C 

C 

C  COMMENCE  INNER  LOOP:  THRESHOLDING/M  ASKING/MAPPING/ ANALYSIS 

C 

c 

WRITE!*. '(A\)')’  STARTING  NEW  SEGMENT  #’ 

WRITEf*.I20I)ICSEG 
F  =0.0 

TRS1  =TRS*(1.0-HYS/100.) 

NCT  =  1 

NSTRT  =  (NSM  +  l)/2 
LLST  =  NSTRT+  1 

179  !F(LLST.EQ.(NSTRT+1))THEN 

SUMW2  =0.0 
DO  4  KK  =  l.NSM 

SUMW2  =  SUMW2  +  SDATAX(KK)  +  SD  ATA  Y(KK) 

4  CONTINUE 

ENDIF 

C 

BTX  =0.0 
BTY  =0.0 
SBTX  =  0.0 
SBTX2=0.0 
SBTY  =0.0 
SBTY2  =0.0 
C 

DO  94  L  =  LLST.NWORDS-NSM 
IF(F.EQ.l.O)  GO  TO  96 
LST=L 

NINC  =  L  +  NSTRT-1 
NDEC  =  L-NSTRT 

SUM  W2  =  SUMW2  +  SD  ATAX(NINC)  +  SDATA  Y(NINC)-SD  ATAX(NDEC)-SDATA  YfNDEC 
LTRS  =  (SUMW2/NSM).GT.TRS 
LTRS1  =(SUMW2/NSM).GT.TRS1 
LTEST  =  LTRS 

IFfNCT.GT.  NMIN)LTEST  =  LTRS I 

LPOS  =  (SDATAX(L).GT.0.0).AND.(SDATAY(L).GT.0.0) 

I  F(LPOS.  AND .  LTEST)THEN 
RDATAXfL)  =  DDATAX(L)/SDATAX(L) 

RDATAY(L)  =  DD  ATA  Y(L)/SDATA  Y(L) 

SBTX  =  SBTX  +  SD  ATAX(L) 

SBTX2  =  SBTX2  +  SD  AT  AX(L)*SD  ATAX(L) 

SBTY  =  SBTY  +  SD  AT  A  Y(L) 

SBTY2  =  SBTY2  +  S  DATA  Y(L)*SD  ATA  Y(L) 

WT  =  RDATAXfL) 

YT=RDATAY(L) 

SQRFLG  =  0 

CALL  INTERPfWT.YT) 

IF(SQRFLG.EQ.  l)GO  TO  II 
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W(NCn  =  WT 
Y!NCT)  =  YT 

SINTX!NCT)  =  SDATAX(L) 

SINTY(NCT)  =  SDATAY(L) 

NCT  =  NCTr  l 

LROUND  =  <SQRT!WT*WT  -r  YT*YT).GT.Q) 
IF!LROUND)THEN 
IF(NCT.LE.NMIN)THEN 
NCT  =  1 

RDATAX(L)  =  0.0 
RDATAY(L)  =  0.0 
ELSE 

NCOUNT  =  NCT-2 
F=  1.0 
ENDIF 
ENDIF 

ELSE  IF  (NCT.LE.NMIN)THEN 
NCT  =  t 

RDATAX!L>  =  0.0 
RDATAYfL)  =0.0 
ELSE 

NCOUNT  =  NCT- 1 

F=  1 .0 

ENDIF 

94  CONTINUE 

NCOUNT  =  NCT- 1 
96  F=0.0 

IF(ICHI.EQ.O)LLLST  =  LST  +• 1 
IF!NCOUNT.GE.  DTHEN 
BTX  =  SBTX2/SBTX 
BTY  =  SBTY2/SBTY 
ENDIF 
NCT  =  1 
C 

IFINCOUNT.LE.NMINVTHEN 
WRITE!*.’(A)’)’  TOO  FEW  POINTS!!!’ 

GO  TO  103 

C 

ELSE  IF(NCOUNT.GE.NMAX)THEN 
WRITE!*, ’(A)’)'  TOO  MANY  POINTS!!!’ 
GO  TO  103 

C 

ELSE 

WRITE!*, ’(I4)’)NCOUNT 
C 

IFdAUTO.NE.  DTHEN 


WRITE!*, ’(A\)’)’  DISPLAY  TRAJECTORY?  (1=YES);’ 
READ!*,’(I1)’)ITD 

IF(ITD.EQ.1)CALL  FTDISP(Y,2,NCOUNT) 
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WRITE(*.'( A')')'  DISPLAY  CORRECTED  X  QUOTIENT*’  (1=YES);- 
READlVdlHlTD 

IFHTD.EQ.  1  )CALL  FTDISPfWJ.NCOUNT) 

C 

C 

WRITE!*. *(A\)')'  DISPLAY  CORRECTED  Y  QUOTIENT  ( I  =  YES);' 
READ(*.‘(I1  )')ITD 

IFdTD.EQ.  1)CALL  FTDISP(Y,3.NCOUNT) 

C 

c 

WRP"'.(*.2I9)QDEF 

WR1  .£(*,*(A\n*  NEW  EDGE  CUTOFF  ?  (1  =YES);* 
READ(V(I1)’)ICH 
IFHCH.EQ.  1)THEN 
READ(*,'(F12.5)‘)QDEF 
Q-QDEF 
ENDIF 
C 
C 

WRITE!V(A\)T  ANOTHER  PASS  THROUGH  THIS  SET  (I  =  YES);' 

READ(*.’(I1)')IR£T 

IFflRET.EQ.  l)GO  TO  8 

C 

ENDIF 


IF(lCHl.EQ.l)  THEN 
IFINCLAST.EQ.NCOUNTKjO  TO  101 
ENDIF 
C 

C  CALL  THREE-COMPONENT  VORTICITY  LEAST-SQUARES  FITTING  ROUTINE 

C 

C 

98  CALL  OMEGA(NEV.IAUTO.ITURN.ICHI.CHIOP) 

C 

C 

C  TEST  FOR  OPD  EDGE  EFFECTS 

C 

101  IFdCHI.EQ.  1)THEN 
NCHI  =  NCHI+ 1 
CHIOPD(NCHI)  =  CHIOP 
IFINCHI.EQ.  I)THEN 

WRITE!*, ’(A)*)’  TAKING  ANOTHER  CUT* 

Q  =  Q-QDECR 
NCLAST  =  NCOUNT 
GO  TO  179 
ENDIF 

IFfNCHI.GE.2)THEN 

RCHI  =  CHIOPD(NCHI- 1  )/CHIOPD(NCHI) 

IF(RCHI.GT.RCUT)THEN 
WRrrE(*,'(A)')’  TAKING  ANOTHER  CUT* 

Q=Q-QDECR 
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NCLAST  =  NCOL'NT 
GO  TO  179 
ELSE 

WRITEf.’tA))’  NO  FURTHER  IMPROVEMENT’ 

I  FfCHIOP .  G  E .  C  UTM  AXITHEN 

WRITE!*.  (A)T  CHI  SQUARE  STILL  TOO  HIGH.  TRAJECTORY  REJECTED’ 

NEV  =  NEV-1 

ENDIF 

ENDIF 

ENDIF 

ENDIF 

ENDIF 

WRITE!*. '(A)')'  LEAVING  OMEGA’ 

NEV  =  NEV  +•  I 
NCOUNT  =  0 
103  ICHI  =  0 

LLST  =  LLLST 
NCHI  =  0 

IF!LLST.LT.(NWORDS-2*NSM))GO  TO  179 
Q  =  QDEF 

1 1  CONTINUE 

C 

C  END  OF  OUTER  LOOP 

C 

107  CLOSEI6) 

CLOSEI7) 

CLOSE(8) 

CLOSED) 

C 

C  WRITS  DATA  TO  DISK 

C 

WRITE(*.’(A\)’)'  ENTER  FILENAME;  ’ 

READ  (*,’(A)’)  ZNAME 

OPENHDFILE.FILE  =  ZNAME. STATUS  =  NEW’.  ACCESS  =  SEQUENTIAL’ . 

1  FORM  =  FORMATTED’) 

DO  109  J=  I, NEV-1 

WRrTE(IDFILE.193)TM(J),WZ(J),WX(J),WY(J) 

WRITE(IDFILE.200)DWZ(J),DWX(J).DWY(J),CHI2(J),NU!J) 

109  CONTINUE 

CLOSE(IDFILE) 

C 

193  F0RMAT(4<1PE12.5)) 

200  F0RMAT(12X,4(1PE12.S),I4) 

211  FORMAT!’  DEFAULT  FIT  IMPROV.  RATIO;’,F12.4) 

212  FORMAT!’  DEFAULT  EDGE  DECREMENT;’, F  12.6) 

2121  FORMAT!’  DEFAULT  TRAJ.  ENDCLIP  FRACTION;’, F12.6) 

219  FORMAT!’  DEFAULT  EDGE  CUTOFF;’,F12.5) 

237  STOP 
END 
C 
C 

c 

PLOTTING  SUBROUTINE 
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<J  u  u> 


CALLS  TO  PLOT88  PLOTTING  ROUTINES-NOT  INCLUDED  HERE 


C 

C 

SUBROUTINE  CAT(JDFILE) 

C 

C  THIS  SUBROUTINE  READS  LECROY  6810  SINGLE  CHANNEL  TIME  DATA 

C 

c 

C 

SUBROUTINE  BASE(NCHAN,OFFSET,RDATA,NEND,TIME) 

C 

C  THIS  SUBROUTINE  SUBTRACTS  THE  MOST  PROBABLE  BASELINE 

C  FROM  THE  RAW  PSD  SIGNALS 

C 
C 
C 

SUBROUTINE  INTERP(X.Z) 

C 

C  THIS  SUBROUTINE  PERFORMS  A  BILINEAR  INTERPOLATION 

C  ON  THE  NON-ORTHOGONAL  DETECTOR  MAP  GRID  TO  FIND  TRUE 

C  REFLECTION  COORDINATES 

C 

C 

C 

C 

C 

C 

C  OMEGA3  VERS.  1  DEC  87 

C 

C 

C  COMPUTES  THETA. PHI  AND  ASSOCIATED  COEFFICIENTS  FROM  X.Z  PAIRS 
C  COMPUTES  dTHETA7dt,dPHI7dt  ....ACCUMULATES  SUMS.  AND  SOLVES 
C  3X3  TO  EXTRACT  Wx.  Wy.  Wz. 

C 

SUBROUTINE  OMEGA(ISEG.IAUTO.ITURN.ICHI.CSUM) 

REALM  WGHTX(1024),WGHTZ(1024) 

REALM  SD{3),W(3) 

REAL*8  SA(5,5),DET,SS1 1,SS22,SS33,SS12,SS13,SS23 

INTEGERS  SQRFLG 

COMMON/FLG/SQRFLG 

COMMON/OPDT/PHP(  1024), THP(  1024), PFUNC(  1024), TFUNC(  1024) 

COMMON  AVDATAVD(  1024),  YD(  1024), SINX(1024),SINY(  1024), NCOUNT 
COMMON/CD  AT/PERIOD, SCALE,  Y,SNOISE,BTX,BTY,CDEFHI,CDEFLO,CUTMAX.Q 
COMMON/OMEG/TM(1024),  WX(  1024),  WY(1024),  WZ(  1024), NU(  1024) 

COMMON/ST  AT/D  WX(1024),DWY(1024),DWZ(1024),CHI2(1024) 
COMMON/TIMR/TOTAL,LST,TSCA,NCSEG 
COMMON/PROX/PROXC,PROXW,PROXP,CLIPR 
PI  =  3. 14159 
C 

WRITE!*, ’(A)’)’  ENTERING  OMEGA’ 

C  LOAD  DATA 
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SY=  I  0 

iF(Y.LT.O.O)SY  =  -1  0 
DT  =  PERIOD 
SNT  =  SNOISE 
AS  =  0  0 
NP  =  3 
ISL’P  =  1 


CLEAN  UP  TRAJECTORY  ENDPOINTS 


edgl  =  SQRT(wd<  l  )*wd<  l)  -yd)  l)*vd(  1)) 

edgn=  SQRT(wd)  ncount )  *wd(  ncount )  +  yd<ncount)*yd(  ncount)) 


11  =  3 

IF(EDG  I .  LT.0. 96*0)1 1=3-  INT(NC0UNT*CL1PR) 

12  =  NCOUNT-2 

IF(EDGN.LT.0.96*Q)I2  =  NCOUNT-2-INT(NCOUNT*CLIPR) 

COMPUTE  SIGNAL-TO-NOISE  STATISTICAL  WEIGHTS 
CORRECT  OPTICAL  DISTORTION 

NMAX  =  I2-I1-1 
DO  3  1  =  I.N’MAX 
II-I  +  IM 
X  =  WD(II) 

Z=  YD(H) 

WGHTXfD  =  2*SINXni)*SINX(II)*Y*Y*DT*DT/(SNT*SNT*(  1  +  X*X)) 
WGHTZm  =  2*SINYnn*SIN  Y(in*Y*Y*DT*DT/(SNT*SNT*(  1  +  Z*Z)) 
IFHTURN.  EQ.  1  )WGHTZm  =  0.0 

IF(ITURN.EQ.2)WGHTX(I)  =  0.0 

RF  =  SQRT(X*X+Z*Z) 

RN  =  1  284*ASIN(RF/1 .284) 

X  =  X*RN/RF 

z=z*rn/rf 

PHP(I)  =  S  Y*PI  '2.0- ATAN(X/Y) 

THPfD  =  PI/2.0-ATAN(Z/SQRT(X*X  +  Y*Y)) 
WRITE(*.399)THP(I).PHP(I) 

N  =  I 

CONTINUE 

CONTINUE 

INVOKE  ORTHOGONAL  POLYNOMIAL  DECOMPOSITION 

CALL  OPD(N.PHP.PFUNC.AP.BP.CP.DPH.BTX.CHP) 

CALL  OPD(N.THP.TFUNC,AT,BT.CT.DTH.BTY.CHT) 

CSUM  =  CHP  +  CHT 

SNR=0.5*(BTX  +BTY)/SNOISE 

NOPT  =  2*INT(0. 2/SQRT((BT*BT  +  BP*BP)*SNR)) 

IF(NOPT.LT.2)NOPT  =  2 

IF(NOPT.GT.N/2)THEN 
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WRITE!*. '(A)')’  WARNING:  NON-OPTIMUM  DIFFERENCING  NEEDED’ 

NOPT  =  N<2 

ENDIF 

C 

DO  33  !=  I.N-1 

33  AS  =  AS  +  SQRT!(PFUNC!I  +  1  )-PFUNC(I))**2  +  <TFUNC(I  i-  1  )-TFUNC(I))**2) 

C 

IFMAUTO.NE.  DTHEN 
C 

WRITE(*.’(A\)T  DISPLAY  TRAJECTORIES?  (1  =  YES);’ 

R£AD(*.'(I1)’)ITD 

IFOTD.EQ.DCALL  FTDISP(PHP,4.N) 

C 

WRJTE(*.502)AP.BP,CP,DPH,BTX,CHP 
WRITE!  •.  503)  AT,  BT.  CT.  DTH.BTY ,  CHT 
WRITE(*.507)NOPT,N 
C 

CCUT  =  CDEFHI 
CCUTLO  =  CDEFLO 

36  WRITE!*.  506 ICDEFHI.CDEFLO.CUTM  AX 

WRITE(*.’(A\)’)'  SET  OPD  CHI2S?  (I  =  YES);’ 

READ!*,'(II)')ITD 
IF(ITD.EQ.  DTHEN 

WRITE!*. ■(A\)’)•  ENTER  MAXIMUM  ACCEPTABLE  OPD  CHI2;’ 
READ!*.*.ERR  =  36)CDEFHI 

WRITE!*. '(A\)')'  ENTER  MINIMUM  SINGLE  PASS  VALUE;’ 

READ(*.*.ERR  =  36)CDEFLO 

WRITE!*. ’(A))’)-  ENTER  MAXIMUM  ACCEPTABLE  CHI2  AFTER  ALL  CUTS;’ 
READ!*.  *.  ERR  =  36)CUTM  AX 
CCUT  =  CDEFHI 
CCUTLO  =  CDEFLO 
ENDIF 
ENDIF 
C 

IF(CSUM.GT.CCUT)THEN 

WRITE!*. ’(A)')’  OPD  TRAJECTORY  REJECTION!’ 

ISEG  =  ISEG-1 

1CHI=0 

RETURN 

ELSEIF!(CSUM. GT.  CCUTLO).  AND.  (CSUM.LT.CCUT))THEN 
ICHI=  1 
ELSE 
ICHI=0 
ENDIF 
C 
C 

C  INITIALIZE  ARRAYS 
C 

S2=0.0 
DO  6  K=  1.3 
SD(K)=0.0 
W(K)=0.0 
DO  5  KK  =  K,3 
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SA(K.KK)=O.DO 
5  CONTINUE 

ft  CONTINUE 

C 

WRITE!*. ‘(A)')’  ACCUMULATING  SUMS' 

C 

C  COMPUTE  DERIVATIVES.  COEFFICIENTS.  AND  ACCUMULATE  SUMS 

C 

D07L  =  I+  NOPT/2.N-NOPT/2 

DTHPDT  =  (THP(L  r  NOPT/2)-THP(L-NOPT/2))/(NOPT*DT) 

DPHPDT  =  (PHPfL-NOPT/2)-PHP(L-NOPT/2))/(NOPT*DT) 

C 

SQRFLG=0 

CALL  TRANS(THP(L).PHP(L),SY,A,B,C,D,E,F) 

IFISQRFLG.EQ.  I  (RETURN 

C 

WGHTX(L)  =  NOPT*NOPT*WGHTX(L)/4.0 
WGHTZ(L)  =  NOPT*NOPT*WGHTZ(L)/4.0 

c 

S2  =S2  +•  DTHPDT*DTHPDT*WGHTZf  L)  +  DPHPDT*DPHPDT*WGHTX(L) 

C 

SD(  1 )  =  SD(  1 )  +  DTHPDT*A*WGHTZ(L)  +  DPHPDT*D*WGHTX(L) 

SD(2)  =  SD(2)  +  DTHPDT*B*WGHTZ(L)  +  DPHPDT*E*WGHTX(L) 

SD(3)  =  SD(3)  DTHPDT*C*WGHTZ(L)  +  DPHPDT*F*WGHTX(L) 

C 

SA(l,l)  =  SA(I,l)*  A*A*WGHTZ(L)  +  D*D*WGHTX(L) 

S  A(  1 .2)  =  S  A(  1 .2)  +  A*B*WGHTZ(L)  +  D*E*WGHTX(L) 

S  A(  1 .3)  =»  S  A(1 .3)  +  A*C*WGHTZ(L)  +  D*F*WGHTX(L) 

SAI2.2) =SA(2.2)  +  B*B*WGHTZ(L)  +  £*E*WGHTX(L) 

S  A(2.3)  =  SA(2.3)  +  B*C*WGHTZ(L)  +  E*F*WGHTX(L) 

SA(3,3)  =  SA(3.3)  +  C*C*WGHTZ(L)  +  F*F*WGHTX(L) 

7  CONTINUE 

IF(NP.EQ.2)THEN 

SA(3.3)=0.0 

SA(I.3)=0.0 

SA(2.3)=0.0 

SD(3)=0.0 

W(3)=0.0 

ENDIF 

C 

SSll=SA(l.l) 

SS22=SA(2,2) 

SS33  =  SA(3,3) 

SS13  =  SA(I,3) 

SS23=SA(2,3) 

SS12  =  SA(1,2) 

C 

C 

C  ...  A  SYMMETRIC  MATRIX 
C 

SA(3,2)  =  SA(2,3) 

SA(3,1)  =  SA(1,3) 

SA(2,1)  =  SA(1,2) 
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DO  20  M  =  l .  3 

20  \VRITEf*.449)SD(M).SA(M.I).SA(M.2).SA(M.3) 


WR1TE(*,'(A)’)’  entering  matinv- 
SQRFLG  =0 

CALL  MATINVfSA.NP.DET) 

IFfSQRFLG.EQ.  HRETURN 
WRITEf*  •(A)’)*  LEAVING  MATINV’ 

1F(DET.LE.0.0)THEN 

WRITE(V(A)’)’  ERROR;  SINGULAR  MATRIX’ 

ISEG  =  ISEG-1 
ICHI  =  0 
RETURN 
ENDIF 

DO  30  M=  1,3 

WRrTE(*.450)SA(M.l).SA(M,2),SA(M.3) 

COMPUTE  Wx.Wy.Wz 

DO  9  1=1. NP 
DO  8  J  =  l.NP 
W(I)  =  W(I)  +  S  A(I,J)*SD(J) 

CONTINUE 

CONTINUE 

COMPUTE  CHI-SQUARE  AND  UNCERTAINTIES 

NU(ISEG)  =  2*N-7 
IFflTURN.GE.  l)NU(ISEG)=N-5 

CHID  =  W(  1  )*W{  1  )*SS  1 1  +  W(2)*W(2)*SS22  +  Wf  3)*Wf  3)*SS33 
CHIC  =  2.  D0*(W(  1)*W(3)*SS  13  +  W(  1)*W(2)*SS  12  +  W(2)*W(3)*SS23) 
CHI  Y  =  -2 .  D0*fW(  1  )*SD(  1 ) + W(2)*SD(2) + W(3)*SD(3)) 

CHI20SEG)  =  (CHID  +  CHIC  +  CHIY  +  S2)/NU(ISEG) 

IF((SS1  l.LE.O.O).OR.(SS22.LE.O.O))THEN 

WRJTE(*,’(A)’)’  ERROR;  NON  POSITIVE-DEFINITE  DIAGONAL’ 

ISEG  =  ISEG-1 

ICHI=0 

RETURN 

ENDIF 

WZfISEC)  =  2.  *W(  1 ) 

WX(ISEG)  =  2.*W(2) 

WY(ISEG)  =  2.*W(3) 

IF(SA(  1, 1).LT.0)THEN 
SQRFLG =1 
RETURN 
ENDIF 

IF(SA(2,2).LT.0)THEN 
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SQRFLG  =  1 
RETURN 
ENDIF 

1F(SA(3.3).LT.0)THEN 
SQRFLG  = l 
RETURN 
ENDIF 

DWZ(ISEG)  =  2.  *SQRT(SA<  1 .  L>j 
DWX(ISEG)  =  2.  -SQRT/SAIZJ)) 

DWY(ISEG)  =  2.*SQRT(SA(3,3)) 

TM(ISEG)= TOTAL -r  PERIOD*LST 

PROXIMITY  TESTS 

AF  LAG  =0.0 
IF(NP.EQ.2)GO  TO  40 
UZ=COS(0.5*(AT  +  PI/2.0)) 

UX  =  SIN(0.5*(AT  -i-  Pl/2.0))*COSC(0.5*AP)-PI/4.0) 

UY  =  SIN(0.5*(AT  f  PI/2.0))*SIN{(0.5*AP)-PI/4.0) 

POLE  Y  =  - 1 

UDOTW=2.  *(UZ*Wf  i )  +  UX*W(2)  +  UY*W(3)) 

WM  AG  =  2.  *SQRTfW(  |)*W(|)+  W(2)*W(2) + W(3)*W(3)) 

TSTMAG  =  2.  *SQRTfW(  1  )-W(  1 )  +  W(2)*W(2)  +  AS*AS*W(3)*W(3)) 
DWMAG  =  SQRT(DWZ(ISEG)*DWZ(!SEG)  +  DWX(ISEG)*DWX(ISEG)  + 
I  AS*AS*DWY(ISEG)*DWY(ISEG)) 

WDOTP=2*W(3)/WMAG 
IF(WMAG.EQ.0)THEN 

WRTTE(*,’(A)')’  DIVIDE  BY  ZERO  ERROR  -  SKIPPING  SEGMENT’ 
ISEG  =  ISEG-1 
ICHI=0 
RETURN 
ENDIF 

THPAR  =  ABS(UDOTW/WM  AG) 

IF(THPAR.GE.  PROXW)THEN 
WRITE(V(A\)’)’  W  NEARLY  PARALLEL  TO  N!’ 

ISEG=ISEG-l 
ICHI=0 
RETURN 
ENDIF 

COSERR = TSTM  AG/SQRT  (TSTM  AG  *TSTM  AG  +  DWM  AG*DWM  AG) 
IF((ABS(COSERR).LE.ABS{THPAR)).OR. 

1  (AS.LE.ACOS(ABS(UY))))AFLAG=  1.0 
IF(AFLAG.EQ.  1.0)THEN 

WRITE(V(A\)y  W  ERROR  CONE  CONTAINS  N!’ 

DZTMPR = DWZ(ISEG) 

DXTMPR  =>  DWX(ISEG) 

DYTMPR = DWY(ISEG) 

ENDIF 
C 

40  IF(IAUTO.NE.  l)THEN 

®*WRITE(*,500)  WZ(ISEG)f  WX(ISEG),  WYfISEG) 

WRITE(*,50 1 )  DWZ(ISEG) ,D WX(ISEG), D WY (ISEG) 
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WRITER, 499)  NL’(1SEG).CHI2(ISEG) 
WRITEC.504)  UZ.UX.UY 
WRITE!*. 505)  THPAR.COSERR.WDOTP 
WRITE!*. 508)  AS 


END1F 

IF(AFLAG.EQ.  l.OlTHEN 

lF((ABS(UY).LE.PROXC).OR.(THPAR.GT.ABS(UY)))THEN 
ISEG  =  ISEG-! 

ICHI  =0 

RETURN 

ENDIF 

C  1F((ISUP.EQ.  1).  AND.!NP.EQ.3))THEN 

NP  =  2 
GO  TO  18 
C  ELSE 

C  1SEG  =  1SEG-1 

C  RETURN 

C  ENDIF 

ENDIF 

1F(N'P.EQ.2)THEN 
DWZ(ISEG)  =  DZTMPR 
DWX(ISEG)  =  DXTMPR 
DWY(ISEG)  =  DYTMPR 
C  1CHI  =  0 

ENDIF 
C 
C 

398  FORMAT(  '0  THETA1  PHI1') 

399  FORMAT!2!IPD12.4)) 

-U9  FORMAT! '0  ’.  1PD12.4,’  >  ’,3(1PD12.4)) 

450  FORMAT!  '0  MATRIX  INVERSE’, 3(  1PD12.4)) 

499  FORMAT! ’0  CHI-SQUARE  PER  \I4.’  DEGREES  OF  FREEDOM  =  ’.  1PE12. 4) 

500  FORMAT!  ’0  Wz=  IPE12.4,’  Wx  =  1PE12.4.’  Wy  = ’.  1PEI2.4) 

501  FORMAT! ’0  +  /-.1PE12.4.’  ’.1PE12.4,’  MPE12.4) 

502  FORMAT! ’0  PHI:  *  .6!  1  PE  1 2. 4)) 

503  FORMAT! ’0  THETA:’,6(IPE12. 4)) 

504  FORMAT! ’0  Uz=’.Fl2.6,'  Ux  =  ’.F12.6,’  Uy  =  ',F12.6) 

505  FORMAT! ’0  W;N  COSINE’,F12.6,’  W;DW  \F12.6.’  W:C  ’.F12.6) 

506  FORMAT!  'DFLTS;  MAX’,F12.6,’  MIN',F12.6,’  MAX  AFTER  CUTS  \F12.6) 

507  FORMAT! ’0  OPTIMAL  SAMLING  NOPT= 14,' OVER ’,14,’ PTS’) 

508  FORMAT!  ’0  TRAJ.  ANGULAR  SIZE’,F12.6) 

RETURN 

END 

C 

C 

SUBROUTINE  TRANS(THP.PHP,SY.A,B.C,D,E,F) 

C 

C  TRANSFORM  TO  UNPRIMED  ANGLES,  COMPUTE  COEFFTS. 

C 

C 

C 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


SUBROUTINE  MATINV( ARR.N.DET) 

MATRIX  INVERSION  ROUTINE 

SUBROUTINE  OPD<N.Y.YFUNC,A.B,C,D,BT,CHI2) 

ORTHOGONAL  POLYNOMIAL  DECOMPOSITION  SUBROUTINE 

END 
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Appendix  B 


SNR  Distribution 


The  signai-u>-noise  ratio  (SNR)  is  determined  by  the  ratio  of  the  total  incident  power  from  a  valid  retleci 
and  total  equivalent  noise  power  Pvy  The  noise  power  which  has  a  variety  of  unconelaied  sources, 
predonunantlv  fluctuating  background  light  Pg  and  thermal  detector  noise  Py,  is 


it  t  (1/2 

PN  =  IPS _ ^  Pi 


B,rms  T.rms, 


rT,rms 


NEPf 


1/2 


uhere  NEP  is  the  noise-equivalent  power  factor  of  the  detector  and  f  is  the  amplifier  band-width.  (Dete 
noise  includes  thermal  noise  and  shot  noise.)  The  noise  voltage  rms,  is 

I  2  2  .2  12 

^N.rms  "  I^B.rms  *  ^T.rms  *  ^A.rms/ 


with 


'  B.rms 


=  P, 


B.rms 


VT.rms  =  NEP 


fl'2 


V  =  V 

A. rms  A, rms 


where  Sq  is  the  detector  sensitivity  and  Ry  is  the  pre-amplifier  transimpedance  and  VA  rms  is  the  pre-ai 
voltage  noise.  The  signal  voltage  is  Vs  =  Pr  sD  RT- 


Both  the  background  fluctuation  noise  and  the  pre-amp  noise  are  typically  small  compared  to  detector  nc 
the  SNR  is  computed  on  the  basis  of  this  value.  (In  the  unusual  case  where  background  fluctuations  are 

errors  are  forced  into  goodness -of-fit  parameter.)  For  NEP  =  IO’^W/^hF,  Sjj  =  0.5  amps/W.  . 

Ry  =  10^  V  amps,  and  the  band-width  f  -  50  kHz.  the  detector  noise  V^,  rmsV,  is  approximately  1  n 

The  SNR  is 
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SNR  =  _ — _ 
^  N.rms 


Ps 

PN.rms 


Since  V^t  is  constant,  the  instantaneous  SNR  is  always  proportional  to  Vv  The  distnbution  of  SNR  over  a  data 
■-et  is  thus  equivalent  to  the  distribution  of  reflected  laser  power. 


The  distnbution  of  reflected  power  is  itself  composed  of  vanations  due  to  random  mirror  positions  within  the 
Gaussian  profile  of  the  incident  beam,  and  the  vam.tion  of  mirror  diameter,  d.  The  instantaneous  reflected 
power  is  Pf 


Pr(r,d)  =  Io-^|5 
V2  l  2 


(B-l) 


where  i  r  -  r .  j  is  the  distance  from  the  ( ax  i -symmetric)  beam  center  and  Og  is  the  beamwidth.  The  factor  of 

J_  is  included  to  account  tor  the  near  45  deg  angle  of  inclination  of  the  mirror-normal  to  the  incident  beam. 

vT 

2 


The  area  of  the  hexagonal  lead  carbonate  mirror  is  approximated  by  ir 


.  In  the  special  case  where  the 


(low  velocity  vector  is  parallel  to  the  incident  beam,  then  Pr(r,d)  is  constant  and  given  by  Eq.  (B-l).  When  the 
incident  beam  is  inclined  to  the  (low.  the  reflection  intensity,  and  thus  the  SNR.  vanes  with  time  and  position 


Pr(y.t.d)  =  I0 


v^2 


(utsinX)-  y 
2  2<4  2a: 


B 


(B-2) 


where  u  is  the  flow  velocity,  X  is  the  angle  between  the  velocity  vector  and  the  incident  beam  direction,  and  y 
is  the  coordinate  perpendicular  to  both  velocity  and  incident  beam.  This  is  illustrated  in  Figure  B-l.  Only  two 
cases  are  considered,  sinX  =  0  and  1,  i.e.,  parallel  and  perpendicular  beam  and  velocity. 
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Definition  Sketch.  Incident  beam  inclined  to  the  flow  gives  nse  to  time  varying  SNR 


A )  s;nA  =  0 

The  incident  light  intensity  in  the  sampled  volume  is 


Hr)  =  IQe 


r* 


But.  p(i)dl  =  p(r)dr.  and  p(r)  =  2tt.  so 


<VD 


eliminating  r  using  Eq.  (8-3) 


PI(I)  =  £  for  0  <1  <10 

A  threshold  will  determine  the  actual  lower  intensity  limit,  and  thus  the  normalization  c. 


Mirror  diameters  are  approximately  log-normally  distributed  with  most  probable  value  dQ  and  log-normal 
width 


(Ind 


•nd. 


p(d)  oc  e 


2(ln<rd)- 
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i  For  ZTX-B  dQ  -  15  ^m,  Oj  *  4<3.)  The  PDF  for  mirror  area  A  is  therefore 


pa(a>  a 


IlnA  -  ioAo)“ 
1  2(lnffd)- 

Ae 


The  reflected  power  Pr  =  IA.  thus  ln(Pr)  =  ln(I)  +  ln(A)  and 


lnl0 

p<lnPr)  =  |  pj(lnI)p^(lnPr  -  InJ)d(lnl) 


Integrated  and  transforming  ln(Pr)  -•  Pr. 


p(SNR)  =  pr(Pr)  oe  J_  erf 

*  r 


lflao^rAo) 
2\Jl  lncrd 


In 

PN 

PN  erf 

SNR  A0 

SNR 

2y/2  In  ad 

(B-4) 


The  threshold  sets  a  minimum  Pr  threshold  determines  the  normalization  constant.  Note  that  the  threshold 
interacts  with  the  beam  diameter  through  Pr  Lfueshoid  ~  ,Athreshoid  Lar£er  mirrors  see.  in  effect,  a  larger 
beam  diameter. 


8)  sin  X  =  I 

The  situation  for  beam  perpendicular  to  velocity  vector  is  more  complex.  The  mirror  diameters  are  log- 
normally  distributed  as  before,  but  now  the  threshold  interacts  with  both  beam  diameter  and  transit  time, 
depending  upon  where  the  particles  traverse  the  beam  profile.  To  simplify  the  discussion,  it  is  assumed  that  od 
is  sufficiently  small  that  mirror  area  A  =  Aq  and  Pr  0  =  I0A0.  From  the  expression  for  the  intensity 
Eq.  (B-3),  the  PDF  of  transit  times  at  given  Pr  threshold 's 


P<Ttr) 


Ttr 


^r,o 

2 

1/2 

2 

u 

^.threshold 

tr 
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This  is  only  the  case  tor  trajectories  passing  close  to  the  detector  center.  Typically,  curvature-related  errors 


will  dominate. 


The  computation  of  5(d\fr/dr)  is  more  involved.  Estimates  of  5^/3 r,  and  6(d\(/ldt)  are 


d V1 

II 

D 

1 

*  [dtf|  .. 

.  *j  n 

rn  ra 

n  rn  ~  ra 

where  ra  is  the  position  of  the  trajectory  centroid,  i.e.,  the  a^1  point  divides  the  trajectory  in  half.  With 


51  (rn  "  ra) 
n=*l.n*a 


but  rn  -  ra  can  be  approximated  by  (n-a)rAt  with  rjj  =  Nrdt 

%  *  "“8e  N) 

thus 

f' ^  ' 
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Appendix  D 


Proximuv  Tests 

The  geometry  "I  traiectones  in  the  detector  plane  shown  in  Figure  17  forms  the  basis  of  another  set  of  : 
particular,  the  vectors  a  y.  and  c.  are  monitored  in  order  to  test  for  conditions  leading  to  erroneous 
interpretations  ot  data  near  the  singularity  in  the  VOP  geometry  ( u  parallel  to  n).  Since  these  tests  use 
magnitudes  i  a: .  yi ,  and  |  cl  ,  they  are  called  proximity  tests.  There  are  three  basic  types: 

1 1  Proximity  to  zero  vorticitv  magnitude. 

For  a  stationary  rejection  on  the  detector  surface  i  yi  -»  0.  A  finite  SNR  invariably  leads  to  the 
misinterpretation  ot  fluctuations  or  dnft  tn  its  position  as  a  rapid  rotation  about  an  axis  nearly  parallel  t< 
mirror-normal  vector.  It  is  highly  improbable,  however,  that  this  rotation  is  real.  Indeed,  this  would  b 
interpreted  as  zero  vorticity  except  for  the  fact  that  the  uncertainty  in  the  vorticity  diverge  near  the  VOI 
singularity. 

6(cijft)  -»  oo  as  oifl  -»  ch 

One  must  interpret  zero  vorticity  data  carefully  since  it  cannot  be  distinguished  from  contaminant  or  but 
Mgnals. 

7)  Test  for  u  parallel  to  n  within  experimental  uncertainty. 

Erroneous  curvature  measurements  for  short,  faster  moving  trajectory  fragments  tend  to  give  |  y|  <  1 
uncertainties  of  the  w  components  can  be  thought  of  as  forming  an  ellipsoidal  volume  with  semi-axes  6u 
6<uv  at  the  tip  of  each  vorticity  vector.  The  projection  of  the  ellipsoid  onto  the  plane  perpendicular  to  t! 
vorticity  vector  defines  a  conical  region  describing  the  range  of  possible  orientations  of  the  measured  v< 
vector  within  experimental  error.  If  the  mirror-normal  vector  falls  within  this  cone,  it  is  likely  that  ern 
curvature  dominates  the  vorticity  measurement.  (This  is  actually  a  special  case  of  type  (1).)  A  decisior 
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then  be  made.  When  i  ai  is  greater  than  a  preset  cutoff,  the  event  is  discarded.  When  |  ai  is  less  than  the 
cutoff,  the  contamination  of  <jj'x  and  u'z  due  to  ui'y  is  within  acceptable  limits  m  which  case  only  two 
components  are  fit  and  there  is  no  measurement  of  oi'y. 

3)  Test  for  i  cl  ■<  1  . 

This  indicates  a  measured  voriicity  consisting  of  almost  pure  Uy.  If  this  occurs,  having  passed  the  other  tests, 
the  trajectory  is  subjected  to  the  same  criteria  as  those  failing  test  type  (2). 
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Appendix  E 


Trajectory  Length  Distribution 

In  the  following,  2f  ,  represents  the  total  angular  size  of  a  trajectory  on  the  unit  sphere  limited  only  by 
time  through  the  sampled  volume,  i.e..  it  is  not  restncted  to  that  piece  which  would  fall  on  the  detector 
i  -  D  2)  is  the  masked  detector  radius,  and  the  random  variable  dfk)(  =  2rf k))  is  the  fractional  chord  lent' 
trajectory  section  on  the  detector  surface,  whose  distribution  is  sought,  f .  R.  and  r  are  used  below  for 
convenience  (without  subscripts). 

The  probability  ot  finding  a  chord  of  (half)  length  r  on  a  circle  of  radius  R  for  a  randomly  placed  segm 
length  « .  in  the  plane  of  the  circle  is  naturally  divided  into  two  regimes:  t  >  R,  and  u^R. 

1)  f  >  R 

This  case  is  itself  divided  into  two  parts:  f_:ll  chords  (spanning  the  circle),  and  partial  chords  (a  segmer 
endpoint  falling  within  the  circle).  The  situation  is  illustrated  in  Figure  E-l. 


y 


FULL 

CHORD 


Figure  E-l 

Definition  Sketch  for  Trajectory  Length  Distribution  Derivation 
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w 


The  hill  chord  distribution  is 


I 


PfU;  r  £R  <f>  =  £ 


The  partial  chord  distribution  is 


pp(r:  r  S  R  < 


dy 


-1 

R2  <-  rf 

sin 

R(r*f) 

1) 


The  total  PDF  for  case  ( 1)  is 


p(r;  r  <  R  <  1)  =  p^r)  +  pp(r)  ,  R  <  i 
which  is  shown  in  Figure  3 E-2  at  several  l/R  ratios. 


Figure  E-2 

PDFs  of  Trajectory  Length  for  1/ R  >  1 
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M  =■! 


1  -*  -> 

.hords  con  lined  to  this  same  region,  partial  chords  ot  length  *'  tor  y  <  yR"  -  i‘  .  and  partial  chords 

region.  r<  t  The  result 

i 

j  p(r;  r  <  l  <  R>  =  p,  s-  p-,  -  p3  -  p4 

is  shown  m  Figure  £-3.  The  only  point  noted  here  is  that  p^,  for  r=i.  is  a  delta-function  with  amplitud 
dependant  on  t  R  which  has  been  artificially  broadened  in  the  figure. 

t 

'Alien  t  is  itself  randomly  distributed,  that  distribution  must  be  convolved  with  the  PDFs  above  to  find  t 

i 

I 

I 

I 

i 

i 
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TRAJECTORY  LENGTH/DETECTOR  DIAMETER 


Figure  E-3 

PDFs  of  Trajectory  Lengths  for  <7R  <>  I 
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